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Abstract. The proton is a composite object with spin one-half, understood 
to contain highly relativistic spin one-half quarks exchanging spin-one gluons, 
each possibly with significant orbital angular momenta. While their fundamental 
interactions are well described by Quantum ChromoDynamics (QCD), our standard 
theory of the strong interaction, nonperturbative calculations of the internal structure 
of the proton based directly on QCD are beginning to provide reliable results. Most 
of our present knowledge of the structure of the proton is based on experimental 
measurements interpreted within the rich framework of QCD. An area presently 
attracting intense interest, both experimental and theoretical, is the relationship 
between the spin of the proton and the spins and orbital angular momenta of its 
constituents. While remarkable progress has been made, especially in the last decade, 
the discovery and investigation of new concepts have revealed that much more remains 
to be learned. This progress is reviewed and an outlook for the future is offered. 
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1. Introduction 

The nucleoli (proton or neutron) is a suitable object for our curiosity in view of the fact 
that it accounts for most of the directly visible mass of the universe. Half a century 
ago, elastic electron-proton scattering provided direct experimental evidence [1] that 
the proton has a finite size of about lfm. (With the benefit of hindsight, we now 
realize that its already-known large anomalous magnetic moment could already have 
been interpreted to indicate that it is not a point-like object [2].) The amazing variety 
of baryons and mesons discovered using new multi-GeV proton accelerators in the 1950's 
and 1960's demanded some explanation in terms of an underlying structure that would 
unify their observed quantum numbers and decay modes. Such a composition in terms 
of more elementary fractionally charged spin-| particles called 'quarks' was postulated 
in 1964 [3,4]. Then the internal structure of the proton came under direct experimental 
observation in the late 1960's at the SLAC laboratory, where an intense electron beam 
became available with an energy large compared to the relevant scale of ~ 1 GeV. In 
these now-famous Deep-Inelastic Scattering (DIS) experiments, a substantial amount 
of energy and momentum is transferred to the proton target, causing it to disintegrate. 
In 1969, it was hypothesized that the energy and momentum was being transferred to 
internal constituents of the proton, called 'partons' [5]. Information about the properties 
of the constituents can be inferred from the distribution of the scattering cross section 
in these measured transfers. In the same exciting year, it was realized that the data 
showed a property called 'scaling' [6]: at large enough momentum transfer, the kinematic 
dependence of the cross section resembles that for elastic scattering off point-like charged 
particles that might be identified with those partons. As the properties of the cross 
section also indicated that the partons have spin |, it became accepted that the charged 
partons could be assumed to be the hypothesized quarks. However, it emerged that the 
total momentum carried by observed quarks accounts for only about half of the proton's 
momentum. This suggested that there must also be electrically neutral partons, possibly 
being exchanged by quarks to provide their binding force, leading to the name 'gluons'. 
Then in the early 1970's, an elegant field theory was devised to explain this binding, 
as well as the observed weakening of this force at short distances or large momentum 
transfers, known as 'asymptotic freedom' [7,8]. This theory became known as Quantum 
Chromodynamics (QCD), the contemporary standard theory of strong interactions. 

The QCD Lagrangian can be written in one line. Nevertheless, as in other fields 
of science, simple rules can give rise to complex behaviour. The unique richness of this 
field theory accounts for the challenge presented by the study of hadron structure. The 
QCD vacuum is believed to resemble a 'dual superconductor', which, unlike conventional 
QED superconductors, repells color electric fields rather than magnetic fields. Unlike 
in QED, where vacuum polarization screens charges as their separation increases, the 
QCD vacuum provides anti-screening that results in a rapidly increasing ('running') 
QCD coupling constant a s at larger distance scales or smaller momentum transfer. 
Moreover, the analogue of the dual Meissner effect in a superconductor leads in the 
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QCD vacuum to the formation of a color flux tube connecting quarks and antiquarks in 
the same way as the magnetic flux connecting two magnetic monopoles in a conventional 
(QED) superconductor would be squeezed into narrow flux tubes. This apparently leads 
to the phenomenon of confinement, as free quarks are not observed, although no proof 
has been found for the theorem that confinement inevitably follows from QCD. In any 
case, the quantitative prediction of hadron properties in QCD is presently subject to 
uncontrolled approximations. Laborious simulations on a Euclidian lattice are now 
producing such predictions for an increasing number of observables. In contrast, the 
asymptotic freedom of QCD at short distances allows precise perturbative treatment 
of hard interactions. Hence such hard interactions with bound partons can be used to 
experimentally probe the structure of hadrons. A good representative for the energies 
that are required to reach the perturbative regime is provided by the scale parameter 
Aqcd- While this parameter appears also in the description of the running of the QCD 
coupling, it provides a rough estimate about the scale of many nonperturbative effects, 
such as the energy density needed to break the superconducting phase. The numerical 
value of Aqcd depends somewhat on the renormalization procedure (scheme) that is 
used in its definition, but a value Aqcd ~ 200-300 MeV provides a good scale for 
estimates of the relevance of effects governed by the nontrivial QCD vacuum structure. 

Quark kinetic energy together with QCD field energy accounts for most of the mass 
of the nucleon, which in turn accounts for most of the mass of the visible universe. Hence 
the nucleon is a worthy subject of investigation. Deep- Inelastic Scattering of leptons has 
continued to provide most of our present knowledge of nucleon structure. As shown in 




Figure 1. Simplified schematic diagram of semi-inclusive deep inelastic scattering. 
The large shaded arrows represent typical helicities, as described later in the text. 

figure 1, a virtual photon 7* emitted by the incident lepton e and absorbed by a bound 
quark of flavour up (u) or down (d) in the nucleon N has space-like four-momentum 
q that can be inferred in inclusive measurements in which only the scattered lepton e' 
is detected. The motion of the mostly low-mass u and d flavours of quarks populating 
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the proton is highly relativistic. Hence virtual quark- ant iquark pairs are abundant and 
may be manifest in suitable processes. (In the general context of DIS, the word 'quarks' 
can be taken to also include antiquarks.) The relativistic motion also implies that DIS 
results are best interpreted in a reference frame where the target proton has 'infinite' 
momentum opposite to that of the incident virtual photon. One may intuitively think of 
this tactic as freezing the transverse motion of the partons during the interaction time 
through time dilatation (while remembering that their intrinsic transverse momentum 
Pt is of course invariant under such a boost). A virtual photon has adequate spatial 
resolution to resolve the target's internal structure [9] if its virtuality Q 2 = —q 2 is 
sufficiently large, a requirement empirically found to be at least Q 2 > 1 GeV 2 . Given 
this condition, a key experimental observable is the distribution of cross sections in the 
Bjorken scaling variable x B = Q 2 / (2P ■ q). In the infinite-momentum frame, x B can be 
interpreted as the fraction x of the proton's four-momentum P that was carried by the 
quark before it was struck by the virtual photon. Thus we distinguish here between the 
experimentally obervable kinematic variable x B and the theoretical concept of z as a 
parton momentum fraction. Only in the parton-model approximation to the DIS process 
can they be equated. In this framework, the dependence of the DIS cross section on 
x B provides information about the x-distribution of partons in the target (often called 
Parton Distribution Functions or PDFs). For a hypothetical hadron consisting of three 
non-relativistic quarks with the same large mass, this distribution would be a <5-function 
at x — |. For systems of light quarks, this is far from being the case. 

PDFs can be considered to be intrinsic properties of hadrons only if they are 
'universal', i.e. if they can be applied in descriptions of various hard processes involving 
the same hadron. Asymptotic freedom of QCD has provided the means to overcome 
its complicating aspects to prove profound 'factorization theorems' showing that cross 
sections for several experimentally important hard processes can be represented as a 
product of a cross section for one or more perturbatively calculable partonic subprocesses 
and one or more non-perturbative but universal parton distributions. For the case 
of a DIS cross section o~(x B ,Q 2 ) in the Bjorken limit Q 2 — > oo at fixed x B , the 
quark distribution qf(x B ,Q 2 ) appears in a sum over quark and anti-quark flavours 
/ = u,d,s,u,d,s . . . weighted by the squares of their electric charges e/ in units of 
the elementary charge: 



A key feature here is the incoherence between the amplitudes for the various quark 
flavours. At least in the Bjorken limit, this may be understood to be associated with the 
possible distinguishability of hadronic final states arising from different quark flavours, so 
that interference can't arise. Parton distributions such as qf(x,Q 2 ) have probabilistic 
interpretations as parton number densities. At finite values of Q 2 , the relationship 
between PDFs and observable cross sections is only a fairly crude approximation at 
leading order in a s as in (1), so the perturbative expansion is continued to Next-to- 
Leading-Order (NLO) or even further. Based on the extensive available data set for DIS 




(1) 



/ 
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as well as other processes, there is now available a precise (except at very small or large 
x) parameterization of x distributions of the u and d 'valence' quarks that determine 
the quantum numbers of the proton, and less precisely of the virtual 'sea' quarks that 
can be produced by gluon splitting, such as anti-up (u), d, and strange (s, s). 

The dependence of parton distributions on the 'hard scale' Q 2 of the process is weak 
(logarithmic), and can be understood to arise from the finer spatial resolution of a probe 
at a harder scale, which allows it to resolve more of the virtual partons that 'dress' each 
quark. These virtual qq or gluon pairs are abundant because their small or zero mass 
makes them energetically inexpensive. With quantum fluctuations so important at short 
distances, the 'appearance' of the nucleon depends on how it is probed. This might be 
considered to be analogous to the familiar change in appearance of a macroscopic object 
when viewed in light of different colours or wavelengths, due in part to the interactions 
between the atoms of its surface. Unlike the x dependence, the Q 2 dependence of 
parton distributions is perturbatively calculable in QCD at large enough Q 2 . Thus a 
distribution that is known over the range xq < x < 1 (or alternatively < x < xq ) for 
one value of Q can be 'evolved' to the distribution over the same x-range at a larger 
(smaller) value of Q 2 , using the so-called DGLAP integro-differential equations [10-12]. 
The evolution kernels for several parton distributions have been computed up to Next- 
to-Next-to-Leading-Order (NNLO) in ct s [13-16], meaning that participation of gluons 
in hard QCD interactions are taken into account. Gluons may be created in these hard 
interactions, or they may be found in the target nucleon. Hence, at order a s or higher, 
the DGLAP equations couple the distributions for quarks and gluons. This feature has 
been exploited in global fits of all available inclusive DIS data, which cover a range in 
Q 2 from 1 to 50GeV 2 , to extract information about the gluon distribution G(x). From 
this and other types of data, G(x) is now quite well determined except at small x with 
moderately small Q 2 . For convenience in the rest of this paper, we suppress the Q 2 
dependence of distributions in our notation. 

Parton distributions depend on the order of perturbation theory relating them to 
experimental observables. Beyond leading order in a s , they lose their probabilistic 
interpretation as number densities, and can even become negative. However, at 
each order, they remain universal among various observable processes in which they 
participate, and hence are still much more than just parameterizations. Furthermore, 
their basic construction is motivated by simple semiclassical parton-model ideas, but 
implemented to account for the intricacies of a quantum system. 

For values of Q 2 below ~ lOGeV 2 , there may be additional contributions to the 
cross section that scale as 1/Q n ,n > 0. These may be corrections arising because 
the mass of the target nucleon is not negligible compared to Q 2 , or they may be due 
to 'higher- twist' subprocesses, in which an additional parton in the target nucleon 
experiences a hard interaction. (Some authors apply the term 'higher twist' for all 
'power-suppressed' terms, even those suppressed by factors of 1/Q n arising only because 
of kinematic relationships. In this case, contributions involving an additional parton are 
designated as 'interaction dependent' or 'dynamical higher twist'.) For DIS processes, 
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leading twist is twist-2, with no such contributions, 

Baryons are the fermionic bound states of the spin-| quarks. All known baryons 
have quantum numbers consistent with those of a set of three quarks of particular 
flavours. (Bound states of 5 or more quarks are theoretically possible, and in fact such 
wave function components may contribute to the ground states of ordinary baryons such 
as the proton.) Soon after the conception of quarks as the constituents of the nucleon, 
nonrelativistic models were constructed in which the mass the the nucleon was directly 
attributed to three massive quarks. One of the successes of these models was considered 
to be the remarkable consistency of the measured magnetic moments of many baryons 
with a combination of Dirac magnetic moments of their three 'constituent quarks'. 

Then in the 1970's, the internal spin structure of the proton came under 
investigation via DIS of helicity-polarized lepton beams on hydrogenous targets 
containing protons polarized with the beam quantization axis, i.e. longitudinal 

polarization. The exchanged virtual photon inherits the helicity of the beam lepton, to a 
degree that depends on the lepton kinematics. Due to conservation of helicities in hard 
or short-distance interactions, only quarks with the opposite spin direction [i.e. same 
helicity) as the spin-1 photon can absorb it in the leading-order process 7*g — > q (see 
figure 2). Hence the 'polarized virtual-photon beam' selects quarks of one helicity. From 




Figure 2. Visualization of scattering of longitudinally polarized leptons and protons. 

the difference in cross sections with same or opposite polarizations of beam and target, 
it is therefore possible to extract the number densities of quarks having the same {qf') 
or opposite (qf 1 ) helicity as the target proton in the infinite momentum frame. Here the 
symbol — designates the helicity of the quark (target proton). The quark helicity 
distribution Agj(x) is then defined as Aqf(x) = qj*(x) — qT(x), while the polarization- 
averaged distribution qf(x) introduced above can be written qf(x) = q"f{x) + q^(x). 
From these definitions, the obvious 'positivity limit' arises: \Aqf(x)\ < qf(x) The 
distribution qj(x) [Aqf(x)] is sometimes written as fi(x) [gf(x)}, although this is 
arguably confusing, because the spin structure functions fi and g x are defined to be 
functions of the observable kinematic quantity x B , while the parton distributions are 
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defined to be functions of the parton momentum fraction x, with the relations between 
them being subject to approximations such as single-photon exchange and the quark- 
parton model. For convenience hereafter, we use the term '(un)polarized' to describe 
measurements, cross sections or parton distributions, meaning (polarization averaged) 
polarization dependent. Polarized PDFs evolve logarithmically with Q 2 according to a 
set of DGLAP equations that differs somewhat from those for the unpolarized case. 

Certain x-moments of parton distributions (integrals over x weighted with x n ) have 
important physical interpretations. First (n = 0) moments can be conveniently written 
as e.g. Aqy = J dx Aqf(x). A quantity that has played a central role in the discussion of 
nucleon spin structure in recent decades is the flavour sum AE = Ylf^Qfi as ^ can be 
interpreted as the net contribution of the quark helicities to the spin of the nucleon. (At 
the present precision of the field, the proton and neutron can be considered to be related 
by isospin symmetry.) Note that contributions of individual quark flavours may tend to 
cancel one another in this sum. A pivotal event in this field in the late 1980's was the 
DIS measurement with longitudinal polarization of both muon beam and hydrogenous 
target by the EMC experiment [17, 18] at the CERN laboratory, which for the first 
time extended over a wide enough range in x B that AE could be inferred, albeit under 
the assumption of SU(3) flavour symmetry (neglecting the effects of the ~ 100 MeV 
mass of the strange quark), which allows use of additional constraints from measured 
weak decay lifetimes of hyperons. The surprisingly small resulting value suggested that 
quark helicities make little net contribution to the nucleon spin, in contrast to naive 
expectations based on the baryon magnetic moments mentioned above. In language 
that with the benefit of hindsight may now seem extravagant, this finding was called 
a 'spin crisis', and inspired more than a thousand theoretical papers attempting its 
interpretation, and several new experiments with the aim of confirming the measurement 
with higher precision, and extending it to include the 'neutron' target. 

Because of the relativistic nature of bound systems of light quarks, there exists 
a third independent flavour-set of quark distributions that have a probabilistic 
interpretation as number densities. These quark transversity distributions Sqf(x) are 
needed to describe the transverse (f) polarization of quarks in a nucleon polarized 
transversely with respect to its infinite momentum: Sqf(x) Their 
probabilistic interpretation applies only in a basis of transverse spin eigenstates, not in 
the helicity basis relevant for the Aqj(x). The distribution Sqf(x) is sometimes written 
as h\(x), although this might be misleading because no spin structure function h\(x B ) 
is defined. There exists little experimental information about transversity, because the 
coupling of hard vector bosons to quarks is insensitive to the transverse polarization 
of the quark, unless the transverse polarization of the struck quark in the final state 
is 'measured'. However, some traces of this transverse polarization appear in the 
structure of the 'jet' of final-state particles produced by the energetic struck quark. 
Information about transversity has begun to emerge from semi-inclusive measurements 
with transverse target polarization, which include the detection of one or more hadrons 
in addition to the scattered lepton (see figure 1). 
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All parton distributions actually have a two-dimensional dependence on both x 
and the intrinsic quark momentum component Pt transverse to the nucleon's 'infinite' 
momentum. This two-dimensional dependence does not factorize. The three flavour-sets 
of distributions mentioned above that have probabilistic interpretations are considered 
basic in the sense that they survive integration over p?, and it is this integral that is 
signified when only the x argument appears. Five other sets of distributions have been 
identified that also have probabilistic interpretations but do not survive such integration. 
Some of these have been found to be particularly interesting, as in DIS they embody 
the effects of a 'final-state interaction' between the departing struck quark and the 
remainder of the target nucleon. The contribution of such a distribution to the cross 
section depends on polarization and momentum vectors in a manner that is called 'T- 
odd', although this does not refer to fundamental time symmetry. The study of these 
T-odd effects has provoked a reconsideration of the meaning of parton distributions and 
their universality. Experimental evidence for T-odd effects has now been observed in the 
same semi-inclusive data from which information about transversity has been extracted. 

The dependence of the elastic lepton-nucleon cross section on the momentum 
transfer is represented by the elastic 'form factor', as its Fourier transform can be 
interpreted as the spatial charge distribution in the nucleon. Up until about a decade 
ago, information about hadron structure embodied in elastic form factors and in 
PDFs were interpreted separately. It was then discovered that these concepts could 
be unified in what became known as Generalized Parton Distributions (GPDs), also 
called off-forward or skewed parton distributions. Even more profound was the finding 
that the dependence of these GPDs on three kinematic variables encodes much more 
detailed information, about e.g. parton orbital angular momentum [19] and correlations 
between transverse position and longitudinal momentum fraction of partons [20-23]. 
Of particular interest is that a second moment of a certain combination of GPDs for 
each parton flavour was found to represent the total (helicity plus orbital) contribution 
of that flavour to the spin of the nucleon [19]. This offers the first opportunity for 
experimental access to parton orbital angular momentum, as it was discovered that the 
additional information could in principle be obtained from experimental studies of hard 
exclusive processes involving additional hard interaction (s), yet 'replacing' the struck 
quark in the target nucleon to leave it intact. Most prominent among these processes 
is Deeply Virtual Compton Scattering (DVCS), the hard exclusive production of an 
energetic real photon. Within the last decade, rapid theoretical progress has been made 
in understanding this new subject, and several DVCS observables have been measured. 

Thus the evolution of our understanding of nucleon spin structure might be 
imaginatively summarized by figure 3. 
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Figure 3. Illustration of nucleon structure. Left: valence quarks only. Middle: Gluons 
and quark-antiquarks in addition. Right: Orbital angular momentum of quarks. 



2. The longitudinal spin structure of the nucleon 

2.1. Introduction 

As mentioned above, parton distributions have been investigated experimentally using 
DIS of lepton beams. In order to measure helicity distributions Aqf(x), both beam and 
target must be longitudinally polarized. The relevant observable is the difference in 
cross sections with same or opposite polarizations of beam and target: 

2 W1. (2) 



0-LL{X B/ 



a 



However, measurements of absolute cross sections are technically difficult with polarized 
targets. Hence the experimental results are typically extracted from the data as 'double- 
spin' asymmetries: 

A \\ M = — 3 

where L (U) designates Longitudinal polarization (Unpolarized) appearing in pairs 
corresponding to beam followed by target. The polarization-averaged cross section 



GUU\?i 



2 L 



a 



(4) 



has been previously measured precisely using unpolarized targets. Using these known 
values and measured values of An, o~ll can be extracted. This approach based on cross 
sections avoids confusions that sometimes arise. For example, it has been proposed 
that, if the values of o~ll are interpreted theoretically in leading order, then the 
values of o~uu should be calculated in leading order from a theoretical representation in 
terms of parton distributions, neglecting the portion involving so-called 'longitudinal' 
virtual photons! \ This misses the key point that auu became involved only because 

| Longitudinal virtual photons arc those having 'magnetic' quantum number m = with respect to 
the quantization axis along the photon momentum direction, while photons with 'magnetic' quantum 
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of the experimental convenience of asymmetry measurements, and is irrelevant to the 
theoretical intepretation of spin observables. Hence it should be treated empirically in 
terms of experimental data, with no influence from theoretical concepts such as PDFs or 
order in a s . Another point of confusion is that the extraction has often been formulated 
in a way influenced by the technical convenience that precisely measured values of o~uu 
are typically represented in terms of an unpolarized structure function called F 2 and 
the ratio R of the cross sections for longitudinal and transverse photons. Unfortunately, 
various values R were used to extract the values of F 2 , so the world data set for F2 
contains inconsistencies. This confusion can be avoided by using in equation 3 the 
original experimental values of o~uui reconstructed if necessary using the values of R 
employed in the original extraction of F 2 . 

If the measurement of A\\ is inclusive, i.e., only the scattered lepton is detected 
while ignoring the produced hadrons, an is linearly related via a kinematic factor that 
includes the QED hard scattering cross section to the spin structure function gi(x B ). 
(There is a small contribution to o~ll from another structure function g 2 , which has no 
probabilistic interpretation, and is measured using a transversely polarized target.) At 
leading order in a s , g\ for the proton (p) or neutron (n) target has a simple probabilistic 
interpretation in terms of quark helicity densities for flavour /: 



/ 

The same linear combination relates the first moments Aqj of the helicity densities to 
the moments r^ ,n of g\. 



Since the quark charges enter only as squares, inclusive DIS measurements cannot 
distinguish quarks from antiquarks, and hence can constrain only combinations of 
Au + Au, Ad + Ad and As + As. As mentioned earlier, charge symmetry is a good 
approximation, implying Au p + Au p = Ad n + Ad n and Au n + Au n = Ad p + Ad p , and 
As + As is the same in the neutron and proton. (By convention, when the superscript 
is absent, the proton is implied.) Hence measurements of and (inferred from data 
on proton and e.g. deuteron targets) provide constraints on two linear combinations of 
these three unknowns. 

A third such constraint can be derived from measurements of the lifetimes of weakly 
decaying hyperons, under the assumption of the SU(3) flavour symmetry among quark 
distributions that would apply in the limit of no mass difference among quark flavours. 
This symmetry has been found to be valid at the few percent level in a study of the 
relationship between the lifetimes of hyperons. The additional third constraint allows 

number m = ±1 are called 'transverse'. Real photons are always transverse. Sometimes the word 
'polarization' is included in this unfortunate nomenclature, even though the distinction has nothing to 
do with the longitudinal or transverse direction of a spin quantization axis. 




(5) 
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all three unknown moments to be extracted, and also their flavour- singlet combination 



As described below, it has been inferred from a series of DIS experiments of ever 
increasing precision that AS = 0.2 . . . 0.4 and As + As is significantly negative. 
This implies that there must be substantial contributions from quark orbital angular 
momentum L q and/or gluon total angular momentum Jq, in order to complete the 
'proton spin budget'. The experimental investigation of these contributions is the 
subject of following sections of this article. Subtle issues about the conceptual basis 
of spin decompositions of a system of interacting constituents are discussed in detail in 
section 5.1. 

2.2. Axial charge of the nucleon 

The axial charge g& of the nucleon is obtained from the forward matrix element of 
the operator ijj^^^ip in the same way as the familiar (vector) charge gy is obtained 
from the forward matrix element of the operator ip^^ifj. In terms of quark helicity 
eigenstates, the contribution of each quark flavour / to the axial charge of the nucleon 
involves the same linear combination qy~ — + qj~ — qj~ that also enters the polarized 
structure function gi(x B ). As a result, the integral of gi(x B ) yields the quark-charge- 
squared-weighted axial charge gA of the nucleon. Since the isovector axial charge of the 
nucleon also describes the /3-decay of the neutron, one can thus relate the integral of 
the isovector polarized structure function g\(x B ) to the neutron /3-decay constant. In 
the scaling limit of infinite Q 2 , this relationship between the spin structure functions of 
the proton and neutron is: 



While we presented here only a simple plausibility argument, this famous Bjorken 
Sum Rule (BSR) was derived [24] in the Bjorken limit of infinite Q 2 even before 
the advent of QCD, using only Gell-Mann's current algebra relations [25-27] together 
with a reasonable assumption about smooth asymptotic high-energy behaviour of the 
virtual Compton amplitude. Within the modern context of QCD where the sum rule 
can be evolved to finite Q 2 , it can be tested precisely by practical experiments. If 
a violation were discovered, the implications would be profound. This is expressed 
eloquently in [28]: "Occasionally it is worth reminding ourselves what it means to 
'understand' something in QCD. In the absence of fundamental understanding we often 
invoke 'effective descriptions' based on symmetries and low-energy expansions. While 
they can be extremely useful, we should not forget that a thorough understanding allows 
us to relate phenomena at very different distance scales to one another. In the case of 
Bjorken's sum rule, the operator product expansion, renormalization group invariance 
and isospin conservation combine to relate deep inelastic scattering at high Q 2 to the 
neutron's /9-decay axial charge measured at very low energy. Even target mass and 
higher twist corrections are relatively well understood." 



AS = (Am + Am) + (Ad + Ad) + (As + As). 



(7) 




(8) 
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2.3. Inclusive DIS experiments 

Experimental information about quark helicity densities in the nucleon has so far come 
only from DIS of polarized high-energy charged leptons on targets containing polarized 
protons, or polarized light nuclei containing neutrons. The technologies for polarizing 
beam and target are challenging, and have been the subject of intense efforts over 
several decades. A high beam polarization of 85% was achieved in even the first 
experiments [29], done at SLAC in the 1970's using a 6-23 GeV beam of electrons with 
polarization inherited in the photo-ionization of lithium atoms selected via the Stern- 
Gerlach method. However, the beam intensity was quite limited, and this approach was 
eventually replaced in modern SLAC experiments [30-35] using electron beam energies 
up to 50 GeV, polarized at a level of typically 85% in the electron source by photoelectric 
emission from a gallium arsenide surface. Beam intensities increased until they were 
limited by the heat tolerance of the cryogenic targets. Measurement times are of order 
a few weeks. 

Meanwhile, a series of even higher energy experiments at CERN [36-43] used 100- 
200 GeV beams of muons naturally polarized to a level of 80% by the weak decay of 
the precursor ir or K mesons. These tertiary beams had limited intensities of typically 
4.5 x 10 7 per spill of 2.4 s duration, with a spill period of 14.4 s, but the more massive 
muons suffer less radiative energy loss in the target materials than electrons, allowing 
much larger target thicknesses of order meters instead of centimeters. Nevertheless, 
measurement times are of order a year. 

Most targets used at SLAC and CERN are cryogenic solids of butanol, ammonia, 
or more recently 6 LiD, in a strong magnetic field, with most protons or deuterons in 
the target polarized through the process of Dynamic Nuclear Polarization (DNP). For 
DNP to proceed, a suitable material doped with a paramagnetic compound is cooled 
to temperatures less than about IK. In the high magnetic field (2-6.5 T), the unpaired 
electron of the dopant is polarized to almost 100%, while the proton or deuteron has a 
polarization of less than 1%. The process of DNP causes most of the electron polarization 
to be transferred to that of the protons by microwave irradiation of the sample. By 
this means, proton polarizations of more than 95% can be achieved. Typical deuteron 
polarizations are 40-55%; recently though, deuteron polarizations of 80% have been 
obtained. However, most of the material encountered by the beam in such solid polarized 
targets is non-hydrogenous. Hence the asymmetry of interest is strongly diluted by the 
resulting background events, by a factor of up to about 7. This constrains the statistical 
precision that can be achieved, and requires careful background corrections. 

As spin structure functions must be determined for both the proton and neutron, 
and targets of free neutrons are still infeasible, measurements are made using targets 
containing either deuterons or 3 He. The structures of these light nuclei are such that 
the spin of the neutron is strongly correlated with that of the nucleus. Protons can 
be replaced by deuterons in solid-state targets without fundamental changes in the 
technology. In contrast, the very different polarized 3 He targets consist of a low pressure 
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gas polarized via spin exchange with alkaline atoms, the polarization of which is pumped 
by an infrared laser. Such targets were used in a pair of experiments at SLAC. The 
dilution in these targets is contributed by the essentially unpolarized protons in 3 He 
and the thin target windows, and hence is smaller than in the cryogenic solids. 

A completely different experimental approach was followed at the DESY laboratory, 
where the 27.5 GeV HERA storage ring contained beams of typically 40 mA of electrons 
or positrons. Such stored high-energy electron beams may become spontaneously 
polarized via a small polarization asymmetry in the emission of synchrotron radiation by 
the beam particles as they are deflected by the magnetic fields of the ring [44]. (Several 
MW of power was radiated via synchrotron radiation, and was continuously replenished 
by a system of accelerating radio- frequency cavities.) The polarization rise time was 
of order | hour, and polarizations as large as 60% were achieved. To avoid disrupting 
the stored beam, only very thin targets were possible. These were tenuous gases of 
nuclear-polarized atoms of hydrogen or deuterium partially contained in an cylindrical 
open-ended storage cell coaxial with the beam. The polarized atoms were injected into 
the cell from an Atomic Beam Source based on Stern-Gerlach polarization filtering and 
radio-frequency transitions between atomic substates in a magnetic field. The atoms 
bounced off the cell walls of order 100 times before escaping from the ends, after which 
they were differentially pumped away by a large system of turbo-pumps. Polarized 
hydrogen target thicknesses of about 10 14 nucleons/cm 2 were achieved, and there was 
no target dilution from non-hydrogenous material. Measurement times were of order a 
year. 

2.4. Quark helicity densities from inclusive DIS data 

The data from the first polarized DIS measurements at SLAC were used to test simple 
theoretical expectations of the time, but did not cover a large enough x B range to permit 
the inference of the first moment of g\. The first opportunity for this appeared in the 
late 1980's with the EMC experiment at CERN, which employed 120-200 GeV beams of 
muons. As already mentioned, the Y\ moment extracted from this data was interpreted 
to imply a surprisingly small net contribution of quark helicities to the proton spin, 
including a negative contribution from the strange sea. The resulting excitement led to 
a series of experiments of ever-higher precision at SLAC [30,31,33-35], CERN [36-40], 
and DESY [45-47]. Figure 4 shows published data for the spin structure functions gi(x B ) 
of the proton and deuteron, while figure 5 shows the results for the neutron inferred from 
the data in figure 4, or from measurements with a 3 He target. This information is best 
displayed as x B ■ gx(x B ) versus loga; B , which conserves contributions to the area under 
the data representing the first moment IV 

A key issue to be addressed using the data sets containing measurements on both 
proton and 'neutron' targets is the BSR, introduced in section 2.2. To test this Rule, 
the first moments and T™ must be extracted from the data at some fixed value of 
Q 2 . However, measurements can never cover the complete range < x B < 1, so the 
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Figure 4. Published experimental results on x B ■ g\ and x B • gf versus x B , shown on 
separate panels, from Hermes [47] at DESY, SMC [40-42] at CERN, E143 [31] and 
E155 [33, 34] at SLAC, and COMPASS [48] at CERN. Error bars represent the sum 
in quadrature of statistical and systematic uncertainties. The Hermes data points 
shown are statistically correlated by unfolding QED radiative and detector smearing 
effects; the statistical uncertainties shown arc obtained from the diagonal elements of 
their covariance matrix, only. The E143 and E155 data points are correlated through 
QED radiative corrections. The lower panel shows the x B dependence of the mean 
values (Q 2 ) for the various experiments. The figure is taken from [47]. 
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Figure 5. Top panel: the structure function x B ■ g r { obtained from g\ and gf, from 
Hermes [47], SMC [40-42], E143 [31] and E155 [33,34]. Second panel from the top: 
x B ■ g" obtained from a 3 He target [45,49-51]. Total error bars are shown, obtained by 
combining statistical and systematic uncertainties in quadrature. The bottom panel 
shows the (Q 2 ) of each data point in the top two panels. E155 data have been averaged 
over their Q 2 bins for visibility. The figure is taken from [47] . 
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inference of I\ requires some estimate of the contributions to the moments from the 
unmeasured regions at the extremes of this range. These estimates are based on some 
assumption or ansatz motivated by reasonable expectations. Another problem is that 
these data were all produced by fixed-target experiments, resulting in a strong kinematic 
correlation between x B and the hard scale Q 2 , so that each data point corresponds to a 
different value of Q 2 , as shown in the bottom panels of figures 4 and 5. This implies that 
it is impossible to directly calculate the first moment of g±(x B ) from experimental data. 
Hence, some assumption or ansatz is typically employed to approximately 'evolve' the 
data from that experiment to a common value of Q 2 . (The resulting model uncertainty 
has not always been estimated and included with the final result.) Applying these 
approaches to the results of individual experiments has resulted in agreement with the 
BSR within the experimental uncertainty, which typically amounts to about 9% of the 
Sum [47]. 

Important information about the polarized parton distributions can be extracted 
from the measured structure functions via (5). The most efficient use of the data 
from the various experiments done at different beam energies and kinematic conditions 
is achieved by incorporating all available data in a global fit based on the evolution 
predicted by QCD in Next-to-Leading Order (NLO) in a s . This approach has the 
important advantage that it also provides information about the polarization of the 
gluons contained in the nucleon. In the order a s version of the DGLAP evolution 
equations, gluons couple to quarks via 'splitting': g — »■ qq. Hence the evolution of the gi 
structure functions with Q 2 contains information on the helicity density of the gluons. 
However, in NLO some choice must be made between various 'factorization schemes' 
for sharing a dependence on Q 2 between the parton distributions and the coefficient 
functions relating them to the structure functions. This introduces an arbitrary scheme 
dependence to the results of the fit, although results obtained in any one scheme can 
be transformed into any other scheme. The results presented here were obtained in the 
so-called MS scheme [52]. 

The fitted parameters appear in some functional form for the x dependence of each 
of several parton helicity densities, including the gluon helicity density AG(x), at some 
arbitrarily chosen 'starting' value Qq. (While this functional form is chosen with the 
benefit of some physical insight, the choice introduces some model dependence into the 
results from the fit.) At each step in the parameter search, the gi(x B ) function defined 
by the current values of the parameters is evolved using the DGLAP equations to the Q 2 
value of each experimental data point, and the discrepancy is taken as a x 2 contribution. 

Since the inclusive cross section is sensitive to only the square of the quark 
charges and not to their signs, inclusive measurements cannot distinguish quarks from 
antiquarks. Because of this as well as the limited precision and Q 2 range of the existing 
data, a useful fit is impossible without further assumptions. The global fits done up to 
now include the previously mentioned constraint based on hyperon beta decay lifetimes 
together with the assumption of SU(3) flavour symmetry, and sometimes also impose 
the BSR. The strong coupling constant a s may be taken from other experimental 
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Figure 6. Helicity densities of the strange sea As(x) + As(x) (left) and of the gluons 
(right), at Q 2 = 2.5 GcV 2 , for three different choices for the shape of the gluon densities, 
all with a similar quality of fit to the g\ data. The figure is taken from [53]. 



information, or may be also fitted as a parameter, in which case it is found to be 
consistent with previous knowledge. Some fits also impose an (arbitrary) symmetry 
condition among the sea quarks: Au = Ad = As = As. This is only a device for 
defining valence quark densities Aqf v = Aqf — Acjf, and does not affect the fitted 
strange sea quark and gluon densities. 

Many groups have presented results of NLO global fits of all available gi data 
available at that time. Relatively recent examples are reported in [34,53-57]. Figure 6 
shows the results of one recent fit [53] that accounts for contributions to gi that scale as 
1/Q 2 , which are significant for the data measured at Q 2 less than a few GeV 2 . As the 
photoabsorption cross section is independent of the sign of the quark charge, Am and 
Ad are unconstrained by inclusive data. Hence for this fit the convention was adopted 
of a flavour symmetric sea: 

Au sea = Au = Ad sea = Ad = As = As. (9) 

It is clearly determined that up quarks have a substantial positive polarization over 
the entire measured x range, while down quarks have a somewhat smaller and negative 
polarization. The polarization of the strange quark sea is found to be small but clearly 
negative. The gluon helicity distribution AG(x) is still poorly constrained, even though 
many new precise data have recently been added to the global set. In fact there are 
two alternative solutions (local minima in \ 2 ) °f equal quality with opposite signs of 
the gluon density, and yet another of similar quality using a bipolar shape for the x 
dependence. The reason for the large uncertainty is that the g\ data set is still far less 
extensive and precise than the unpolarized set, and covers a much more limited range 
in Q 2 (see section 2.8). The first moments of the helicity densities determined from 
the fit are shown in table 1. (These moments have some sensitivity to the assumed 
functional form of the densities, especially in the unmeasured small-x region.) The 
flavour- singlet moment AS is found to be in the range 0.2 to 0.25, implying that the 
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Table 1. First moments of the helicity densities resulting from the NLO QCD fit to all 
available g\ data of reference [53] in the (MS) factorization scheme, at the input scale of 
Q 2 = 1 GcV 2 . The uncertainties shown are total (statistical and systematic combined 
in quadrature). The parameters marked by (*) are fixed by the assumption of SU(3) 
flavour symmetry mentioned in the text. The Bjorken sum rule is also enforced. The 
table is taken from [53]. 







AG > 


AG < 


DoF 




826 - 20 


826 - 20 


x 2 




721.7 


722.9 


X 2 /DoF 




0.895 


0.897 


Au v = Au — 


Am 


0.926* 


0.926* 


Ad v = Ad - 


Ad 


-0.341* 


-0.341* 


AS 




0.207 ±0.040 


0.243 ± 0.065 


As + As 




-0.063 ±0.005 


-0.057 ±0.010 


AG 




0.129 ±0.166 


-0.200 ±0.414 



total |(1 - AH) of the remaining contributions to the spin of the nucleon from the 
gluons and/or quark orbital angular momenta is substantial. The gluon moment AG is 
still poorly determined, but it appears that large values are disfavoured. 

This collection of increasingly precise data from all of the experiments also offers 
the opportunity to more severely test the BSR. It is difficult to efficiently employ all of 
the data while computing the Sum directly from measured values of g\ and g™ (or gf) 
using (8), because kinematic interpolation and evolution to a common value of Q 2 is 
needed. This in turn requires a paramerization of the x B dependence at some 'starting' 
value of Q 2 , just as in the fits described above. In fact, those fits can be used to test 
the BSR, simply by omitting the usual BSR constraint from the fit, and computing the 
Sum from the resulting Aq(x) distributions produced by the fit (while accounting for 
uncertainty correlations). The relevant combination providing the structure function 
difference is the 'flavour non-singlet' distribution: 

Aq NS (x) = $>J/(e 2 > - 1) (Aq f (x) + Aq f (x)) . (10) 
/ 

The removal of the BSR constraint from the fit may result in large uncertainties 
for individual densities, but their uncertainty correlations should result in a smaller 
uncertainty for Aq^s- This approach might at first seem bizarre, because the sum rule 
is believed to be much more reliable than the assumption of SU(3) flavour symmetry used 
in the fits! However, it can be expected that Aq^s will be insensitive to the constraint 
from hyperon decay data if the proton and deuteron data have comparable precision 
and kinematic coverage, which is the case. (This can easily be verified by varying this 
input to the fit to study the sensitivity.) Hence for this purpose, the fitting process 
provides mainly the evolution of the data points to a common value of Q 2 , as well as 
the model used to extrapolate to x — 0, 1. This approach leads to good agreement with 
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the BSR within the total uncertainty from the fit, which amounted to about 10% of the 
sum, using the data available in 1998 [58]. Much more precise data are available today. 

There are several technical issues arising in the QCD fits that can be kept in mind 
when considering their results. First, typically the values taken from experiment are 
those for g\/ F\, even though the unpolarized structure function Fx has no relevance, 
and the systematic uncertainties for g\/Fi include a significant contribution from the 
poorly known ratio R = Cl/ct of cross sections for longitudinal and transverse photons, 
which has no effect on g\ itself. The reason given for this choice is the fact that 
experimental groups may have extracted g\ from the measured asymmetries effectively 
using values for the unpolarized cross section that were calculated using various values 
of R. However, as it is the ratio g\jF\ that is sensitive to R and not gi, it appears 
that this choice aggravates the problem. (Fortunately, the effects of this may be 
small.) Secondly, it can be difficult to obtain reliable results from multidimensional 
fits of data with limited statistical accuracy. Evidence for such difficulty here is the 
irregular shape of the dependence of \ 2 on some parameters that are reported. Up 
to three minima have recently been seen with similar values of \ 2 [53]. A common 
method of testing the stability and accuracy of algorithms for extracting results from 
experimental data is to analyze Monte Carlo pseudo-data. The situation here could 
be clarified if g\ pseudo-data were generated from a chosen set of polarized parton 
distributions, which could be compared to the results of the fit. Thirdly, two groups 
report statistical uncertainties that are inflated by approximately the square root of the 
number of fitted parameters (assuming that the \ 2 surface is quadratic in a substantial 
region near the minimum), so that their uncertainties cannot be directly compared 
to those of other groups. These inflated uncertainties do not correspond to what 
is understood by most of the community as a statistical uncertainty. (One way of 
defining this conventional uncertainty is the standard deviation of the distribution of the 
results that would be derived from similarly fitting a large number of Monte Carlo data 
sets resembling the global experimental data set, but with each data point fluctuating 
independently according to its experimental statistical uncertainty.) Fourthly, some 
groups account for the (substantial) experimental normalization uncertainty that is 
common to an entire data set from one experiment (arising from e.g. uncertainty in 
polarization of beam or target) by adding that common component of the systematic 
uncertainty in quadrature to the statistical uncertainty of each data point. This 
incorrectly deprecates the information available in those data. It is straight-forward 
to more appropriately account for the normalization uncertainty by means of a \ 2 
penalty term driven by a renormalization factor for that data set. Fifthly, few groups 
(see e.g. [55,58]) estimate 'theoretical' uncertainties for their results, arising from e.g. 
the somewhat arbitrary choice of functional forms to parameterize the helicity densities, 
and the dependence of the results on factorization and renormalization scales, indicating 
incomplete convergence of the perturbative expansion. 

We have seen that the global data set for the spin structure function g\ on the proton 
and neutron target has firmly established that quark helicities account for only a fraction 
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of the proton spin, leaving substantial contributions to be made by gluons and quark 
orbital angular momenta. The large uncertainty in the moment AG from global fits to 
gi data has motivated several major efforts to directly measure AG(x) using processes 
in which gluons enter in leading order, as described in section 2.8. Furthermore, 
information about quark orbital angular momentum is beginning to emerge from both 
calculations on the lattice and measurements of cross section differences or asymmetries 
for the deeply virtual Compton scattering process, as described in sections 4 and 5.2. 

2. 5. Helicity densities of sea quarks from semi-inclusive DIS data 

As mentioned above, quarks and antiquarks of the same flavour cannot be distinguished 
in inclusive DIS measurements, since the contributions of the various quark flavours to 
the cross sections scale as the square of the quark charges and hence are independent 
of the sign. Furthermore, the strange sea can be separated only in terms of first x- 
moments, because of the necessary constraint from hyperon beta decay data. In order 
to obtain information about the x dependence of the helicity densities of the various 
flavours, particularly of sea quarks, it is necessary to use additional types of experimental 
information. In semi-inclusive scattering (see figure 1), an energetic hadron is detected 
in coincidence with the scattered lepton. If suitable kinematic criteria are applied in 
the selection of the hadrons, they are likely to be members of the 'jet' of hadrons 
produced in the fragmentation of the struck quark, the non-perturbative process in 
which the energy of the quark is shared among hadrons created through the excitation 
of quark-antiquark pairs from the vacuum. For our purposes here, a key point is that 
the identities of the hadrons are statistically correlated with the flavour of the struck 
quark. For example, a ir + meson having valence composition ud is more likely to appear 
in the fragmentation of a u or d struck quark, and kaons are more likely to arise from 
strange quark fragmentation. However, these probabilities are also modulated by the 
square of the charges of the various flavours and their (known) spin-averaged quark 
densities. Hence the exploitation of this method is limited to some degree by the 'it- 
quark dominance' arising from both of these effects. For example, ir~ = du mesons 
and even K~ = us mesons, which provide good sensitivity to sea quarks, arise more 
often from fragmentation of u quarks. Hence a quark flavour decomposition based on 
data for various identified hadrons requires the solution of a strongly coupled system of 
relationships, which in turn requires primary data of high precision. 

Given an adequate understanding of the fragmentation process, a complete flavour 
decomposition of the quark and antiquark helicity distributions can be extracted from 
sufficiently precise measurements on both proton and 'neutron' targets of double-spin 
asymmetries in the cross sections for leptoproduction of various types of hadrons. 
The extraction requires knowledge of the probabilities for the various types h of 
hadrons to appear in the fragmentation of a struck quark of a given flavour /. These 
probabilities are embodied in the (familiar unpolarized) fragmentation functions Df(z), 
where z = E^jv and v and Eh are the energies in the target rest frame of the absorbed 



Spin-polarized high- energy scattering of charged leptons on nucleons 



22 



virtual photon (and hence of the struck quark) and of the detected hadron. While 
these fragmentation functions have typically been extracted [59, 60] from mostly high 
energy e + e~ collider data, a recent extraction [61] includes preliminary semi-inclusive 
DIS multiplicities of identified hadrons [62]. 

For leptoproduction of a hadron h, the semi-inclusive cross section difference that 
is analogous to the inclusive case of (2) can be factorized in leading order in a s as 



If, for example, there are available experimental measurements of these cross section 
differences for production of 7r + , 7r~, K + and K~~ mesons on both proton and 'neutron' 
targets, this system of eight over-constrained equations can in principle be solved for 
the helicity densities of up to six light quark flavours. As mentioned above, in practice 
w-quark dominance limits the sensitivity to some sea quark flavours. 

As in the case of inclusive measurements, absolute cross section differences are 
difficult to measure, especially with polarized targets. Hence cross section asymmetries 
are extracted from the data. Figure 7 shows all of the measured semi-inclusive double- 
spin asymmetries on both hydrogen and deuterium targets. It includes data from SMC 
for unidentified hadrons [43], from Hermes for production of identified 7r + , 7r~, K + and 
K~ mesons [63], and from Compass for unidentified hadrons [64]. Each collaboration 
extracted helicity densities in leading order in a s . In the case of Hermes, all three sea 
quark flavours could be distinguished, while fragmentation functions were effectively 
fitted to hadron multiplicities measured by the same experiment. It was found that 
the helicity densities of all three sea quark flavours are consistent with zero within the 
measured x range. Further analysis of all these data were later done in NLO, including 
data from other processes. This is described in the next section. 

2.6. Helicity densities from a global analysis in NLO 

The most complete extraction of helicity densities from all relevant published data 
appears in [65]. It is a global NLO QCD analysis like that described in section 2.4 for 
inclusive DIS data, but with the crucial inclusion of semi-inclusive DIS data to provide 
sensitivity to individual sea quark flavours, and for the first time RHIC + data 
for inclusive production of neutral pions and jets, providing much stronger constraints 
on the gluon helicity. Another new feature of the fit is the use of the above-mentioned 
set of fragmentation functions [61] based on not only e + e~ data but also preliminary 
semi-inclusive DIS multiplicities of identified hadrons. 

The extensive data set allows some freedom in the application of the constraint 
of neutron [hyperon] beta decay data through the assumption of SU(2) [SU(3)] flavour 
symmetry. Penalty terms were added to x 2 to limit violations of SU(2) and SU(3) 
symmetries, but their values at the solution are negligible compared to unity, suggesting 
that, at least in terms of first moments, the data are fully compatible with these 
symmetries. Because the data are still unable to distinguish s and s flavours, the 
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Figure 7. Semi-inclusive asymmetries for virtual-photon-nucleon DIS on both the 
proton (p) and deuteron (d) , producing the type of hadron labelled in the superscript 
on the symbols in the figure ('h' means unidentified hadrons). These data from 
SMC [43], Hermes [63] and Compass [64] are compared to curves computed from 
hclicity densities fitted to these and other data in NLO [65]. The acronym 'DNS' 
refers to the previous fit by this group [56, 66] using older information [67] about 
fragmentation functions. The figure is taken from [65]. 



assumption is made that As = As. (A previous analysis [63] in LO found that only As 
is constrained by the Hermes data, and the above assumption has negligible effect on 
that analysis.) 

Well-defined unique minima appear in the dependence of x 2 on each of the three 
sea quark flavours, and the quality of the fit is excellent: x 2 / DoF = 0.88. Figure 7 
shows a comparison of the resulting fits with the available semi-inclusive asymmetries. 
Asymmetries calculated from the same fragmentation functions but with helicity 
densities from a previous fit based on an older parameterization [67] of fragmentation 
functions are also shown, illustrating the sensitivity of the kaon asymmetries and hence 
sea densities on this input to the fit. 
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Figure 8. The continuous and dashed curves represent parton helicity densities 
multiplied by x, extracted in NLO from double-spin DIS and p^ + p^ data 
and evaluated at Q 2 = 10 GeV 2 , with their inner and outer uncertainty bands 
corresponding to S\ 2 = 1 and 5\ 2 = 2% of \ 2 respectively. The dotted and long- 
dashed curves represent the previously known unpolarized densities, which act as 
positivity limits to the helicity density limits through the requirement that all cross 
sections be non-negative. The dashed (dot-dashed) curves represent the results of the 
previous fits of [56] ( [54]). The figure is taken from [65]. 



Figure 8 shows the resulting parton helicity densities with uncertainty bands, while 
table 2 shows their (truncated) first moments. The densities Am + Am and Ad + Ad 
are not shown because they were already well constrained by the inclusive data. The 
semi-inclusive data constrain the sea densities, especially As, where a solution to a long- 
standing mystery appears. Leading-order analysis of semi-inclusive data has indicated 
that the strange helicity density is small or even positive at large x, whereas all analyses 
of only inclusive data imposing SU(3) symmetry clearly favour negative values of the 
first moment of As + As. These findings are now reconciled by a bipolar shape of As 
that is negative at small x and positive at large x. However, figure 8 shows that the first 
moment As is negative only because of the substantial contribution from the region 
x < 0.01, and table 2 shows that even x < 0.001 contributes strongly. It may seem 
mysterious that the fit forces As to be negative in a region of small x devoid of data, 
while the penalty term for SU(3) violation is so small at the solution point. However, 
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Table 2. Truncated first moments J x dx Aq(x, Q 2 ) at Q 2 = lOGcV 2 , taken from 
[65]. It is important to note the factor of two relating As in this table (assuming 
As = As) to As + As in table 1. The uncertainties for A^ 2 = 1 account only for 
experimental statistical and systematic uncertainties added in quadrature, and not for 
'theoretical' or model uncertainties. 
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+0.028 
-0.031 


AG 
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0.013 


+0.106 
-0.120 


0.013 


+0.702 
-0.314 


AS 


0.242 


0.366 


+0.015 
-0.018 


0.366 


+0.042 
-0.062 



the explanation may be that the As distribution is completely unconstrained in this x 
range by anything but the penalty term, so a very small value of the penalty is sufficient. 

The gluon helicity density is constrained, primarily by the addition of the RHIC 
data, to be significantly smaller than previous uncertainties, in the experimentally 
accessible x range. Table 2 shows that the first moment AG of the possibly bipolar 
distribution is consistent with zero, still with rather large uncertainties. 

A quantity of particular interest that is accessible to semi-inclusive DIS data is the 
flavour asymmetry Au(x) — Ad(x). Naively, one might expect that since u and d quarks 
have similarly small masses and they both are produced by gluon 'splitting' (g — > qq), 
both their densities and polarizations would be similar. However, the difference in 
their unpolarized densities u(x) — d{x) is experimentally well established [68], and is 
explained by various non-perturbative models that also predict an asymmetry in the 
helicity densities. Such a delicate imbalance in nucleon structure provides an excellent 
test of such models. The flavour asymmetry Au — Ad in the first moments can be 
compared to the predictions of two representative models. One is the chiral quark 
soliton model (xQSM) [69], which is an effective theory where baryons appear as soliton 
solutions of the chiral Lagrangian, and the other is a 'meson cloud' model [70] that 
considers quantum fluctuations of the spin-| nucleon into a virtual pion with spin 1 
plus a A hadron with spin |. The value Au — Ad = 0.15^0^ from the fit is consistent 
with the predictions of both models, but with a preference for the xQSM. 

2.7. Comparison of first moments of helicity densities 

Two previous sections discussed three significantly different but not independent sources 
of information about first moments of helicity distributions: the direct interpretation 
of the first moment of the spin structure function gf of the deuteron, NLO QCD fits 
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of parameterized helicity densities to all g\ data (see section 2.4), and global fits of 
helicity densities to all available relevant data, including double-spin asymmetries for 
both semi-inclusive DIS and inclusive production of neutral pions and jets in p^ + p^ 
(section 2.6). The first moments from these analyses are compared in table 3 on page 104. 
The comparison of results from the NLO QCD fits should be undertaken with caution, 
as the analyses differ in several respects. The LSS fit uses only inclusive DIS data, 
including the Compass data on the deuteron, which leads to two equally probable 
solutions. They are merged into a range of values given in the table rather than central 
values with standard deviations. The AAC fit does not yet include the Compass data, 
but includes some double-spin asymmetries for inclusive production of neutral pions in 
p=> -\. m Finally the DSSV fit includes all inclusive and semi- inclusive DIS data, as 
well as double-spin asymmetries for both semi-inclusive DIS and inclusive production of 
neutral pions and jets in p^ + . Another important difference between the analyses 
is the treatment of experimental statistical uncertainties, indicated in the first column. 
The LSS uncertainties are based on the standard treatment specified by the PDG [71], 
while the AAC group adopts a different treatment that increases their uncertainties. 
The DSSV group provides the results from two treatments. Finally, the common 
elements of the data sets used for the various results imply that the uncertainties are 
correlated, so that they should in principle be in better agreement than indicated by 
their uncertainties. 

The values for AS are all consistent, with the possible exception that the DSSV 
result for the 'standard' treatment of statistical uncertainties (A% 2 = 1) appears to 
be somewhat smaller. However it is important to bear in mind that the contribution 
—0.124 to that result from the unmeasured region x < 0.001 is much larger than this 
'discrepancy'. 

Having accounted for the factor of two relating the DSSV value for As (assuming 
As = As) to the values for As + As from the other two fits, all of the results again are 
consistent, again with the possible exception that the DSSV result for A% 2 = 1 appears 
to be somewhat more negative. However, essentially all of this value is contributed by 
the unmeasured region x < 0.001. The uncertainty of the DSSV result is much larger 
than that of AAC, in spite of the additional data used in the DSSV fit, presumably 
because of the relaxation of SU(3) flavour symmetry. 

Finally, all the results for AG are consistent with each other and with zero. The 
additional RHIC data included in the DSSV fit resolves the ambiguity arising in the 
LSS fit, and much improves the precision. Here the contribution to that result from the 
unmeasured region is about —0.07. Other (leading order) experimental constraints on 
AG are discussed in the next subsection. 

2.8. Gluon helicity distribution 

The nucleon spin can be decomposed into contributions of quarks and gluons, as will 
be explained in section 5.1. A decomposition in which all terms have a probabilistic 
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a b c 

Figure 9. Feynman diagrams for hard subprocesses in lepton-nucleon scattering: 
a) 0(a°) DIS (7**7 -> q), b) 0{a\) Photon-Gluon Fusion (7*5 -> gq), c) O(aJ) QCD 
Compton scattering (7*<3 I — ► 93). 

interpretation as parton densities in the helicity basis can be written as: 

i = iAS + AG + £, (12) 

where C is the total orbital angular momentum of partons. The intrinsic gluon 
contribution AG is accessible to experiment. Until recently the only existing knowledge 
about the gluon helicity distribution AG(x) was derived as its first moment AG = 
J dx AG(x), from NLO pQCD fits to the spin-dependent structure function gx(x, Q 2 ) 
measured in double-spin lepton-nucleon scattering experiments (see section 2.3). 
However, this method is subject to large uncertainties, as gi(x,Q 2 ) is only weakly 
sensitive to AG(x) via the Q 2 evolution of parton distribution functions. Very recently, 
better precision was achieved in a global analysis of DIS and pp data, as explained in 
section 2.6. 

More experimental data are now emerging from the study of the subprocess in 
which gluons appear in leading order in double-spin lepton-nucleon scattering: Photon- 
Gluon Fusion (PGF), see figure 9(b). Here, the gluon emits a qq pair in which each 
quark inherits the helicity of the gluon. As explained in section 1, the (anti-) quark 
can absorb the photon only if they have the same helicity. Given positive (negative) 
values of AG(x), the PGF process has a negative (positive) asymmetry, while the major 
background processes DIS (figure 9(a)) and QCD Compton (figure 9(c)) have always a 
positive asymmetry. 

The PGF process produces a qq pair resulting in 2 jets of hadrons. At fixed-target 
energies, one has to use a leading hadron to represent a jet, plus possibly another 
hadron separated by large transverse momentum to represent the other jet. In order 
to suppress background processes, final states are selected which contain either large 
transverse momenta or charm. These signatures are also exploited to access AG(x) in 
double-spin proton-proton scattering in the experiments Star and Phenix at Rhic. A 
cleaner channel in these experiments will be single direct photon production. 

The gluon polarization ^=r{x) is extracted in a leading-order analysis from the 
measured double-spin asymmetry A\\ in the hadron production cross section. Results 
on 4r(x) are available from SMC, Compass, and Hermes. In particular, Compass 
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Figure 10. Asymmetry measured at Hermes for positive (left) and negative (right) 
inclusive charged hadrons on the proton (top) and deuteron (bottom), as a function of 
PT(beam) ■ The errors are statistical only. There is an overall experimental normalization 
of 5.2% (3.9%) for proton (deuteron) data. The curves show the asymmetries 
as expected from the Pythia [73, 74] Monte Carlo simulation for three different 
assumptions on the gluon polarization. This figure is taken from [75]. 



uses charmed meson decay products and pairs of high-fc^ hadrons, while Hermes uses 
single charged high-Zcy hadrons. The hard scale is given by the charm quark mass in 
the former, and by kx in the latter case. Using as an example Hermes data for h ± - 
production on the proton and deuteron, respectively [75], A\\ is shown in figure 10 as 
a function of hadron transverse momentum PT(beam) that is measured with respect to 
the beam axis. Also shown is the dependence expected from a Monte Carlo simulation 
for several values of (0, ±1), which demonstrates that for PT(beam) > 1 GeV the 
asymmetry clearly starts to depend on the underlying value for ^p. Comparing it to 
the data reveals at the level of the asymmetry, that small values of ^f- are prefered. 
In order to enhance the sensitivity, the 'detected' fragmenting quark should have the 
largest possible transverse momentum. This is then reflected in large PT(beam) values for 
the outgoing hadron(s). 

In figure 11 the above mentioned results on -gr(x) are shown. The three statistically 
most precise ones are from HERMES on single hadron production on the deuteron 
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Figure 11. Leading-order results on ^-{x), published or preliminary, from 
Compass [76-78], Hermes [75,79], and SMC [80]. The sequence of references 
corresponds to the order of measurements shown in the legend, where also the 
analysed final state is indicated. Outer error bars, if shown, represent statistical and 
systematic uncertainties combined in quadrature. Inner error bars represent statistical 
uncertainties only. Additionally, two equivalently good QCD NLO fits of the structure 
function gf are shown [72], obtained with the constraints AG < and AG > 0, 
respectively The figure is taken from [76], updated with [78]. 

at very low Q 2 [75] and from COMPASS on high-transverse-momentum hadron-pair 
production on the deuteron at Q 2 < 1 GeV 2 [77] as well as at Q 2 > 1 GeV 2 [78]. 
For the Hermes result, very detailed studies were performed to estimate systematic 
uncertainties associated with Monte Carlo models of the background processes. The 
major contribution here stems from the uncertainties in the fractional contribution from 
soft processes, which are poorly known. Among the Compass results on high-p t hadron 
pairs, the most recent also includes a study of systematic uncertainties associated with 
the Monte Carlo model. 

From the three high-p t hadron results, it can be concluded that ^p(x) at medium- 
x (0.06. ..0.3) is about 0.1 with a total uncertainty of the same magnitude. This 
leading-order result is in fair agreement with results of NLO fits of g\ from [72] shown 
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as curves in figure 11. The value of the unpolarized gluon distribution from the 
CTEQ6 set of PDFs [81] is G(x)\ x= o,i ~ 10, so that for the gluon helicity distribution 
AG(x)|a; = o.i ~ 1 ± 1 is obtained. Expressed as x ■ AG(x), the LO result from high-p t 
hadrons is 0.1 ± 0.1, in agreement with the less model-dependent result from the NLO 
QCD fits of [65] discussed in the previous section (see bottom right panel of figure 8). 
We note that any comparison of next-to-leading-order to leading-order results has to be 
done with caution since k factors, defined as their ratio, are known to be large. In the 
kinematic region considered here, k factors were found [82] of up to 5 for cross sections 
and up to 2 for asymmetries that are used to extract information on the polarized gluon 
distribution. 

If and when it becomes possible to experimentally constrain the full first moment 
AG with sufficient accuracy, it will be possible to determine the total orbital angular 
momentum C of partons according to (12). 

3. The transverse spin structure of the nucleon 

3.1. Transverse spin in deep-inelastic scattering 

The quark transversity distribution Sqf(x) was defined in section 1 as the difference 
between number densities of quarks with transverse polarization in the same and 
opposite direction as that of the nucleon spin, in a nucleon that is polarized transverse 
to its 'infinite' momentum. It may at first appear mysterious that it is necessary to 
distinguish between Sqf(x) and the helicity distribution Aqf(x), since there would be 
no need to distinguish between them in a nonrelativistic system. 

All twist-2 parton distributions are defined as light-cone Fourier transforms of 
products of quark fields ipf (or gluon fields). The only difference between qj(x), Aqf(x), 
and Sqf(x) are the Dirac matrices that appear between xjjl and ifjf. For Aqj the 
Dirac matrices are (in Dirac representation, where the lower components vanish in the 
nonrelativistic limit) 

*v-0 °0U «)■ <i3) 

while for 5qf they read 

(:;)-(»-)■ (i4 » 

So the only difference between Aqf and Sqj is the minus sign in front of the lower 
components. In the nonrelativistic limit, the lower components vanish like p/m, so the 
difference between Aqf and 5qj vanishes like v 2 /c 2 . 

Transversity loses its probabilistic interpretation as a difference in number density 
distributions when expanded in a helicity basis. Every observable in quantum mechanics 
has a probabilistic interpretation as the sum of eigenvalues of the observable weighted by 
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the probability to be in each eigenstate.f However, if basis functions are used in which 
the observable is not diagonal, the expectation value of the observable becomes a sum 
over products of amplitudes and the probabilistic interpretation is 'lost', or at least not 
apparent. In the present context this elementary principle appears most clearly using 
Pauli matrices and identifying eigenstates of a z or a x with helicity eigenstates | <->•) 
or transversity eigenstates | |) = -j= (| —>) + e %ip \ <—)), respectively, where (p specifies 
the orientation of the transverse polarization direction. Obviously o~ x is diagonal when 
expressed in terms of its own eigenvectors, but loses this property when expressed in 
terms of the eigenvectors of a z : 

*» = ITXTI - 14X41 = l<->H + HXH- 

Therefore the expectation value of the 'transversity operator' o~ x in some state ip becomes 
a sum of products of amplitudes in the helicity basis: 

= M T)| 2 - \(rf>\ |)| 2 = M + M — > <^l «->*■ (16) 

The right hand side of (16) is not a sum over absolute squares and hence the probabilistic 
interpretation for (ip\a x \ip) is no longer apparent. 

Despite the loss of a direct probabilistic interpretation, there are (at least) two 
reasons for discussing transversity in a helicity basis. First, that basis has already been 
used to study the more familiar longitudinal spin polarization, and secondly only helicity 
and not transversity is conserved in hard processes. The second statement is based on 
the observation that for negligible quark masses, both the QCD and QED Lagrangians 
are invariant under a 'chiral' transformation of the quark spinors: 

ifj f (x) — ► e is ^^ f (x). (17) 

As a consequence of this 'chiral symmetry', the chirality of a massless quark (the 
eigenvalue of 7 5 ) is conserved. It turns out that for a massless quark its chirality 
agrees with its helicity, but for an antiquark with the opposite sign. J Chiral symmetry 
thus implies the conservation of the helicity of a massless quark in hard scattering 
processes. However, even though helicity is conserved, transversity is not, as the 
scattering amplitudes for the two helicity states may differ in their phase. An analogue 
from optics is the propagation of light through an optically active medium where the 
refraction indices for the two circular polarization states are different. While the circular 
polarization does not change as the beam propagates through this medium, the linear 
polarization, which is related to the longitudinal polarization through a simple linear 
combination, precesses around the beam axis since the waves describing the two circular 
polarization states acquire different phases in such a medium. 

An important distinction must be made between transversity and transverse spin. 
Transverse spin is the expectation value of the operator ^7 75^ and is chiral-even, 

f This is one of the postulates of quantum mechanics. 

\ This minus sign for antiquarks also illustrates that, even for massless quarks, chirality and helicity 
are not the same in the sense of being described by the same operator, which is one reason why it is 
necessary to introduce the concept of chirality in addition to helicity. Another is that the helicity of a 
massive particle is frame dependent, whereas chirality is always Lorentz- invariant. 
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whereas chiral-odd transversity is the expectation value of the operator ip'a^^y^ip, which 
has matrix elements between states of opposite quark helicity in the proton rest frame. 
Transverse spin does not have a partonic interpretation in the light-cone framework 
and is associated with the (twist-3) spin structure function g 2 , whose effects vanish like 
1/ \/Q^ in the Bjorken limit, whereas transversity does not vanish in this limit. For a 
detailed discussion of these subtle issues see [83]. 

Since hard electromagnetic (or weak) processes conserve the helicity of the quark, 
inclusive DIS experiments cannot be used to measure the transversity of the quarks. 
In order to see why, we first use the optical theorem in a transversity basis and then 
express transversity eigenstates in terms of helicity eigenstates. The difference between 
the contribution to the cross section from quarks with the same and with opposite 
transverse polarization as the nucleon would be proportional to the imaginary part of a 
forward scattering amplitude involving quark helicity flip (see figure 12). However, since 





oc Im 








Figure 12. When expressed in terms of quark helicity states, the dependence of the 
DIS cross section on the transversity of the quarks is proportional to the imaginary 
part of the quark helicity- flip forward Compton amplitude which vanishes at large Q 2 . 



absorption of a hard vector boson (e.g. photon) cannot flip the helicity of a massless 
quark, this forward scattering amplitude must be zero and therefore inclusive DIS cannot 
be sensitive to the transverse spin of the quarks. Since this argument does not involve 
the polarization of the photon in DIS, this result also illustrates why a transversely 
polarized charged lepton beam does not help to determine the transverse polarization 
of the quarks in inclusive DIS. In semi-inclusive DIS, the meson production vertex can 
flip quark helicity, which is why the latter process can be sensitive to transverse quark 
polarization. 

It is instructive to compare DIS mediated by vector boson exchange, where quark 
helicity is conserved, with a Gedanken DIS experiment involving scalar bosons, where 
quark helicity flips at the hard quark-scalar vertex. If only scalar bosons are exchanged, 
then the quark helicity flips at both vertices in the forward Compton amplitude 
analogous to figure 12 and such a process would not be sensitive to quark transversity 
either. Only the interference term between vector and scalar boson exchange provides 
the single helicity flip that is needed to make this amplitude non- vanishing, i.e. the 
interference between photon exchange and the exchange of a hypothetical scalar boson 
(such as the Higgs boson) would in principle be suitable for accessing transversity 
experimentally using inclusive DIS [84]. 
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3.2. Tensor Charge of the Nucleon 

The tensor charge is obtained from the forward matrix element of the operator 
ipj.a^ipf, analogously to the vector and axial vector charge (see section 2.2). In terms 
of quark transversity eigenstates, the tensor charge involves the linear combination 



the antiquarks. The fact that the tensor charge contains the difference between quark 
and antiquark transversity is convenient as the contributions from short-lived virtual 
qq pairs cancel, and it can be identified with the transversity of valence quarks only. 
Moreover, the tensor charge can be computed in lattice gauge theory, and with better 
precision than the axial charge (see section 4.4). 

In a helicity basis, transversity involves both a quark and a nucleon helicity flip, 
i.e. it involves the overlap between a nucleon 'state' that had positive helicity and from 
which a quark with positive helicity has been removed and another nucleon 'state' that 
had negative helicity and from which a quark with negative helicity has been removed. 
The Cauchy-Schwarz inequality tells us that this 'scalar product' is less than the norms 
of the states involved. Since the latter represent the probability to find a quark with 
the same helicity as the nucleon spin, one thus finds the Soffer inequality [85] 



which holds independently for each quark and antiquark flavor. It provides useful 
constraints on parameterizations of transversity. 

Attempts have been made to develop a 'transversity sum rule' [86, 87] that is 
analogous to the 'Ji decomposition' of the nucleon spin (39) but involves the net 
transversity of quarks (plus antiquarks). However, one problem with any transversity 
sum rule is that quarks and anti-quarks contribute with opposite signs to the tensor 
charge (first moment of transversity), while they contribute to either longitudinal or 
transverse nucleon spin with the same sign. Therefore, transversity cannot contribute 
to such a sum rule in terms of the tensor charge, but rather in terms of the tensor 
charge for the quarks plus the anti-quarks — a combination that cannot be represented 
in terms of a local operator. Therefore lattice calculations will be unable to calculate the 
transversity contribution to a transversity sum rule. Furthermore, in the transversity 
sum rule in [86,87], neither an operator definition nor an experimental procedure for 
determining the orbital angular momentum contribution has been provided. Thus this 
sum rule cannot be tested and is unfalsifiable. 

3.3. Transversity in the Drell-Yan reaction 

The conceptually simplest practical way to measure transversity is to exploit the Drell- 
Yan reaction in doubly transversely polarized hadron-hadron collisions h[ +h 2 — > 
H + /J.~ + X, involving the hard scattering process 




the transversity of the quarks minus the transversity of 




(18) 



q + q 




A 4 A 4 



(19) 
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The high-energy cross section for this reaction is obtained by convoluting the quark 
distribution in one of the hadrons with the antiquark distribution in the other hadron 
and with the cross section for q + q — > , where the latter can be calculated in 
perturbation theory (one photon in the intermediate state plus higher-order corrections). 
Transversity can be studied in Drell-Yan experiments with transversely polarized (anti-) 
protons because the quark and antiquark preferentially annihilate into a vector boson 
when their transverse spins are parallel. Hence the double-spin asymmetry Att is 
capable of measuring the correlation between the spin of the quark in one of the hadrons 
with that of the antiquark in the other hadron. This double-spin asymmetry is therefore 
sensitive to the product of the transversity distributions for quarks and antiquarks 
carrying momentum fractions X\ and Xi respectively: 

a ff _ a U 

A T t(x 1 ,x 2 ) = oc 22 e } 5( lf/hi(xi)Sq f/h2 (x2). (20) 

/ 

The momentum fractions of the quark and antiquark can be reconstructed from the 
invariant mass and net longitudinal momentum of the muon pair into which they 
annihilate. The struck q in a Drell-Yan reaction can be either a sea or valence quark 
in an incident antiproton (or meson or proton). In the interpretation of the data for 
p+p — > fi + fi~ + X, one x value can be associated with the q by considering asymmetric 
events where the momentum fraction of one of the initial quarks was so high that it is 
very unlikely to have been an antiquark since the latter carry only small momentum 
fractions. As the Drell-Yan cross section depends on products of PDFs, determining the 
overall normalization may in general be ambiguous. In the case of p + p — > fi + fi~ + X 
there is no such ambiguity as the same valence PDF of the proton appears as a product 
with itself. Once this normalization has been determined, it can be used to determine the 
normalization of sea quarks in pp scattering. Plans are being laid at several laboratories 
to exploit this reaction to measure transversity [88,89]. 

All existing experimental constraints on transversity are based on the sensitivity of 
fragmentation to the transverse polarization of the struck quark in processes other than 
Drell-Yan, as manifested in single-spin asymmetries. These will be discussed in Section 
3.5. 

3.4- Intrinsic parton transverse momentum 

Partons can also have momentum px transverse to the infinite momentum of the hadron. 
Hence parton distributions can be introduced that depend not only on x but also on 
Pt [90,91] and have interpretations as number densities or differences between number 
densities [92]. Only three of these distributions survive integration over p? leading to 
the familiar PDFs that depend only on x. These have already been discussed in section 
1. They are the spin-averaged distribution g/(x, pr), the helicity distribution Aqy (x, p T ) 
and the transversity distribution Sqf(x,pT)- However, many other twist-2 distributions 
exist and they have probabilistic interpretations similar to that of integrated PDFs [92]. 
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These distributions§ are sometimes known as Transverse-Momentum-Dependent (TMD) 
parton distributions (TMD-PDFs, or simply TMDs). In the rest of this subsection we 
take the unpolarized case as an example. 

When a quark in a nucleon absorbs a hard virtual photon, the resulting momentum 
of the quark is the sum of its initial momentum and the momentum transfered from the 
lepton. In leading order a s , any momentum component transverse to the virtual photon 
direction thus originates from the intrinsic transverse momentum p? of the parton before 
the scattering, described by qf(x,pr)- 

The distribution of hadrons h in a jet originating from a fast quark of flavor / 
is described by fragmentation functions Df/h(z,kr)- Here &t = Ph,t — zpr is the 
transverse momentum of the hadron relative to the momentum of the fragmenting 
quark. In DIS, the factorized semi-inclusive cross section for producing a hadron with 
transverse momentum p^T with respect to the virtual photon momentum is obtained by 
convoluting|| the transverse-momentum-dependent parton distributions Qf(x,pr) with 
these fragmentation functions: 

cr ep ^ ehX (x B ,z,p h>T ) cc^2q f (x B ,p T )®D f/h (z, k T ) ® 6 (zp T +k T -p h>T ). (21) 

f 

These fragmentation functions are expected to be universal, i. e. independent of both 
process and target [93]. For example, the distribution of hadrons in jets produced in 
e + e~ — > hadrons is described by the same fragmentation functions that describe jets 
in DIS. However, as discussed in section 3.6, initial-state interactions (ISIs) and final- 
state interactions (FSIs) need to be considered before the concept of universality can be 
literally applied. 

3.5. Semi-inclusive DIS and the Collins effect 

When a transversely polarized quark fragments into hadrons, the structure of the jet is 
sensitive to the polarization of the quark. In particular, the orientation of the transverse 
spin direction together with the quark's momentum defines a plane relative to which 
the momentum distribution of hadrons may not be left-right symmetric. To understand 
how this may occur, we consider a simple model for the case of fragmentation into 
pions: in this process for example a quark with spin pointing up needs to pick up from 
the vacuum an anti-quark with spin down in order to form a spin zero state. The 3 Pq 
model for quark anti-quark pair creation [94] postulates that the qq pair is created with 
'vacuum quantum numbers' J p = + , the simplest configuration of which is S = 1, 
L = 1 and no gluon. This model is based on the idea that as the QCD string between 
the fragmenting quark and the target remnant is stretched and stores energy, a qq pair 
can tunnel from the vacuum to remove some of the stored energy and the string breaks. 
Due to angular momentum conservation, this model suggests that the qq pair is thus 

§ In the 'Amsterdam notation', distributions that do not 'survive' integration over d 2 px are signified 
by the superscript _L. 

|| Only for Gaussian dependences can these convolutions be factorized. 
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Figure 13. Behind the leading quark (qo) moving along z a string forms, which 
eventually breaks through the formation of the qiqi pair. In the 3 Po model, the qq 
pair has orbital angular momentum in the opposite direction as the spin of the qq pair, 
i.e. in the case of fragmentation into pions (spin 0) in the same direction as the spin 
of the leading quark. Since the antiquark from the pair is expected to be produced 
closer to the leading quark as the produced quark, this implies that the produced 
antiquark inherits a positive transverse momentum kq ± x in the y direction from the 
orbital motion (for a leading quark polarized in the y-direction). 



produced with orbital angular momentum in the same direction (i.e. up in this case) as 
the spin of the fragmenting quark in the produced pion (see figure 13). The antiquark 
that merges with the fragmenting quark should retain some of that transverse orbital 
angular momentum, causing it to move to the left (when looking into the direction of 
motion of the fragmenting quark with spin up) [95]. The pion inherits the transverse 
momentum carried by the antiquark. 

The fragmentation function H^, h (z, hr) describing the left-right asymmetry in the 
fragmentation of a transversely polarized quark of flavor / into a hadron h was identified 
by Collins [96], and found to be chirally odd. Its combination with transversity Sqf(x) 
that describes the distribution of transversely polarized quarks thus gives rise to a 
correlation between the transverse target spin and the transverse momentum of the 
produced hadron(s). The associated asymmetry is proportional to sin (0 — (f) q ), where 
(f) q is the azimuthal angle of the polarization vector of the fragmenting quark about the 
virtual photon direction and (p is the azimuthal angle of the trajectory of the produced 
hadron h. However, the polarization direction of the struck quark and that of the 
fragmenting quark are not identical. In QED, when the electron emits a high-momentum 
virtual photon, the electron helicity is conserved and therefore the orientation of the 
electron spin is rotated in the lepton scattering plane. The virtual photon inherits the 
change in angular momentum and is thus linearly polarized in the lepton scattering 
plane. When a quark absorbs this linearly polarized photon, its spin component along 
the polarization direction of the photon tends to flip. For a transversely polarized 
quark, this implies that its polarization direction is on average tilted symmetrically with 
respect to the normal of the lepton scattering plane. For azimuthal angles measured 
relative to the lepton scattering plane, as illustrated in figure 14, the tilt results in 
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Figure 14. The definitions of the azimuthal angles of the hadron production plane 
and the axis of the relevant component S±_ of the target spin, relative to the scattering 
plane that contains the momenta k and k' of the incident and scattered leptons. 

cf) q = it — 0s, where s denotes the polarization direction of the target nucleon, since the 
polarization of the quark before the photoabsorption is proportional to the product of the 
target nucleon polarization and the transversity distribution. The resulting sin (0 + 0s) 
angular dependence can thus be used to extract the product of transversity and the 
Collins fragmentation function. 

However, this 'Collins effect' is not the only mechanism that can give rise to 
a correlation between the target spin and the momenta of final state hadrons. A 
correlation between transverse momenta of quarks in the target and the target nucleon 
spin may give rise to a sin (0 — 0s) asymmetry, as we will now discuss. 

3.6. Sivers effect 

In a target that is polarized transverse to the virtual-photon direction {e.g. up), the 
quarks in the target nucleon can exhibit a (left/right) asymmetry of the distribution 
fq/pi ( x i Pt) in their transverse momentum 

U/A^Pt) = ff&Pr) - fiT q (^PT) {P X m ] ' S > ( 22 ) 

where S is the spin of the target nucleon and P is a unit vector opposite to the 
direction of the virtual-photon momentum. The fact that such a term may be present 
in (22) is known as the Sivers effect and the function f^(x,pj<) is known as the Sivers 
function [97]. The latter vanishes in a naive parton picture since (P x p T ) ■ S is odd 
under naive time reversal (a property known as naive- T-odd), where one merely reverses 
the direction of all momenta and spins without interchanging the initial and final states. 
In fact, until a few years ago, it was believed that full time reversal invariance of QCD 
forbids the existence of the Sivers effect. However, then it was realized that a FSI 
mediated by soft gluon(s) between the struck quark and the remainder of the target 
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nucleon, before the quark fragments into hadrons, could avoid this restriction [98]. This 
led to a generalization of our understanding of what is meant by a parton distribution 
and its universality, and how it may appear in expressions for matrix elements of 
processes. The momentum fraction x, which is equal to x B in DIS experiments, 
represents the longitudinal momentum of the quark before it absorbs the virtual photon, 
as it is determined solely from the kinematic properties of the virtual photon and the 
target nucleon. In this new view, in DIS the transverse momentum represents the 
asymptotic transverse momentum of the active quark after it has absorbed the virtual 
photon and then left the target but still before it fragments into hadrons. Thus the 
Sivers function for semi-inclusive DIS includes the FSI between struck quark and target 
remnant, and time reversal invariance no longer requires that it vanishes. Indeed, time 
reversal not only reverses the signs of all spins and momenta, but also transforms FSIs 
into ISIs. It has been shown that the Sivers function relevant for SIDIS and that relevant 
for Drell-Yan (DY) processes must have opposite signs [98], 

fir { x iPt) sidis = -fvr( x ^Pr)DY, (23) 

where the asymmetry in DY arises from the ISI between the incoming antiquark and 
the target. The experimental verification of this relation would provide a test of the 
current understanding of the Sivers effect within QCD. While the argument leading to 
(23) is nonperturbative, it is instructive to elucidate its physical origin in the context 
of a perturbative picture: for instance, when the virtual photon in a DIS process hits a 
red quark, the spectators must be collectively anti-red in order to form a color-neutral 
bound state, and thus attract the struck quark. In DY, when an anti-red antiquark 
annihilates with a target quark, the target quark must be red in order to merge into a 
photon, which carries no color. Since the proton was colorless before the scattering, the 
spectators must be anti-red and thus repel the approaching antiquark [99] . 

Another reason why the Sivers effect is attracting attention is that it relates to 
the correlation between the orbital angular momentum of partons and the spin of the 
parent nucleon. A semi-classical argument can help to understand how. Should quarks 
have orbital angular momentum, it would be correlated with the spin of the nucleon, as 
the latter defines the only preferred direction. When viewed in the nucleon rest frame, 
quarks orbiting in a transversely polarized nucleon move towards the virtual photon 
on one side of the nucleon but away from it on the opposite side. Upon boosting the 
nucleon to the infinite-momentum frame, the quarks on one side of the nucleon with 
a particular momentum fraction x thus appear to have larger longitudinal momenta 
than on the other side. Since parton distributions are rapidly falling functions of x, 
shifting them towards larger (smaller) x results in an increased (decreased) probability 
of interaction when viewed at a fixed value of x. The result is that, at a given value 
of x, more quarks absorb photons on the left or right side of the nucleon when the 
nucleon spin is pointing e.g. up. As the struck quark moves away, it experiences an 
attractive FSI that tends to 'focus' it towards the axis through the center of the nucleon 
along the direction of the virtual photon (see figure 15). Thus the excess of interactions 
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Figure 15. A correlation between the nucleoli spin and the quark orbital angular 
momentum, in combination with an attractive final state interaction, gives rise to a 
Sivers effect for it (d) quarks imparting a transverse momentum that is on average 
directed upward (downward). 

on one side results in a tendency of the struck quark and hence the hadron jet to be 
deflected towards the opposite side, producing a sinusoidal azimuthal distribution of the 
hadrons relative to the transverse polarization direction of the target nucleon. Quarks 
of different flavours may have different spin-orbit correlations, and also tend to produce 
different populations of hadron types. Hence different hadron types may have different 
azimuthal distributions. 

The type of deflection shown in figure 15 also arises in optics when a convex lens 
is illuminated asymmetrically, e.g. when only one side of the lens is illuminated with a 
laser beam parallel to the axis of the lens. The beam is then deflected towards the lens 
axis and the outgoing beam has a left-right 'momentum asymmetry' [100]. Because of 
this analogy, the above mechanism has been called 'chromodynamic lensing' [101]. 

As the asymmetry is caused by the orbital motion of the quarks, it is a (left-right) 
asymmetry relative to the direction of the (vertical) transverse target polarization, i.e. 
proportional to sin(0— 0$) in SIDIS. This differs from the sm(<t> +(/)$) angular dependence 
caused by the Collins effect, thus permitting the two contributions to be separated 
through the use of transverse target polarization. 

The Sivers function not only provides a clue about the orbital angular momentum 
structure of the nucleon wave function, but in addition carries important information 
about the complex phase structure of the FSI [102]. 



3.7. Boer-Mulders effect 



Even when the target nucleon is unpolarized, there can exist a correlation between 
the distribution in transverse momentum and transverse polarization of its quarks. 
The Boer-Mulders effect is a possible dependence of the distribution f q i / p (x,Pt) of 
transversely polarized quarks in an unpolarized target on the transverse polarization S q 
of the quarks: 



f q y P (x,p T ) 



fi(x,pr) - hf q (x,pr) 



[P x pr) ■ gg 

M 



(24) 
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Figure 16. Illustration for the interplay between transverse momenta from the 
Boer- Mulders and the Collins effect, leading to a cos2^i asymmetry: a) 'primordial' 
transverse spin distribution in an unpolarized target, b) transverse momentum from 
the Collins effect for hypothetical SIDIS without transverse spin flip, c) transverse 
momentum from the attractive FSI, d) after the virtual photon absorption, the 
transverse spin component in the (horizontal) lepton scattering plane has flipped, e) 
transverse momentum from the Collins effect including the transverse spin flip, f) net 
transverse momentum from FSI and Collins effect. 

and the Boer-Mulders function hf q (x,Pr) embodies this correlation [103], which can 
naturally arise from non-perturbative spin-orbit interactions. 

In the case of the Sivers effect, the relevant spin is that of the target nucleon, 
whereas here it is that of the struck quark. However, both effects can be understood 
to arise in part from the preference of the virtual photon to absorb on a quark with 
excess momentum fraction associated with its orbital motion, because of the monotonic 
decrease with x of the unpolarized PDF. Thus a quark found on the top (bottom) [left] 
{right} side of the nucleon, as seen by the virtual photon, will tend to have orbital 
angular momentum pointing right (left) [up] {down} (on first reading it may be helpful 
to consider only the first option and compare to figure 16). If the quark spin is correlated 
with its orbital angular momentum as represented by the Boer-Mulders function, then its 
polarization will have a similar (or opposite) spatial pattern (figure 16a). As described in 
section 3.5, the quark polarization affects the transverse momentum of produced hadrons 
via the Collins effect. Quarks with polarization pointing right (left) [up] {down} may 
be deflected e.g. up (down) [left] {right} , or the opposite depending on the sign of the 
Collins fragmentation function. We first consider the simpler case of exchange of a scalar 
boson rather than the vector photon, where the former cannot carry any information 
about the orientation of the lepton scattering plane so that the initial quark transverse 
polarization direction persists through the photoabsorption process. Then the combined 
effect would be a net 'focussing' or 'defocussing' of the produced hadrons relative to the 
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direction of the exchanged boson, according to the product of the Collins function with 
the Boer-Mulders function representing the original spin-orbit correlation of the quark 
(figure 16b). However, we must beware that this term '(de)focussing' can be misleading 
in this context. Since one side of the nucleon is not experimentally distinguished from 
the other, and the target nucleon is uniformly illuminated by the parallel 'beam' of 
virtual photons, the effect only randomly scatters or smears the direction of the hadron. 
At each azimuthal angle of the hadron direction about the virtual photon direction there 
is a net smearing effect, while the mean change of direction due to this '(de)focussing' 
is zero. 

An additional independent axially symmetric 'focussing' action is provided by the 
attractive FSI that is crucial also for the Sivers effect (figure 16c). (The Sivers effect 
itself plays no role here because the target nucleon is unpolarized.) This FSI effect 
either enhances or competes with that in figure 16b to produce a net '(de)focussing' 
action on the direction of the produced meson relative to that of the virtual photon. 
Thus for scalar exchange, the net effect may be a change in the cross section for hadron 
production that is independent of azimuth. For example, for observed values of the 
meson momentum component Ph,t transverse to the direction of the virtual photon that 
are larger than the most probable value, this smearing would increase the cross section. 

On the other hand, as also described in section 3.5, for vector (photon) exchange, the 
initial transverse quark polarization does not always persist through the photoabsorption 
process. Because of the linear polarization of the virtual photon in the lepton scattering 
plane, the component of the quark polarization in that plane can be flipped (figure 16d), 
inverting the effect in figure 16b only for hadron transverse momentum components 
normal to that plane (figure 16e). Now this effect tends to enhance the FSI effect in 
one plane and cancel it in the orthogonal plane (figure 16f). Since the net smearing 
effect is proportional to the magnitude of the net '(de) focussing' action, the smearing 
is different in the lepton scattering plane from that in the orthogonal plane. Hence the 
effect on the cross section can be different for hadrons produced in the lepton scattering 
plane compared to the orthogonal plane. This amounts to a cos(2</>) dependence of the 
hadron azimuthal distibution, relative to the lepton scattering plane. In the absence of 
the FSI, the inversion of the Collins effect in the lepton scattering plane (illustrated in 
figure 16e) would not change the '(de)focussing' magnitude, which is all that matters 
for the cross section. The necessary role of the FSI is associated with the T-odd nature 
of this Boer- Mulders effect. 

A more formal argument is also included here for the sake of completeness, as it 
may not be easily available in the literature. The Collins effect acts as a polarimeter 
with analyzing power oc H^~(x, /ct)(-P x &t) • S' q , which tags the transverse polarization 
S' of the quark after it has absorbed the virtual photon and its component in the lepton 
scattering plane has flipped compared to S q . Here &t is the transverse momentum of 
the outgoing hadron relative to the momentum direction of the fragmenting quark. The 
resulting asymmetry in the hadron distribution relative to the lepton scattering plane is 
obtained by convoluting the Boer-Mulders function with the Collins function to obtain 
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the distribution of the transverse (with respect to the virtual photon direction) hadron 
momentum p^x in the final state. (For simplicity, the dependence of these functions on 
the longitudinal momentum is not shown here.) In addition, the result must be averaged 
over the (unmeasured) transverse polarization direction 05 of the quark in the initial 
state: 

d0 5 / d 2 PT I d 2 k T hf q (p 2 T )( PT x S q ) z (25) 

x Ht(k 2 T )(k T x S' q )J(p T + k T - p h>T ). 

As Ph,r is the only vector on which the transverse momentum integral can depend, 
rotational invariance implies 

J d 2 p T J d 2 k T h{\p 2 T ) k T>i Hi(k 2 T ) p T>j 5(p T + k T - p h>T ) (26) 
with some suitable functions a{p 2 j-) anc ^ Kp1t)i yielding 



A(p h>T ) oc / d0 s Upl, T )S q ■ S' q + b{p 2 h)T ) (S q x p KT ) z (S' q x p h>T ) 1 . (27) 
Jo 1 J 

The transverse polarization of the active quark does not always flip, but events where it 

does not flip do not produce any asymmetry. For the determination of the asymmetry we 

thus focus on those events where the transverse spin component in the lepton scattering 

plane flips. Taking this plane to lie in the x-z plane, the transverse polarization vectors 

take on the form S q = (cos 0^, sin 0^) and S' q = (— cos 0s, sin 0s), yielding 

rd<PsS q -S' q =0 (28) 

Jo 

d0 5 (S q x p KT ) z (S' q x p h>T ) z = - {p 2 h Tx - p 2 hTy ) = -^ T co S 20, (29) 

where specifies the transverse direction of the outgoing hadron momentum Ph,r- In 
combination with (27) this implies a cos 20 asymmetry for the transverse momentum 
distribution of the produced hadrons 

Mpkt) oc p 2 hiT b{p 2 h)T ) cos 20 (30) 

which involves a convolution b(p\ T ) of the Boer- Mulders and Collins functions. We note 
here that a flip of the quark polarization direction relative to the normal of the lepton 
scattering plane is essential to provide a reference plane (the lepton scattering plane) for 
the transverse momentum of the outgoing hadrons and thus making the Boer-Mulders 
effect observable, just as it was essential for separating the Collins from the Sivers effect 
in SIDIS. We also note that the Boer-Mulders effect is a spin effect, even though no spin 
is being measured. 

Alternatively, one can also study the cos 20 asymmetry in the unpolarized Drell- 
Yan process, in which case is the azimuthal angle of the plane about the virtual 

photon axis with respect to the incident proton trajectory. 
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While ISI/FSI are obviously essential features for the Sivers and Boer-Mulders 
functions, they affect all twist-2 pj-dependent parton distributions. However, this does 
not spoil standard universality, as all other p^-dependent PDFs are (naive) T-even 
and thus the ISI in DY has the same effect as the FSI in SIDIS. Nevertheless, their 
Pt dependence is expected to differ from the 'intrinsic' px dependence that one would 
for example obtain theoretically through the evaluation of a matrix element of a local 
operator, as is usually done in lattice gauge theory (see section 4.4). 

3. 8. Recent semi-inclusive data with transverse target polarization 

Almost 20 years ago, surprisingly large single-spin asymmetries were observed in 
inclusive production of pions by beams of polarized protons and anti-protons [104]. 
(Recently this phenomenon was also observed at the higher energies of RHIC [105, 106], 
where the cross section for these reactions can be understood within the framework 
of perturbative QCD.) Various theoretical interpretations were proposed, including the 
Collins and Sivers effects. However, no way of disentangling these effects on those data 
could be found. A similar ambiguity plagued the interpretation of the first substantial 
single-spin asymmetries, with respect to longitudinal target polarization, observed in 
the semi-inclusive leptoproduction of mesons by the Hermes experiment [107]. When 
the target is polarized along the lepton beam axis, both Collins and Sivers mechanisms 
produce a common sinusoidal behavior in the azimuthal angle <fi of the detected hadron 
about the direction of the virtual photon, with respect to the lepton scattering plane. 
On the other hand, when the target polarization is orthogonal to the beam axis 
(more precisely, when it is orthogonal to the direction of the virtual photon), the 
azimuthal angle 05 of the axis of transverse target polarization resolves the effects of 
the two mechanisms into distinctive signatures: sin(0 — for the Sivers mechanism, 
and sin(</> + <ps) for the Collins mechanism [103]. In DIS, the Collins effect serves 
as a 'polarimeter' sensitive to transverse quark polarization to reveal its correlation 
with transverse target polarization, which is known as transversity. As described in 
section 3.4, the Collins mechanism depends on the influence of the struck quark's 
polarization on the transverse momentum kx acquired by the produced hadron in the 
fragmentation process, in particular on the component of fey that is orthogonal to the 
quark's transverse polarization. This 'measurement' of the struck quark's polarization 
occurs after its spin component in the lepton scattering plane has been flipped by the 
photo-absorption. The orientation of the lepton scattering plane is 'remembered' by the 
virtual photon because that is the plane of its linear polarization. Hence the Collins 
mechanism is sensitive to the orientation of the lepton scattering plane from which 
the azimuthal angles are measured, as indicated by the dependence on <fi + In 
contrast, as described in section 3.6, the Sivers effect can be understood in terms of 
retention by the struck quark of the Pt that it had in the target due to a correlation 
with transverse target polarization. This transverse momentum tends to be inherited 
by a forward hadron that may 'contain' this quark. Hence in either case, the hadron 



Spin-polarized high- energy scattering of charged leptons on nucleons 



44 



Ph,T is correlated with transverse target polarization. However, the different azimuthal 
dependences allow their experimental separation, and semi-inclusive DIS measurements 
with a transversely polarized target were eagerly anticipated. 

The first such data were reported a few years ago by Hermes [114]. The 
cross section asymmetry with respect to target polarization was extracted as a two- 



This simultaneous extraction of both contributions was shown by detailed Monte Carlo 
simulations to avoid significant cross-contamination, even when they have very different 
magnitudes in the context of a limited detector acceptance. These initial data on the 
proton target provided evidence for both Collins and Sivers signals. On the other hand, 
asymmetries measured by Compass using higher energy muons on a deuteron target 
are consistent with zero [115]. As illustrated in figures 17 and 21, the nonzero results 
for the proton are now firmly established by recent more precise preliminary data [108], 
while the most recent results for the deuteron remain consistent with zero [109]. For the 
proton target, the Collins azimuthal amplitude for n + is clearly positive and non-zero, 
while it is negative and at least as large for 7r~ . The latter finding was a surprise because 
of the contrast with the double-spin longitudinal asymmetries shown in figure 7. In that 
case asymmetries for pions of both charges are positive, even though Am is positive and 
Ad is negative, albeit smaller in magnitude. The reason there is that both charges are 
produced mostly from photoabsorption by u quarks [63], but not as dominantly for 7r~. 
Lattice calculations and models predict that transversity densities resemble the helicity 
densities, at least at large x where valence quarks dominate, and where the tt + /tt~ 
contrast is most apparent in figure 17. Hence efforts to find an explanation focused on 
the Collins fragmentation function. 

The explanation that emerged for the larger negative ir~ amplitudes is a substantial 
magnitude for the disfavoured Collins function describing e.g. the fragmentation of u 
quarks to n~ mesons, and with a sign opposite to that of the favoured function. Opposite 
signs could be expected in the light of the string model of fragmentation discussed in 
section 3.5. If a favoured pion forms as the string end created by the first break, 
a disfavoured pion from the next break will inherit transverse momentum from the 
first break in the opposite direction from that acquired by the first pion. Such a p Wt T 
anticorrelation between favoured and disfavoured pions is demonstrated by the Jetset 
simulation [116], which is based on a string fragmentation model. Hence any correlation 
between p n< T an d another kinematic or spin observable should have the opposite sign 
for favoured and disfavoured pions. 

This expectation for the Collins function has recently been confirmed by results 
for both favoured and disfavoured functions from a simultaneous fit of these DIS data 
together with information about the azimuthal structure of jets in e + e~ collisions from 



dimensional distribution in <f> versus <ps- The Collins azimuthal moment ( sm^+^s)) UT 
and Sivers moment ( sin(0 — (J) S ))^ T of the virtual-photon asymmetry were extracted 
in the fit 





= (sm((f) + 4> s )) UT sm((f) + 4> s ) + (sin(0 - <f> s )) UT san(<f> - fa) .(31) 
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Figure 17. Collins sinusoidal amplitude of the lepton-beam asymmetry with respect to 
transverse polarization of a proton [108] (deutcron [109]) target, for the semi-inclusive 
production by positrons (muons) of identified charged mesons as labelled in the upper 
(lower) panel, as a function of x, z or Ph,t (the latter appearing as Pt in the figure). 
The sign of the amplitudes is according to the Trento Conventions [110], although for 
the deuteron target the sign is reversed because of the extra value of ir added to the 
argument of the sin function. The error bars represent the statistical uncertainties, 
while for the proton target there is an additional common 8% scale uncertainty. The 
curves are fits to the data, as described in the text [111]. The error bands on these 
curves arise from the statistical uncertainties of the data, but the procedure chosen 
implies that they are several times larger than the conventional definition given by, 
e.g., the Particle Data Group [71]. On the other hand, uncertainties related to model 
assumptions are not estimated. The figure is taken from [111]. 
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Figure 18. Two alternative jet asymmetries extracted from Belle data [112] for 
e + e~ — > 2 jets, compared to a fit [111] of parameterized Collins functions. The figure 
is taken from [111]. 



the Belle collaboration [112]. Two alternative jet asymmetries extracted from the 
Belle data are shown in figure 18. The most recent such fit [111] is shown as the 
curves in figures 17 and 18. Such fits of data from these three experiments produced 
the first experimental determination of the transversity distributions of u and d quarks, 
the most recent example of which is shown in figure 19. At this stage, the quark sea is 
neglected. As is the case for the helicity distributions, transversity is found to be positive 
for u quarks and negative and smaller in magnitude for d quarks. For both flavours, 
the magnitudes are of order half of the Soffer positivity bound [85]. Especially for u 
quarks and for larger values of x (bearing in mind that there are no data for x > 0.3), 
the extracted density is smaller than existing model predictions [117-119]. This fit also 
yields values for the Collins fragmentation function, shown in figure 20. They confirm 
that indeed the disfavoured Collins function is opposite in sign and larger in magnitude 
than the favoured function. The initially surprising features of the Collins asymmetry 
data now appear to be well understood. A further confirmation of this recently appeared 
in the form of preliminary single-spin asymmetry data for unidentified hadrons from the 
proton by the Compass collaboration [120]. These data for the Collins asymmetry are 
consistent with predictions based on this fit of earlier data from Hermes on the proton 
and from Compass on the deuteron. 

The Sivers azimuthal amplitudes for production of both pions and kaons from the 
proton shown in the upper panels of figure 21 are clearly positive and nonzero for the 
positive charges. This constitutes the first evidence from leptoproduction for a nonzero 
T-odd parton distribution function. We note here that these Sivers amplitudes receive 
no contributions that are suppressed by 1/Q (twist-3) [91]. Since the production of 




Figure 19. The (red) curves within the error bands represent values of the transversity 
distributions at Q 2 = 2.4 GeV 2 for u (upper panels) and d (lower panels) quarks as a 
function of x and pr (labelled k± in the figure), from the fit [111] to the pion data shown 
in figure 17 together with e + e~ data from Belle [112]. The pr dependence in the 
right panels was chosen to be the same as that used for the unpolarized distributions, 
which was derived from fitting the azimuthal dependence of unpolarized SIDIS cross 
sections. The 'Torino notation' used in the figure for the transversity distributions is 
equivalent to the 'Amsterdam notation' used elsewhere in this paper. (This is not true 
for, e.g., TMDs [110].) The outermost (blue) curves represent the Soffer positivity 
bounds [85]. The error bands have the same meaning as in figure 17. The figure is 
taken from [111]. 



both 7r + and K + is dominated by u quarks, this positive value with the definition of 
azimuthal angles used here implies a negative value for the Sivers function of this 
flavour [110]. The amplitudes for the negatively charged mesons are consistent with zero. 
Their interpretation is complicated by competition between fragmentation of quarks of 
at least two flavours. Remarkably, for the proton target, the K + amplitudes are on 
average much larger than those for ir + , even though one might suppose that they are 
similarly dominated by photoabsorption by u quarks. The amplitudes for the deuteron 
target measured by Compass are all consistent with zero (as are their preliminary Sivers 
amplitudes on the proton for unidentified hadrons [120]). 

The curves on these figures from [113] represent an 11-parameter fit to data for 
pions and kaons from Hermes on the proton and from COMPASS on the deuteron, 
extracting representations of the Sivers functions for all six quark flavours. The shapes 
of all four sea flavours u, d, s and s were constrained to be the same. The success 
in explaining the large K + amplitudes for the proton depends on the use of a recent 
fit of unpolarized fragmentation functions [121], which differs from all previous such 
fits in that Df + (z) > D* + (z) over the whole z range. This feature implies that K + 
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Figure 20. The (red) curves within the error bands represent values of the favoured 
and disfavoured Collins fragmentation functions at Q 2 = 2.4 GeV 2 as functions of 
z or kx (the latter appearing in the figure as p±), from the fit [111] to the pion 
data shown in figure 17, together with e + e~ data from Belle [112]. In the left 
panels, the /c^-integrated functions are shown as ratios to twice the corresponding 
unpolarizcd fragmentation functions, while in the right panels, the functions arc 
evaluated at fixed values of z. The Collins function is represented in the Torino 
notation, which is related to the Amsterdam notation used elsewhere in this paper by 
A N D h / q ]{z,p 2 hT ) = 2 ^/^ ffli~ g ( z )Pfa,r) [HO]. The outermost (blue) curves represent 
the positivity limits. The error bands have the same meaning as in figure 17. The 
figure is taken from [111]. 



production need not be dominated by u quarks, allowing s quark distributions to be 
probed. 

Values of the Sivers functions resulting from the fit are shown in figure 22. The 
Sivers function for d quarks is found to be at least as large in magnitude as that for u 
quarks. Their signs agree with expectations based on an attractive final-state interaction 
between the struck quark and spectator diquark (see section 3.4). Such a relatively large 
magnitude of the Sivers function for d quarks could be explained by a similar trend in 
the known relative contributions of u and d quarks to the anomalous magnetic moment 
of the proton, as explained in section 4.3. Hence the data appear to support a profound 
connection between static and partonic properties of the proton. Furthermore, there is 
evidence for a significant role of sea quarks. In considering the statistical significance of 
this finding, it is useful to bear in mind that the method of assigning the error bands 
results in values much larger than those that would correspond to the definition of 
e.g., the Particle Data Group [71]. On the other hand, uncertainties related to model 
assumptions are not estimated. 

We now look forward to measurements of single-spin asymmetries for the Drell-Yan 
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Figure 21. Sivers sinusoidal amplitude of the asymmetry with respect to transverse 
polarization of a proton [108] (deuteron [109]) target, for the semi-inclusive production 
by positrons (muons) of identified charged mesons as labelled in the upper (lower) 
panels, as a function of x, z or Ph,T (the latter appearing as Pt in the figure). The 
sign of the amplitudes is according to the Trento Conventions [110]. The error bars 
represent the statistical uncertainties, while for the proton target there is an additional 
common 8% scale uncertainty. The curves are fits to the data, as described in the 
text [113]. The error bands on these curves arise from the statistical uncertainties of 
the data, but the procedure chosen implies that they are several times larger than 
the conventional definition given by, e.g., the Particle Data Group [71]. On the other 
hand, uncertainties related to model assumptions are not estimated. The figure is 
taken from [113]. 



process, from which values of the Sivers function can be extracted that are predicted 
on the basis of fundamental properties of QCD to be equal in magnitude but opposite 
in sign to these from DIS, as in (23). This will test our basic understanding of T-odd 
parton distribution functions within the context of QCD. Thus continues the exciting 
story of a feature of proton structure that only a few years ago was (incorrectly) asserted 
to be forbidden by fundamental time reversal invariance [98]. 

4. Towards a 3-dimensional picture of the nucleon 

4-1. Generalized parton distributions 

In section 1 it was explained that factorization can be applied to the theoretical 
interpretation of DIS cross sections. Using the optical theorem, the inclusive DIS 
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Figure 22. The left [right] panels show as continuous (red) curves within the error 
bands values of the first px [x] moment of the Sivers function as functions of x [pt, 
appearing in the figure as k±], from a fit [113] to the data shown in figure 21. The 
Sivers function is represented in the Torino notation, which in this case is opposite in 
sign and equal in magnitude to that in the Amsterdam notation used elsewhere in this 
paper. The error bands have the same meaning as in figure 21. The figure is taken 
from [113]. 



cross section can be related to the imaginary part of the forward Compton amplitude, 
which can also be factorized into a 'hard part' describing the vertices of the active 
quark with the photons and the quark propagating between these two vertices, and 
a soft part describing the quark correlation function in the hadron (see figure 23a). 




Figure 23. Comparison between a) the Compton amplitude with two photon vertices 
and a quark propagating between them and b) a form factor with just one photon 
vertex. 



A similar factorization of the scattering amplitude occurs in Deeply Virtual Compton 
Scattering (DVCS) [122], 7*p — > 7p, where 7* is a space-like photon with virtuality 
q 2 = —Q 2 < 0. As in DIS, the process is best interpreted in a frame in which the 
nucleon has infinite momentum opposite to that of the virtual photon. In the Bjorken 
limit, and for a momentum transfer t to the proton much less than Q 2 , the DVCS 
amplitude factorizes into a convolution of the Compton amplitude for scattering off a 
quark constituting the hard part, and a quark correlation function, constituting the soft 
part of the amplitude. After factorization has been applied to the DVCS amplitude, 
its matrix element has much in common with the matrix element describing the form 
factor (figure 23b). Indeed, the soft part entering the DVCS amplitude becomes identical 
to the soft part entering the form factor. However, the DVCS process provides more 
'surgical' access. In the case of the DVCS amplitude, the hard part has a dependence on 
the quark momentum fraction appearing in the quark propagator, which is not present 
in the form factor. The additional kinematic degree of freedom provided by the hard 
part in the DVCS process provides access to the dependence of the soft part on the 
momentum fraction x of the active quark. This motivates the concept of Generalized 
Parton Distributions (GPDs) as a description of the soft part that enters the DVCS 
amplitude [123]. What is gained by introducing this concept is that universality as 
introduced in section 1 applies with an extended scope, now also encompassing form 
factors and hard exclusive processes — the same GPDs that describe elastic form factors 
as well as the soft part of DVCS should also appear in other hard exclusive reactions. In 
each of these reactions, they appear convoluted with the appropriate coefficient function 
representing the hard part. In section 4.5ff, it is discussed how and to what extent 
exclusive reactions can be used to constrain GPDs experimentally. 

In elastic scattering, the photon can couple to any of the quark lines connecting 
the two 'soft blobs' in figure 23b, leaving the nucleon intact. The form factor is the 
coherent superposition of such diagrams for all the quarks of all flavors in the nucleon. 
The experimental decomposition of the nucleon form factor into the contributions from 
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various quark flavors has been partially accomplished through the comparison of proton 
and neutron form factors using isospin symmetry, and the interpretation of parity- 
violating electron scattering. Two flavor combinations (2u — d — s) and 2{d — u — s) 
are obtained from the proton and neutron electromagnetic form factors using charge 
symmetry. A further combination is accessible in parity violating electron scattering, 
but the corresponding constraints are not very strong yet [124]. 

However, what is still needed for a more comprehensive understanding of the 
microscopic physics that underlies form factors is their decomposition with respect 
to the momentum fraction of the quark that absorbs the photon. Generalized 
parton distributions (GPDs) embody this information. By definition, they provide a 
decomposition of form factors evaluated at a given value of the invariant momentum 
transfer t = [pjin —Pini) 2 on the target, with respect to the average momentum fraction 
x = \ (xi n i + Xfi n ) of the active quark of flavor /: 

J dxHf{x,U) = F{(t) f dxEf(x,U) = Fl® (32) 

J dx H f (x, £, t) = G f A {t) J dx E f (x, f , t) = G f P (t), (33) 

where Xi n i and Xfi n are the initial and final momentum fractions of the active quark 
with respect to the average momentum | (p ini + Pfi n ), and the dimensionless 'skewness' 
£ = \{.Xini — Xfin) represents their difference. Thus the linear combinations x ±£ 
represent the momentum fraction carried by the quark before/after the hard interactions. 
The functions F({t), F((t), G f A {t), and G f P (t) are the contributions to the Dirac, Pauli, 
axial, and pseudoscalar form factors from quark flavor /, respectively, and the functions 
, E*, H* , and E* are the quark GPDs providing their momentum decomposition. 
One can thus imagine the form factor being the result of adding (or integrating over) 
the contributions from all possible quark momentum fractions x (here x is not to be 
identified with the Bjorken variable x B defined in section 1). 

Just like the ordinary (forward) PDFs, GPDs are subject to higher-order corrections 
in a s and are thus also dependent on the scale Q 2 (and implicitly also on the factorization 
scheme used) at which they are probed. Nevertheless, as F( (t) and F^if) are defined 
as matrix elements of a conserved current, they are Q 2 independent and hence the Q 2 
dependence must disappear in the lowest moment of the corresponding GPDs (32). 

Decomposing the form factors with respect to the quark momentum fraction x 
involves the infinite-momentum direction. Hence it can make a difference whether the 
momentum transfer to the target is parallel, or perpendicular to this momentum. The 
GPDs must therefore depend on an additional variable, here called £. Its relationship 
to t/M 2 characterizes the direction of the momentum transfer relative to the infinite 
momentum direction. The case £ = corresponds to the case where those two momenta 
are perpendicular to each other (which is not directly experimentally accessible in 
DVCS), while nonzero £ corresponds to the situation where the momentum transfer 
to the target has a component parallel to the virtual photon direction. Very little is 
known about the full x vs. £ dependence of GPDs, except that they might possibly be 
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Figure 24. Model calculation for H q (x, £, 0) from [125]. The thick (red) line represents 
the forward limit £ = 0. Accessible to the DVCS process is the 'outer' or 'DGLAP' 
region |x| > £ that represents quarks for x > £ and antiquarks for x < — £. In the 
'inner' or 'ERBL' region |x| < £, the quark GPD behaves like a meson distribution 
amplitude in deeply virtual meson production. The figure is taken from [125]. 



quite complex, as indicated for example by the model calculation shown in figure 24. 
In (32) and (33), the integration over x removes all reference to the infinite momentum 
direction, with respect to which x and £ are defined. Lorentz invariance thus implies 
that the £ dependence disappears in these integrals. 

One constraint on quark GPDs comes from the fact that in the limit of vanishing 
t and £, the nucleon helicity non-flip quark GPDs must reduce to the respective PDFs: 

H f (x,0,0) = q f (x), H f (x, 0, 0) = Aq f (x) for x > (34) 
H f (x,0,0) = -q f {-x), H f (x,0,0) = Aq f (-x) for x < 0. (35) 

GPDs became the focus of intense interest after the discovery of the 'Ji-relation', 
connecting a second moment of GPDs with the total angular momentum of quarks with 
flavor / or gluons, which is the respective contribution to the nucleon spin | [19]: 

/ 

where 

i r 1 

J f = -gm y dxx [Hf(x,U) + E f (x,U)] , (37) 

J G = \ Jim / dx [H G (x, £, t) + E G (x, £, t)} . (38) 
^ Jo 

Note that no ^-dependence remains in (37) and (38) after integration over x, although 
the two unpolarized quark GPDs H* and E$ have to be evaluated at some common 
value of £. Unlike the lowest moments (32), these second moments and hence Jf t c 
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depend on Q 2 (as well as order in ct s and factorization scheme), but this dependence is 
usually suppressed for simplicity. 

Equation (37) is usually applied to the z component of the quark total angular 
momentum in a longitudinally polarized nucleon. However, in principle (37) can be 
applied for any polarization direction of the nucleon target. Nevertheless, (37) finds its 
most important application in the helicity sum rule through 

J f = ±(Aq f + Aq f )+Lf, (39) 

since in that case at least the quark helicity contribution Aqf + Aqf in Jf has a partonic 
interpretation. (An attempt to associate a partonic interpretation of (37) in the case of 
transversely polarized nucleons can be found in [126].) 

The gauge invariant orbital angular momentum Lf of quarks is defined through the 
expectation value of ipff x (p~ gAj ipf. No process has been identified to measure Lf 
directly, but it can be evaluated indirectly through (39) as 

L f = J f -\ (Aq f + Aq f ) . (40) 

The issues arising in a decomposition of the nucleon spin are further discussed in section 
5. 



4-2. General properties of GPDs 

GPDs must satisfy several relations arising from fundamental symmetries or 
conservation laws: 

• Invariance under time reversal (for both quark and gluon GPDs): 

F(x,-£ 1 t) = F(x,£,t), (41) 

• Polynomiality property of GPDs [127] as a consequence of Lorentz invariance: 

/l n-1 
dxx n ^Hf(x,^t) = ^(20X fc W + mod(n-l,2)(20"^(t), (42) 

1 fc=0 

even 

/l n-1 
&xx n - l Ef{x,t„t) = X)(20*<*(*) - mod(n - l,2)(2£) B C2(t). (43) 
-i 



fc=0 
even 



Polynomiality means that the entire £ dependence of the n'th x-moment (n > 0) 
of a GPD is given by a polynomial in £ of order at most n. Polynomiality is 
satisfied, e.g., by the representation of GPDs in terms of double distributions (DDs) 
F } DD ((3,a,t) [128-130], which can be considered as spectral functions of nonforward 
matrix elements. In terms of a DD, a generic quark GPD F can always be written 
as: 

•l pi- 



nl ri-\P\ 

F f (x,^,t) = d/3 da 6(x - (3 - fa) F[, D (J3, a, t) 

J-l J-l+\8\ 



+ 9(£-\x\)Df(^t). (44) 
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The step function 8(() is zero for ( < and unity for £ > 0. Thus the so- 
called 'D-term' D* [127] contributes only in the 'ERBL' region — £ < x < £, 
where quark GPDs have the character of distribution amplitudes for the creation 
of a quark/antiquark pair. It does not contribute in the complementary so- 
called 'DGLAP' region \x\ > £, where quark GPDs describe the emission and 
reabsorption of an (anti-)quark in the infinite momentum frame, thereby having 
properties analogous to the familiar (anti-)quark distribution functions. The D- 
term provides a convenient means of representing this profound difference in GPD 
properties between the two regions, which would otherwise require singularities in 
the DD having a severity beyond represention by delta functions or their derivatives. 
Another perspective is that the maximum power of £ provided by the DD itself is 
one less than that allowed by Lorentz invariance. Expressions analogous to (44) 
hold for the polarised quark GPDs H and E, for which the D-term vanishes. The 
mapping between GPDs and DDs is not one to one; each GPD can be represented 
by any of an infinite family of DDs. The quantity C? n (t) in (42,43) is the n'th 
x-moment of -D^(|, t), which has the same magnitude for the quark GPDs H* and 
E* but enters with the opposite sign. Hence the second moment of the sum of the 
unpolarized quark GPDs Hf(x,£,t) + Ef(x,£,t) ('Ji integral' (37)) is independent 
of£. 

Several sum rules apply at £ = and t — [131]. Suppressing the kinematic arguments: 
• Quark and gluon GPDs are related by the conservation of momentum, 



1= / dxxJ^H 1 + I dxH G , 
J-i f Jo 



I 

and of the proton's angular momentum, 



- = - [ dxxJ2(H f + E f ) + - [ dx(H G + E G ) 
2 2 J_ x 2 J 



(45) 



(46) 



Combining them leads to: 



= y dxx^E f + J dxE c 



(47) 

/ 

For each quark flavour, the GPD E* is constrained by the contribution Kf of flavour 
/ to the proton anomalous magnetic moment k: 

dxE f = F}(t = 0) = K f (48) 

-l 

As quark GPDs are not generally even or odd in x, it is often useful [131] to consider 
the GPD combinations 

H^\x, C, t) = Hf(x, f , t) T H*(-x, £, t) (49) 
H^\x, 1 1) ee H f (x, e, t) ± H f (-x, £, t), (50) 
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and their analogues for and E> '. Here the upper (lower) sign corresponds to the 
exchange of a C = +1 (C = —1) object in the t channel, and denotes the 'singlet' 
('non-singlet' [or valence]) combination. These names are taken from the corresponding 
forward limits given by: 

# /(±) (x,0,0) = q f (x) ±q f (x) (51) 

H fi±) (x,0,0) = Aq f (x) ± Aq f (x), (52) 

for x > 0. While the singlet (C-even) combination of quark GPDs is probed in DVCS 
(see (68)) or exclusive vector meson production, the non-singlet (C-odd) one is selected 
in pseudoscalar meson production. 

Gluons are their own antiparticles so that gluon GPDs are either even in x (H G , E ) 
or odd (H G ,E G ). Their forward limits involve an additional factor of x, so that (for 
x > 0): 

H G (x,0,0) = xG(x), H G (x,0,0) = xAG(x). (53) 
4-3. Impact-parameter- dependent parton distributions 

A particularly simple physical interpretation for GPDs can be developed in the limiting 
case £ = 0, where the quark momentum fraction is the same in the initial and final state. 
So far we have introduced GPDs as momentum-disected form factors. Hence a GPD 
evaluated at a certain value of x can be considered to be the form factor for only those 
quarks in the nucleon that carry momentum fraction x. In analogy with nonrelativistic 
form factors, whose Fourier transform provides information on the charge distribution 
in position space, one would thus expect that the Fourier transform of the t dependence 
of GPDs for fixed x provides information on the position space distribution of quarks 
carrying a certain momentum fraction x. Indeed, it can be shown that the distribution 
of partons in impact parameter space can be obtained from GPDs via [20]: 

b±) = J Wf e ~ tA± ' b±Hf{x ' °' ~ Ai) ' (54) 
where — A.\ = t for £ = 0. Analogously, the parton helicity distribution in a 
longitudinally polarized target is obtained by replacing — > Aqf and H* — > H*. Note 
that the Heisenberg uncertainty principle forbids a simultaneous exact measurement of 
x and the longitudinal (i.e. along the direction of the virtual photon) position of the 
quarks [22]. Therefore, only the impact parameter b±, i.e. the position transverse to the 
nucleon's momentum can be determined simultaneously with x. For the determination 
of the b± distribution, the momentum transfer should be purely perpendicular, which 
implies the restriction to £ = 0. The reference point for the impact parameter b± is 
the transverse center of longitudinal momentum R±, which plays a role similar to the 
center of mass in a nonrelativistic system. In analogy to that case, it is obtained as an 
average of the positions r_|_ ti of all the partons, except that the weight factor is not the 
mass fraction, but the momentum fraction Xi carried by each parton [132]: 




Figure 25. Qualitative shape of q(x, b±_). 



For x — > 1, the active quark dominates the sum in (55) and q(x, 6j_) must necessarily 
become narrow in this limit (see figure 25). 

In analogy to the inclusive PDFs that are measured in DIS, given that we have 
restricted ourselves to £ = 0, these Impact-Parameter-Dependent Parton Distributions 
(IPD-PDFs or simply IPDs) have a probabilistic interpretation and hence they satisfy 
a variety of inequalities, such as [133,134]: 

Qf(x,bx) > \Aq f (x,b ± )\ > for x > 0, (56) 
-qf(x,b±)> \Aq f (x, b±)\ > for x < 0. 

Remarkably, (54) is not limited to the nonrelativistic domain of \t\ < m 2 , in contrast 
to the interpretation of form factors as i/iree-dimensional Fourier transforms of charge 
distributions. Its only limitation is \t\ <C Q 2 as 1/Q represents the spatial resolution of 
the virtual photon. Since the electromagnetic form factors are obtained as integrals of 
GPDs (see (32) and (33)), this feature of (54) implies that in the infinite momentum 
frame, the identification of the two-dimensional Fourier transform of form factors with 
the charge distribution in the transverse plane is also not limited to \t\ < m 2 . 

For a longitudinally polarized or unpolarized target, the distribution of quarks in 
the transverse plane is axially symmetric. However, when the target is transversely 
polarized, the polarization direction of the quark in the transverse plane breaks this 
axial symmetry. For example, for a target polarized in the x direction, the quark 
distribution qf(x, b±) is deformed along the y direction and the sign and magnitude of 
this deformation are embodied in the GPD E*{x, 0, — A^_) [135]: 

q f (x, b ± ) = (57) 
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Such deformations can be interpreted semi-classically as an enhanced probability of 
photoabsorption occuring on the left or right of the nucleon as seen by the incident 
photon, when the nucleon spin is pointing e.g. up. As was explained in section 3.6, this 
suggests an intimate relationship between these deformations, parton orbital angular 
momentum, and the Sivers effect. 

Experimentally, only the x-integral of the GPD E*(x, 0, — Aj_), i.e. the form factor 
F2(— Aj_), is known (see (32)). Any additional constraints come only from lattice 
calculations, which will be discussed in the following section, as well as model studies. 
As a result, the details of the above deformation are not very well known. In the 
considered case of a transversely polarized target, this deformation implies the existence 
of a transverse flavor dipole moment. The integrated deformation can be related model- 
independently to the anomalous magnetic moment of the nucleon. For example, when 
the target is polarized in the x (transverse) direction, the average deformation T>f of 
the distribution for quarks of flavor / in the y direction reads [135]: 

V f = Jjx J d 2 b ± q f (x, b ± ) \ = ^- (58) 

The Kf/ p are the contributions from quark flavor / to the anomalous magnetic moment of 
the proton with the quark charge factored out: n p = \n u /p — \>^d/p + ■ ■ ■ ~ 1-79. A flavor 
decomposition of the anomalous magnetic moment, using n n pa —1.91, charge symmetry, 
and neglecting the small contribution from s quarks yields k u / p ~ 1.67 and Kd/ P ~ —2.03, 
which translates into average deformations of T> u pa +0.16 fm and T>d pa —0.20 fm. 

Such a significant transverse deformation of quark distributions should have 
observable consequences, e.g. in SIDIS on a transversely polarized target. Armed with 
the above quantitative information, the semiclassical arguments of section 3.6 about the 
transverse deformation of PDFs as a cause of single-spin asymmetries can be made more 
specific. When a virtual photon collides with a proton that has transverse spin down, it 
'sees' more u quarks on the right side of the proton and more d quarks on the left side, 
with the signs being determined by the signs of K q / P (see figure 15). Since the final-state 
interaction acting on the escaping quark is on average expected to be attractive, more 
ejected u quarks are expected to be deflected to the left than to the right. For d quarks, 
the opposite is predicted. The predicted negative sign of the Sivers function (22) for 
u quarks fffl p [100], has been confirmed by a recent analysis of SIDIS data for pion 
production on the proton [108,114,115]. Here the experimental result for f^ p is less 
clear as its contribution to pion production is strongly suppressed. The asymmetry on 
the deuteron is consistent with zero [136,137]. A small Sivers asymmetry could arise 
from cancellation between the contributions from the proton and neutron. Using charge 
symmetry this implies ~ — /{j^ p , which is consistent with the prediction based 

on the transverse deformation of IPDs. However, it is also possible that at very small 
x, both and f^ p are small. 
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4-4- QCD simulations on a space-time lattice 

Ideally, in order to describe properties and interactions of hadrons theoretically, hadronic 
matrix elements should be evaluated in terms of QCD wave functions. However, this is 
not possible as long as wave functions in QCD cannot be calculated from first principles. 
QCD calculations on a Euclidean lattice (also known as Lattice QCD or Lattice Gauge 
Theory) are currently the most promising tool to overcome this problem. Presently there 
are a large number of workers active in this field — resulting in many types and flavors 
of lattice QCD. The proceedings of the yearly lattice symposia (currently published in 
PoS) provide a good overview of the state of the field; specifically concerning hadron 
structure, the two most recent summaries were reported in [138,139]. 

The inevitable truncations for numerical computations in QCD are less problematic 
on a space-time lattice, where space rather than momentum is discretized. After 
discretizing space-time on a 4-dimensional grid, approximate solutions for hadron 
masses and matrix elements are obtained numerically using Monte-Carlo techniques, 
because the latter are capable of dealing with the truly astronomical numbers of 
degrees of freedom (Contemporary lattices typically have about 32 3 x 64 discretization 
points. Including spin and color degrees of freedom, this implies a total of more than 
> 10 7 degrees of freedom for the entire lattice) that are necessary for an adequate 
approximation to the continuum. These Monte-Carlo techniques are used to evaluate 
a so-called 'path integral', which sidesteps a direct computation of wave functions. For 
example, a pion form factor can be computed without having to know what a pion 
is, only its quantum numbers [140-143]. The evaluation of these path integrals would 
be impractical working with a real time variable, as amplitudes in quantum mechanics 
oscillate as a function of time and would thus lead to poor statistical precision. These 
oscillations are avoided by introducing an imaginary time t — > —ir, which renders the 
rime (sic) evolution operator real and positive: exp(-iHt) — > exp(-Hr), and thus 
suitable for Monte Carlo sampling. The variable r is often referred to as 'Euclidean time', 
as the Minkowskian metric becomes Euclidean when expressed in terms of (r, r) . This 
last step requires an analytical continuation to translate the results into a Minkowskian 
space-time. 

It is impossible to apply lattice QCD to directly calculate DIS cross sections, since 
this would require a summation over all possible final states. Instead one can use the 
optical theorem and evaluate the imaginary part of the forward Compton amplitude. 
This requires the calculation of a correlation function (^P ip^(§)'~f + ipf(z~) Pj between 
nucleon states, where xjjf is the quark field operator. Here the light-cone variables are 
defined as z ± = (z° ± z 3 ), where the presence of z° makes the correlation function 
non-local in time. Such a space-time structure arises since the Compton amplitude 
involves two vertices in figure 23a coupling the nucleon to the virtual photon, separated 
in a light-like direction, as the struck quark moves with nearly the speed of light. This 
causes additional difficulties in the translation from imaginary rime to real time. This 
problem restricts lattice gauge theory calculations of PDFs to only their Mellin moments. 
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Such moments are accessible because they can be represented by matrix elements of 
local (i.e. containing products of field operators and their derivatives evaluated at the 
same space-time point) operators for which this analytical continuation is trivial. For 
PDFs, which can apparently be parameterized by relatively simple functions, only a 
few moments may be sufficient to constrain these functions, except very close to the 
endpoints x — ► 0,1. In higher moments the region x — ► is suppressed and they thus 
only indirectly help constrain this region. Even though higher moments emphasize the 
x — > 1 region, PDFs are very small for very large x, and only very high moments yield 
useful constraints on the x — > 1 behavior of PDFs. 

Parameterizing GPDs is significantly more difficult as there are two longitudinal 
variables instead of one, and the t-dependence adds yet another variable. Furthermore, 
it is not yet known if GPDs can be described by similarly simple functions as PDFs 
and hence it is not obvious how many moments will be needed to generate a reliable 
parameterization for GPDs. Indeed, both phenomenological model calculations (see 
figure 24) as well as perturbative QCD evolution [144] suggest a 'kink', i.e. a 
discontinuous first derivative, for GPDs at x — ±£, whose adequate parameterization 
may therefore demand a larger number of moments than was the case for PDFs. 

Remarkable progress in lattice calculations has been made over the last three 
decades. One source of systematic uncertainties is the coarseness (spacing) of the 
lattice. Improved algorithms now permit improved precision at the same lattice spacing 
[145-149]. Another source of uncertainties arises from the finite size of the lattice 
(physical volume). Ideally, it should be large enough to contain the hadron state 
completely. Both because improved algorithms allow larger lattice spacings at the same 
precision and because increased computing power allows more discretization points, 
large lattices with several fm diameter are now possible. In the past, the omission 
of sea quarks in 'quenched calculations' gave rise to additional uncertainties. Steady 
increases in computational resources and improved algorithms make it now possible 
to also include sea quarks ('full' lattice QCD) and unquenched calculations are now 
routine. Nevertheless, calculating observables that require inserting operators into 
disconnected quark loops is at present still prohibitively expensive in most cases; some 
exceptions which are tractable have been studied, e.g., in [150]. Usually, contributions 
from disconnected diagrams are omitted, resulting in an uncontrolled approximation 
that affects most flavour- singlet observables. Lattice calculations of flavour non-singlet 
observables are not affected by this omission and hence are much more reliable. For the 
Dirac form factor, the contribution from disconnected quark loops has been estimated to 
be small [151-154], consistent with the small influence from strange quarks on the charge 
radius of the nucleon [155,156]. For most other observables, it is not clear whether 
the disconnected quark loop contribution is small, but steady advances in computer 
technology will allow including those contributions in the near future. 

The masses of light quarks (/ = u, d) are of the order of a few MeV. Calculations are 
presently restricted to much higher values, corresponding to pion masses of m, « 250 



Spin-polarized high- energy scattering of charged leptons on nucleons 



61 



MeV or even higher. This limitation is due not only to the increase of statistical 
fluctuations in the Monte Carlo algorithms for small quark masses but also to the 
increase in the required overall lattice volume: a proper description of the virtual 
pion fluctuations around a hadron demands lattices that are several pion Compton 
wavelengths in diameter. As a result of the restriction to performing calculations in 
an 'unphysical universe' in which the light quarks are too heavy, another significant 
source of uncertainty arises from the extrapolation to small quark masses and the 
difficulties involved in restoring the chiral symmetry (17) that was broken by an 
expedient necessary in discretizing the Dirac equation. Recently, new methods for 
approximating fermion fields on a lattice have been applied which provide significantly 
better approximations of chiral symmetry [157-159]. The latter advances now allow 
more meaningful calculations at smaller quark masses, where chiral symmetry becomes 
increasingly important. Fortunately, it is not necessary to perform lattice calculations 
at quark masses for u and d quarks of a few MeV. Chiral Perturbation Theory (xPT) is 
an effective field theory for QCD at small quark masses and momentum transfers. It can 
be used to extrapolate results from larger quark masses down to these values. However, 
a reliable extrapolation first requires lattice simulations for quark masses where %PT 
can be trusted, a matter of controversy. 

Lattice QCD calculations have now been used to evaluate a large variety of 
observables, and due to lack of space we present only a few selected cases. Figure 
26 shows lattice results for the nucleon's axial charge [160], which is related to the 
spin content of the nucleon by Bjorken's sum rule (8). These results illustrate not only 
the wide range of quark masses that is now accessible, but also the accuracy that can 
be achieved when extrapolating results to the physical value m n = 139 MeV. For the 
tensor charge one finds after extrapolation to the chiral limit 5u — 5u = 0.857 ± 0.013 
and 5d-5d = 0.212 ± 0.005 [164]. 

Using lattice QCD calculations to determine the second moment of GPDs entering 
the Ji-relation (37) yields J u = 0.214 ± 0.027 and J d = -0.001 ± 0.027 [165, 166]. 
Combining these results with lattice calculations for Aqf + Aq~f one thus finds 

jlatt. = jlatt. _ _ ( Au latt. + fajlatt.j = _Q.195 ± 0.044 (59) 
jlatt. = jlatt. _ ( Ad latt. + A jlatt.^ = Q 200 ± 0.044 (60) 

at the physical pion mass where the uncertainties are only statistical since some of the 
systematic uncertainties are difficult to quantify. In these calculations, disconnected 
quark loops have been omitted in the calculation of the matrix elements. The perhaps 
surprising result that the net orbital angular momentum carried by the quarks is close 
to zero may thus be subject to large corrections. However, the isovector combination 
jlatt. _ jjatt. _q 4 j s no t a ff ec ted by the omission of disconnected quark loops, and 
therefore more reliable. Note that in all quark models, where the orbital angular 

% Since the quark mass itself is not an observable, results are usually presented as a function of the 
pion mass, which strongly depends on the masses of the quarks. 
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Figure 26. Comparison of all presently existing unqucnched lattice QCD calculations 
of the nucleon's axial charge g A (LHPC/MILC [160], LHPC/SESAM [161], RBCK 
[162], QCDSF/UKQCD [163]). The solid/dotted line and error band denote a chiral 
perturbation theory fit to the red squares only. The error bars reflect only statistical 
uncertainties and thus the differences between results from different collaborations 
gives an indication of the systematic uncertainty in these calculations. The figure is 
taken from [160]. 

momentum arises from the lower component of a confined relativistic quark wave 
function, such as the MIT bag model, u quarks have a positive orbital angular 
momentum, while for d quarks the orbital angular momentum is negative, yielding 
L u — Ld > 0. Moreover, the sign of the Sivers function, both in SIDIS and in purely 
hadronic reactions strongly suggests L u > L d . QCD evolution of the quark orbital 
angular momentum from very low Q 2 , where the quark models are appropriate, to 
the lattice scale of about 4 GeV 2 has been suggested as a possible solution for this 
discrepancy [167]. However, in considering this apparent discrepancy, one should bear in 
mind that there are different possibilities for defining what is orbital angular momentum, 
which is discussed in detail in section 5.1. 

Figure 27 shows results for the lowest moment J dx q(x,b±) of IPDs. The 
magnitude of the deformation in these panels (i.e. the deviation from axial symmetry) 
is a measure for the correlation between the nucleon spin and the quark orbital angular 
momentum (left) and between the quark spin and quark orbital angular momentum 
(right). Notice that the average density on the left is twice as large as on the right 
since both quark polarizations are summed over on the left. The figures clearly indicate 
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Figure 27. Lowest moment of the IPDs for vmpolarized quarks in a transversely 
polarized nucleon (left) and transversely polarized quarks in an unpolarizcd nucleon 
(right) for u (top) and d (bottom) quarks. In the inserts, the quark spins (inner arrows) 
and nucleon spins (outer arrows) are oriented in the transverse plane as indicated. The 
figure is taken from [168]. 



that the correlation between quark angular momentum and the nucleon spin (left) has 
opposite signs for the two quark flavors, while they have the same sign for the correlation 
with the quark spin (right). Applying the chromo dynamic lensing mechanism (section 
4.3), this is consistent with opposite sign Sivers functions /{^ (see section 3.6) for u and 
d quarks. Using the same reasoning for transversely polarized quarks in an unpolarized 
target, this leads to same-sign Boer-Mulders functions hi q (see section 3.7) and the sign 
is the same as the sign of the Sivers function for u quarks (negative). The correlations 



hf q > f : 



e±q 
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on the right also appear somewhat stronger than on the left, suggesting 

Lattice QCD calculations play an important role in complementing both 
experiments and phenomenological models. For selected hadron structure observables, 
lattice calculations are reaching a precision that will allow the comparison of 
QCD predictions to experimental data. Lattice calculations also complement 
experiments because there are many interesting observables which are difficult to probe 
experimentally, but which are easily accessible on a lattice. Examples include most 
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chirally odd distribution functions and form factors. Nevertheless, calculations on the 
lattice will never be able to fully replace experiments in the exploration of hadron 
structure. For example, no local manifestly gauge invariant operator + exists for AG, 
which is therefore inaccessible in lattice QCD calculations. Simulation of scattering 
reactions and investigation of unstable particles is enormously difficult. Already the 
determination of low energy scattering phase shifts requires a sequence of calculations 
on lattices of various volumes [169-171], thus multiplying the computational time. Only 
those resonances (unstable particles) the mass of which is close to their decay threshold 
are accessible by this technique. Furthermore, the necessary use of Euclidean rime in 
lattice calculations makes those observables inaccessible that require calculations with a 
real time variable. For those high-energy scattering processes where initial and/or final 
state interactions are essential, such as SIDIS, no useful algorithm has been developed 
yet and the interpretation of experiments still has to rely on phenomenological models. 
Even after factorization has been applied in SIDIS, calculating the final-state interaction 
phase that is crucial for TMDs would again require a real time variable. This is the 
main reason why no one has found a direct way for calculating the Sivers (22) or Boer- 
Mulders (24) functions on the lattice. However, ISI/FSI affect the pr dependence of all 
PDFs and since the lattice presently can only calculate intrinsic distributions (without 
ISI/FSI) a direct comparison to experimentally measured p^-dependent distributions 
(which include ISI/FSI) thus requires modeling of those effects. Likewise, no algorithm 
exists for calculating fragmentation functions on the lattice. 

As a result of these limitations, there has also been a significant synergy between 
lattice calculations and phenomenological models in recent years. Even though lattice 
calculations cannot be used for direct calculations in dynamical processes, they can 
play a useful role in providing phenomenological models with input parameters, thus 
effectively increasing the predictive power of these models. The above-mentioned 
predictions for the Boer-Mulders functions is an excellent example for such a synergy 
between lattice and phenomenology: in the 'chromodynamic lensing' model for SSAs 
presented in section 4.3, IPDs are needed as input. For the Boer-Mulders functions, the 
relevant GPDs are chirally odd and no experimental constraints on them exist, but the 
lattice calculations of these GPDs [168] leading to the deformations shown in figure 27 
allowed predictions to be made for the signs of the Boer-Mulders functions for u and d 
quarks [126,172]. 

4-5. Access to quark GPDs through exclusive lepton-nucleon scattering 

Presently the experimentally most accessible GPDs are the chirally-even quark GPDs 
F*(x,£,t) (F = H, H, E, E and f — u,d). In order to constrain their non-forward 
(£ 7^ 0) behaviour, measurements can be performed of hard exclusive leptoproduction 
of a photon or meson, leaving the target nucleon intact. Useful observables are cross 

+ An operator that can be written as the same expression in all gauges is called manifestly gauge 
invariant. 



Spin-polarized high- energy scattering of charged leptons on nucleons 



65 



sections and cross section differences or asymmetries with respect to the polarization 
of the beam or target, or the beam charge. In these observables, GPDs typically 
appear convoluted with hard scattering amplitudes, so that they cannot be extracted 
directly from the measured data. Instead, their functional dependence on x, £, and t is 
parameterized in theoretical models (see section 4.6), the parameters of which have to 
be constrained by comparison to experiment (see section 4.7). For various final states, 
the convolution integrals appear combined additively or multiplicatively in various ways, 
with various kinematic prefactors. Hence eventually the extraction of quark GPDs will 
require a (presently not yet existing) global fit to data from various experiments, as has 
been customary for decades in the extraction of ordinary PDFs. 

4-5.1. Deeply virtual Compton scattering The exploitation of exclusive production of a 
real photon, i.e., Deeply Virtual Compton Scattering e p — > e p 7, has two advantages: 

i) among the hard exclusive processes that can be accessed experimentally in the 
foreseeable future, it is considered to be the one with the most reliable theoretical 
interpretation 

ii) effects of next-to-leading order [173-175] and sub-leading twist (see section 1) 
corrections [176-178] are under theoretical control. For suitably chosen observables 
the lowest correction is of order 1/Q 2 instead of 1/Q. 

In the generalized Bjorken limit of large photon virtuality Q 2 at fixed xb and t, 
the dominant pQCD subprocess of DVCS is described by the 'handbag' diagram shown 
in the left panel of figure 28. The skewness parameter £ is given by £ ~ 2-x B * n ^ s 




Figure 28. Leading-order processes for leptoproduction of real photons. Left: Deeply 
Virtual Compton Scattering. Middle and right: Bethe-Heitler process. 



limit. Pictorially, the parton (of flavor q) taken out of the proton carries the longitudinal 
momentum fraction x + £ and the one put back into the proton carries the fraction x — £. 
The quark GPD F f (x,£,t) can then be considered as describing the correlation between 
these two quarks at the given value(s) of t (and Q 2 ). 

The ^-dependence of quark GPDs is directly accessible in DVCS although high 
experimental precision, i.e., high statistical accuracy in conjunction with sufficient 
resolution of the experimental apparatus, is required to extrapolate to the limit t — > 0. 
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The latter is of particular importance for the evaluation of the second x-moment of the 
sum of the two 'unpolarized' quark GPDs, H* + E* , which is related to the total angular 
momentum Jf of the parton species / = (u,d,g), at a given value of Q 2 [19,123] (see 
(37)). 

The Bethe-Heitler (BH) process, or radiative elastic scattering, is illustrated in the 
two right panels of figure 28. Its final state is indistinguishable from that of the DVCS 
process and these two mechanisms have to be combined on the level of the process 
amplitudes. The differential cross section for leptoproduction of real photons is written 
in terms of the BH and DVCS process amplitudes as: 
do~(ep —> epj) 



dxBdQ 2 d\t\d(j) 



oc r 



BH I 



\ T DVCS\ ^ T DVCS T BH + T DVCS T BH • 



(61) 



I [ere <j) is the azimuthal angle between the scattering plane spanned by the trajectories 
of the incoming and outgoing leptons and the production plane spanned by the virtual 
photon and the produced real photon (see figure 29). The BH process amplitude is 
calculable using the electromagnetic proton form factors, which are well measured at 
small t (see, e.g., [179]). The pure DVCS contribution, \t dvcs \ 2 , can then be extracted 
by integrating over the azimuthal dependence of the cross section. In this situation the 
interference term I vanishes to leading order in 1/Q; its total contribution at collider 
kinematics was estimated to be at the percent level [178]. 




Figure 29. Definition of the azimuthal angle </> in DVCS in the target rest frame. 

The interference term I is of high interest, as the measurement of its azimuthal 
dependence opens experimental access to the complex-valued DVCS process amplitude. 
In this way, both their real and imaginary parts (magnitude and phase) become 
accessible [180]. This method to study DVCS can be compared to holography [181] 
in the sense that the phase of the process amplitude of DVCS is measured against 
the known 'reference phase' of the BH process. Moreover, at fixed-target energies, the 
more abundant BH process can be considered as an 'amplifier' to study the rare DVCS 
process. 
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An expansion in 1/Q of the electroproduction cross section (61) reveals distinctive 
kinematic signatures of DVCS [178, 180, 182]. For the interference term I this expansion 
reads: 

I oc ei ^ k n Re ^n cos n<p + P beam ^ K Im ^n sin n0 + 0(1 IQ) , (62) 

_n=l,3 n=l,2 

where e% = ±1 is the charge of the lepton beam and the k% s are kinematic factors. 
The M n represent linear combinations of helicity (Compton) amplitudes describing the 
DVCS process. The subscript n = 1,2,3 labels the helicities for initial and final state 
photons: n = 1 (++) means helicity conservation, n = 2 (+0) single-helicity flip and 
n = 3 (H — ) double- helicity flip. 

In general, the portion of the interference term containing real parts of Compton 
amplitudes can be isolated by forming a cross section difference or asymmetry with 
respect to the charge of the lepton beam [183], while that containing imaginary parts 
can be accessed by forming single-spin differences or asymmetries with respect to the 
polarization of the lepton beam [184] or of the target [131, 178]. In the interpretation 
of experimental results in the subsequent sections, always only the helicity-conserving 
twist-2 Compton amplitude Mi = M ++ is considered. In the convention of Ref. [178], 
this amplitude can be written as a linear combination of products of the Dirac and 
Pauli elastic nucleon form factors F\ and F 2 with Compton form factors (CFFs) T 
[T = H,£,H,£). For the three different target polarizations (U for unpolarized, L (T) 
for longitudinally (transversely) polarized target) the expressions are: 

t 



Mu cx 
M L oc 
M T cx 



Vh=t 

2m 



t 



2m 

tn-t 



F l K + i{F l + F 2 ) H 
£ (F 1 + F 2 )H + F 1 H 



4m 2 



F,E 



F 2 H-F 1 £ + £(F 1 + F 2 ) (tS) 1 



t 



Am 2 



(63) 
,(64) 
(65) 



4m 2 

where terms proportional to £ 2 and higher are neglected in (64) and (65) [131,178]. Here 
m is the nucleon mass and -t = 4£ 2 m 2 /(l - £ 2 ) is the minimum possible value of —t 
at a given £. 

The complex-valued CFFs J 7 ^, t, Q 2 ) are flavor sums of convolutions of the 
corresponding leading-twist quark GPDs with the functions CJ describing the hard 
7*g Compton scattering. The latter are available up to NLO in pQCD [173-175]: 



t,Q 2 ) = J2 f dx C J log g 2 ) (x, t , g 2 



(66) 



Here the — (+) sign in the superscript applies to the CFFs T = H,£ (TC,S), 
corresponding to the GPDs F* = H^E^ (H-^,E^). At 0(a s ), also contributions from 
F G (F G ) exist in (66) [131]. Except when explicitly needed, the dependence on Q 2 is 
omitted in the following. 
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(67) 



Only at leading order in a s there exist relatively simple relationships between the 
real or imaginary parts of a CFF and (the flavor sum over) the respective quark GPDs 
that it embodies [178]. Using the same sign convention, the real part of (66), 

Re Hi, t) L =° £ e\ \v £ dx F f(x, £, t) {-^ T 

becomes at given skewness £ a flavor sum of integrals over — 1 < x < 1 of the respective 
quark GPDs, where V denotes Cauchy's principal value. The imaginary part of (66), 

Im F(£, t) = tt ^ e) (F'{£, £, t) T F'(-£, £, *)) , (68) 
/ 

becomes a direct flavor sum of quark GPD singlet combinations Ff(£, £, t)^pFf(—£, £, t). 
Singlet combinations are also probed in exclusive vector meson production while 
in exclusive pseudoscalar meson production nonsinglet combinations Ff(£,£,t) ± 
F? t) are accessed [131]. 

A recent theoretical development exploits dispersion relations for CFFs. A CFF 
F satisfies a single-variable dispersion relation showing that the entire £ dependence 
of the real part is given by aterm that can be constructed from only the imaginary 
part [185-187]: 

ReF(U) = \vj d£'Im.F(£^) (-L^ T +Cjr(t). (69) 

At leading order, the imaginary part depends only on the GPDs on the trajectories 
x = ±£. The D-term form factor is then given by: 



where the D-term D* was introduced in (44) of section 4.2. This form factor vanishes 
for the CFFs H and £ and is identical up to the sign for TC and £ [186, 187]. Equation 
(69) together with (68) imply that at leading order and without considering effects of 
evolution, the real part of the CFF can be decomposed into a term that shares the 
well-known property of the imaginary part in depending on the value of the GPDs only 
along the one-dimensional trajectories x = ±£, and another term that depends only on 
t and Q 2 , and not on £. Thus, in a leading-order analysis of DVCS data, parametric 
models for quark GPDs can be constrained only in the subspace (x = ±£, £). Beyond 
this scenario, a full (x,£) mapping of quark GPDs is still possible, at least in principle: 

i) when a large enough range in Q 2 is covered by DVCS measurements, the known 
<5 2 -evolution of quark GPDs can be used to constrain their x-dependence. 

ii) in hard exclusive leptoproduction of a virtual photon (double DVCS or DDVCS), its 
virtuality, i.e. the invariant mass of the produced lepton pair, is an additional variable 
that facilitates a complete mapping of quark GPDs [188,189]. However, the DDVCS 
cross section is suppressed by an additional factor a em , thereby making this reaction 
practically inaccessible using present facilities. 
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Even if complete DVCS data sets were available for both proton and neutron targets, 
it would still be a very intricate theoretical problem to deconvolute the quark GPDs F* 
from the CFF T . The variable t does not pose any problem as it does not appear in 
the kernel; i.e., the deconvolution would have to be done at a given value of t. 

4-5.2. Hard exclusive meson production Additional information on GPDs is available 
from hard exclusive lepto-production of a pseudoscalar or vector meson M: e p — > e p M. 
To leading order in 1/Q and a s , the meson structure itself enters only through a few 
constants [131]. (In higher order the meson distribution amplitude has to be explicitly 
taken into account.) The meson production cross section receives contributions from 
both longitudinal (L) and transverse (T) virtual photons: 

1 1 spins spins 

where e is the ratio of longitudinal and transverse fluxes of the virtual photon. The 
dependences of the scattering amplitudes A on the hadron polarizations is summed or 
averaged over in ^2 spins . A factorization theorem has been proved to date only for the 
longitudinal part [190]. 

The contribution by transverse photons to the cross section for meson production 
appears at subleading order in 1/Q [131, 190]. In case of production of a vector meson 
its spin-density matrix can be measured, thus facilitating the longitudinal-transverse 
separation. In case of pseudoscalar meson production only data at sufficiently large Q 2 
may be used where the contribution of transverse virtual photons can be neglected. 

It appears that for meson production by longitudinal photons, at leading order in 
1/Q and for an unpolarized target [131,190], only one pair of quark GPDs is involved in 
the relevant combination, controlled by the type of spin-parity exchange in the meson 
production, so that by measuring various meson types different pairs of quark GPDs are 
filtered out. For natural-parity exchange, i.e. for the production of e.g. vector mesons 
such as p°, the quark-helicity-conserving GPDs H* and E* are involved: 

\ E \M 2 = (1 - ai^l 2 - (£ 2 + ^)|£| 2 - 2eRe(S*H). (72) 

spins 

For unnatural-parity exchange, i.e. for the production of e.g. pseudoscalar mesons such 
as pions, the quark-helicity-flip GPDs H* and E* are involved: 

\ £ \m 2 = (i - em 2 - e^\£\ 2 - 2eRe(£*n). (73) 

A'A 

Recently, detailed numerical studies were presented of the perturbative convergence 
of the longitudinal cross section for the production of p°, u, (j) and 7r + mesons [191]. At 
intermediate to large x B , typical for fixed-target experiments, NLO corrections of up to 
100% were found, which somewhat decrease in size when going from Q 2 = 4 GeV 2 to 
9 GeV 2 . At lower x B , typical for collider kinematics, NLO corrections are huge even 
for Q 2 well above 10 GeV 2 . Such an increased sensitivity to higher-order corrections is 
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due to the above shown quadratic dependence on GPDs. It also exists in DVCS (see 
(61)) for the \tbvcs\ 2 term, but not for the interference term I. We note that huge NLO 
corrections up to 200-300% appear also in other, mainly collider processes. Here in most 
cases the full NLO result is available and used, and typically the NNLO corrections are 
not again 100% of the NLO ones. 

4-6. GPD parameterizations 

The application of the Ji relation (37) to evaluate the total angular momentum of 
quarks of flavour / is a challenging task as it requires, at some fixed common value 
of £, knowledge of both H* and E* over the entire x range in the limit of small t. 
In section 4.5.1, it was explained that the best presently practical probe to constrain 
quark GPDs is DVCS, which unfortunately does not directly provide all of the required 
information because of its kinematic limitations. While H* (x, 0, 0) = q/{x) is well 
known, E'(x, £,0) is poorly constrained at £ = 0, because there it has no influence 
on experimental observables aside from its first moment, i.e. the Pauli form factor. 
Hence the Ji relation must be evaluated in the nonforward regime. Furthermore, 
the integral (37) must be evaluated in the limit t — > 0, even though the crucial 
dependences of DVCS cross sections on beam charge and polarisations of the beam 
and target disappear in this limit. Thus GPD parameterizations constrained by DVCS 
data must be used to extrapolate as t — * 0, which might seem to introduce some 
degree of arbitrariness. However, GPDs are far from arbitrary functions of the three 
variables x,£,t (or four, when including Q 2 ), as they must satisfy several relations as 
described in section 4.2. Considerable ingenuity is therefore required in finding analytic 
representations or phenomenological parameterizations of quark GPDs that embody the 
constraints of various symmetries and principles, while conveniently parameterizing the 
remaining degrees of freedom. The discussion of dynamical models, a recent review of 
which is given in [192], is beyond the scope of this paper. 

The simplest ansatz that one might consider for the t dependence of quark GPDs 
is a 'factorized' form: 



However, there is evidence from phenomenological considerations [125, 193, 194] and 
lattice QCD calculations [165,168] that the t dependence is entangled with both x and 
£, so that factorized ansatze are disfavoured. As discussed in section 4.3, t is conjugate to 
the impact parameter b±, so that this entanglement corresponds to the interplay between 
longitudinal and transverse degrees of freedom in the nucleon, sometimes called 'nucleon 
tomography'. Measured electromagnetic form factors together with some physically 
motivated functional forms for the limits x — > 0, 1 have provided some information 
about the entanglement of x and t for 'valence' quark flavours in the limited regime 
£ = [193]. In this work, an exponential ansatz was chosen for this entanglement in 
the valence GPDs: 



Hf(x,Z,t) = h f (x,Z)F{(t). 



(74) 



H f (-\x,0,t) =q fv (x)exp[tf f (x)] 



(75) 
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Figure 30. Impact parameter relationship between a struck quark with longitudinal 
momentum fraction x and the collection of 'spectator' partons. The dashed line 
indicates the center of momentum of the spectators and the thick continuous line 
the center of momentum of the proton. The figure is taken from [193]. 



where the unpolarized valence PDF qf v {x) = q/{x) — q~f(x) is known. Although the 
dependence on Q 2 has been suppressed here for convenience, this exponential form is 
approximately stable under QCD evolution, more so at small x. It corresponds to 
an IPDPDF with guaranteed positivity and a plausible Gaussian shape with an in- 
dependent width: 



f . , i qfv(x) 



<bl)= 4f f (x). (76) 



4//(x) 

The choice for the function //(x) is motivated by simple physical ideas. Partons with 
small x can be considered to arise from a cascade of branching processes, sometimes 
called 'Gribov diffusion' [195]. This sort of random walk corresponds to the mean-square 
impact parameter (b 2 ± ) growing at small x as log(l/x), a behavior that is provided by 
the 'Regge-motivated' form 

HH-\x,0,t) oc e^' 2 * nX\ fv {x) = e a ' tlo ^+^ t q fv (x) . (77) 

The choice for the function ff(x) is also influenced by plausible expectations for the 
opposite situation with x approaching unity. As was mentioned in section 4.3, a struck 
quark with large x is expected to have a small impact parameter b± with respect to 
the nucleon centre of momentum given by (55). With the center of momentum of the 
spectator collection given by — b±x/(l — x), the separation b±/(l — x) between quark 
and spectators can be interpreted as a lower limit on the transverse size of the nucleon, 
as illustrated in figure 30. If the size of this confined system is not to diverge as x — > 1, 
then (b 2 ± ) must vanish at least as rapidly as (1 — x) 2 [133, 196]. This implies that 
ff(x) — > Af(l — x) n as x — > 1, with n > 2. With n — 2, this as well as the above 
constraint for small x are respectively satisfied by the third and first terms in the form 

f f ( x ) = a ' [l - xf log - + B f {\ - xf + A f x{l - x) 2 , (78) 

where the second term provides a smooth interpolation between them, and a', Af and 
Bf are free parameters to be fitted. Using in (75) and (32) this form for H^~\x,0,t), 
or a similar form for E^~\x, 0, t) together with a simple ansatz for its unknown forward 
distribution constrained by the known magnetic moments of the proton and neutron, a 
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Figure 31. Square root of the mean-square impact parameters (6j_) for the sum 
and the difference of u and d quark distributions from the fit to elastic form factors, 
compared with lattice QCD results from [197]. The figure is taken from [193]. 

good fit is obtained of existing data for both Dirac and Pauli elastic form factors of the 
nucleon. For x < 0.8, H u (~\x,0,t) is well determined, thereby determining the mean 
square impact parameter b± of u quarks as a function of x in this range. As illustrated 
in figure 31, there are clear indications that the transverse distribution of d quarks is 
broader than that of u quarks at large x. While the fits for E^~'(x, 0, t) are ambiguous 
due to the unknown forward distributions, their second x-moments appearing in the Ji 
relation (37) are reasonably well determined within the context of the assumed ansatz. 
It is found that the contributions of the orbital angular momenta of u and d valence 
quarks to the nucleon spin almost cancel, leaving a small net contribution. In summary, 
a remarkable amount of information about 'nucleon tomography' was gleaned from only 
data for nucleon elastic form factors. 

In the following, two specific parametric models for quark GPDs that have been 
most frequently used (in several variants) for comparisons to DVCS data are briefly 
described as examples. Both are based on quark GPD parameterizations that might 
in principle be employed at any order in a s , while the phenomenology has been done 
at only leading order. The remainder of this section has limited goals as a reference 
tool for continuing students of the field, and might be omitted on a first reading. The 
first example, the 'VGG-model' [125,198], is based on an ansatz that factorizes the 
dependences on £ and t in the framework of double distributions including the earlier 
mentioned D-term (see section 4.2). Here the flavour dependence of the D-term is taken 
from the chiral quark soliton model [199], resulting in a somewhat modified expression 
for a generic unpolarised quark GPD F (F = H, E): 

Ff(x,£,t)= f d/3 f m da6(x-P-£a)Fl )D ((3,a,t)±-^9(t-\x\)D(^t),(79) 
J-i J-i+\P\ iV / £ 



Spin-polarized high- energy scattering of charged leptons on nucleons 



73 



where the upper (lower) sign applies for the GPD H (E). The model ansatz proposed 
in [200] is invoked to write the DD in (79) as a product 

F f DD {{3,a,t) = h^(P,a) F f { rw (/3, 0, t) (80) 

of a 'forward' (^-independent) quark GPD F/ orw (j3,0,t) with a 'profile function' 

ft(t >(ft")<x VW)^'' (81) 

that controls the skewness dependence through the (free) profile parameter b, with 
b — > oo approaching the forward limits (51) and (52) for H{j D (x,0,t) and H{y D (x,0,t), 
respectively. Often profile parameters are chosen separately for valence and sea quarks. 
We note that for b — > oo this ansatz for Fdd(x, £, t) is singular at x = for finite £, which 
is physically implausible. While the clear advantage of the ansatz (80,81) is its simplicity, 
it is considered to be too restrictive for flexible modelling of the (£, t) dependence of 
quark GPDs [201,202], especially at small values of £. The (x,t) dependence is based 
on either a factorized or the above-mentioned Regge-inspired ansatz. 

The quark GPD E* is evaluated from the 'pion pole', which is an ansatz based on 
the assumption that the dominant contribution at small \t\ arises from the emission of 
a virtual pion by the target nucleon [203,204], which for small |t| is not far off-shell 
because of the small mass of the pion. The 'pole' in the pion-exchange t dependence 
\j{rn\ — t) appears in the unphysical region at t — m^. This ansatz contributes to only 
the real part and only in the ERBL kinematic region. Also the pion quantum numbers 
limit this contribution to only the C-even part of the isovector combination E u and 
to the corresponding proton-neutron transition GPDs. For the unknown forward (spin- 
flip) PDF of the GPD E* , a parametric combination of valence and sea quark 'PDFs' is 
used with a functional form inspired by results obtained in [125] using the chiral quark 
soliton model [205]. In this parametric combination, the quark total angular momenta 
Jf of quarks and antiquarks of flavour / (/ = u, d) enter as model parameters, so that 
the derivation of mo del- dependent constraints on Jf becomes possible by comparing 
predictions to experimental data, as described in section 4.7.6. The VGG model, after its 
implementation in the 'VGG code' [198], was widely compared to polarised cross sections 
and asymmetries in beam spin and charge at fixed-target energies and found to describe 
essential features of the data, especially when using the Regge-inspired ansatz for the 
(x,t) dependence, (see section 4.7.2). As this model inherited the above mentioned 
limitations of the adopted DD ansatz, it cannot describe unpolarized cross sections 
at collider energies where the skewness is small. It was not studied how seriously its 
limitations may affect the various interpretations of experimental data in terms of this 
model. For a similar model [178], the effects of NLO corrections to the quark CFFs 
were found to be significant but moderate, so that they need to be included when 
extracting model parameters from precise data, while LO calculations were considered 
to be sufficient for the purpose of estimating observables. 

The so-called 'dual-parameterization' model, presented as the second example, is 
based on the idea of duality in soft hadron-hadron interactions, which is the assumption 
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that the 2^2 scattering amplitude in the s-channel can be represented as an infinite 
sum of t-channel exchanges [206]. In the case of quark GPDs, this idea is implemented 
in the form of a t-channel partial wave expansion of the Mellin moments [207]. In 
addition, by using Gegenbauer moments instead of Mellin moments, one can diagonalize 
the leading-order QCD evolution kernel, which is a major advantage of this approach. 
After resumming the partial wave expansion for the Gegenbauer moments, one arrives 
at an expansion for a quark GPD F* in which each subsequent term is suppressed by 
two powers of the skewness parameter £: 

Ff{x, e, t) = K Q f Q (x, t) + eK 2 Q f 2 (x, t) + ?K 4 Ql(x, t) + . . . (82) 

Here the Kk represent parameter-free ^-dependent (but t-independent) integral 
transformations acting on the variable x. The functions Q{(x, t) are called 'forward-like' 
functions as their LO evolution is governed by the DGLAP evolution equations, just as 
for PDFs. The parameters of the model appear in these functions. In the (t-dependent) 
'forward' limit £ = 0, (82) simplifies to only the first term and the K becomes simple 
enough to invert. In the case of the quark GPD , the functions Q Q (x,t) can thus be 
related to q f (x,t) = # /(+) (x, 0, t): 

- 1 dC 



Qi(x,t) 



(q f (x, t) + q f (x, t)) - | J ^(q f (C, t) + q f ((, t)), 



3) 



and Qq(x,0) is thus constrained in terms of the PDFs. An analogous expression holds 
for the GPDs E^ + \x,^,t) with (t-dependent) 'forward' distributions ef(x,t). The 
expansion in (82) is introduced in such a way that as k increases their terms are 
increasingly suppressed for small skewness £ and fixed x. It remains to be shown that 
this expansion in £ gives a good approximation for fixed x/£ when truncated at finite 
k [208]. 

The leading-order DVCS amplitude relates to quark GPDs through the convolution 
integral: 



LO 



dx F /(+) (x,£,t) 



x - £ + iO 



+ 



x 



£-z0 



(84) 



Another major advantage of the dual parameterization is that this singular integration 
can be performed analytically in terms of the Q[ [201]: 



1m A f (£,t)= -2 C —N f (x,t) 

Jx x 



f-x 2 -l 



55) 



dx 



RsA f (£,t)= -2 / —N f (x,t) 
Jo x 

° x dx 



l-J^f + x 2 l + J^f + X 



2£ i ™2 1 + X 2 



2 f ^Nf(x,t) 

Jx x 



1 



4Vf(t), (86) 
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where x = 1 ^ , N^(x,t) = XlfcLo xk Q{( x , t) and T>f(t) is the .D-term form factor. 
The integral transform in (85) can be inverted analytically, yielding 



-Invt'&t)-^— InU'«,t) 
2 d£ 



87) 



7T (1 + x 2 ) 3/2 7,o C 3/2 V 7 ^^ 

with xo = • This result can be summarized in the following statements that all hold 
to leading order only: 

• The imaginary part of the DVCS amplitude alone (but for all £) is sufficient to 
determine N*(x,t). If data are available only above a certain value of £, one can 
still reconstruct N^(x,t) above a certain value of x. In principle, i.e. given ideal 
data (negligible uncertainties), this reconstruction can be done exactly. 

• Upon inserting N-f(x,t) from (87) into (86), data for the real part of the Compton 
amplitude can be used to determine the D-term form factor. Given ideal data with 
full coverage in £, the only new information contained in ReA^(£,t) that was not 
already contained in Im v 4^(^, t) is the D-term form factor. Of course, as ideal data 
are unfortunately unavailable, (86) also provides useful constraints on N*(x,t) as 
the D-term form factor does not depend on £. 

• At fixed Q 2 , N-f(x,t) and the D-term form factor constitute the only information 
about quark GPDs that can be obtained from DVCS. However, QCD evolution can 
be a viable tool for the extraction of quark GPDs, as the LO evolution of the Q^ k is 
known to be governed by the LO DGLAP evolution equations. Since the Q2 k (x,t) 
enter N^(x,t) multiplied by x 2k , the n'th Mellin moment of N^(x,t) will depend 
on the (n + 2fc)'th Mellin moment of the Q{ k . Since the anomalous dimensions 
that govern the LO evolution depend on the order of the moment (and only on 
the order), this implies that at different Q 2 the Mellin moments of Q 2k (at some 
renormalization scale Qq) enter N*(x,t) with different relative weights. At least 
in principle, DVCS data taken at a set of different values of Q 2 thus provide a 
system of linear equations for the Mellin moments of the Q 2k that (after suitable 
truncation) can be solved, which in turn can be used to reconstruct the quark GPD. 

Unfortunately, not only is there a lack of ideal data that are required for the above 
procedure to work, but also LO evolution is likely to be insufficient at presently accessible 
values of Q 2 . An approach that is closely related to [201] is that of [209,210], which uses 
conformal symmetry to diagonalize the evolution up to NLO and may thus overcome 
the above discussed limitations of the LO treatment. 

The dual parameterization does not model the t dependence of the GPDs, which 
must be specified separately. Present applications of the model do not include the 
quark GPDs H and nor higher-twist contributions. The model is worked out only 
to leading-order accuracy, as only then the DVCS amplitude has a simple form in terms 
of the Q k {t), and also their evolution can be handled easily. The generalization from 
LO to NLO is far from obvious. On the other hand, an advantage of this approach 
is the clear separation of the (^-independent) 'forward' distributions qj(x,t) [ef(x,t)] 
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of the quark GPDs H* (x, £, t) [E? (x, £, t)} from genuine non-forward effects encoded 
in the distribution of strength among the functions QI,Q{,... in (82). Restricting 
oneself to the two lowest functions (k = 0, 2) yields the so-called 'minimal version' of 
the model [207]. 

Also assuming in this 'minimal version' Q^XjO) oc Qq(x,0), it was initially found 
that the Q 2 and W dependences of the low- a; HERA collider data were fairly well 
described [211], until a missing factor was noticed, leading to a cross section magnitude 
discrepancy of about a factor of four [212]. In [213], calculations in the minimal version 
using only the lowest function (called 'zero-step' approach) as well as calculations of 
the VGG model were recently compared to JLab data on cross sections and asymmetries 
(see section 4.7.3, e.g. figure 39). These comparisons are unaffected by that missing 
factor. 

No global fits to all existing DVCS data yet exist. Future algorithms must 
necessarily include enough orders beyond LO to ensure convergence. The only existing 
approach fulfilling this requirement is a recently developed formalism [186] based on 
dispersion relations discussed in section 4.5.1. In this new formalism, dispersion relation 
and operator expansion techniques are combined leading to convergence of the conformal 
partial wave expansion of the VCS amplitude. This formalism was successfully applied 
to fit DVCS observables in NNLO at low xb [186], as described in section 4.7.1. In a 
first attempt, this fit was very recently extended to include also data from fixed-target 
experiments [214]. 

4-7. DVCS experimental results 

4-7.1. DVCS cross section at high energies The DVCS cross section (61) integrated 
over its azimuthal dependence has been measured in hard exclusive electroproduction 
of photons by the experiments HI and Zeus at the Hera collider. For the x B range 
accessible at collider kinematics (O(10~ 2 . . . 10 -4 )), two-gluon exchange plays a major 
role besides the quark-exchange handbag diagram of figure 28, so that both gluon and 
quark GPDs are probed simultaneously. This is however limited to small skewness 
values below 1CT 2 . 

The analysis method used by the two experiments is similar due to their similar 
geometries. As the outgoing proton remains undetected in the beam pipe, the event 
topology is defined by two electromagnetic clusters, the outgoing lepton and the 
produced real photon, at most one of which is associated with a charged track. Two 
events samples are selected: 

i) in the 'DVCS-enriched' sample an energetic photon is required. With respect to 
the incoming lepton it is required to be detected at a large scattering angle while the 
outgoing lepton is measured at a small angle. Still, a high enough virtuality Q 2 > 4 GeV 2 
is ensured by requiring a large energy of the scattered lepton (> 15 GeV). 

ii) in the Bethe-Heitler dominated 'reference sample', the radiatively produced photon 
is emitted at a small angle with respect to the incoming lepton, and the outgoing lepton 
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at a large angle. 

A Monte Carlo simulation of the completely known BH process, validated by the 
reference sample, is used to subtract the BH contribution from the DVCS-enriched 
sample. The remainder of the spectrum is due to DVCS and possible additional 
background; no contribution from the interference term exists at leading twist, as the 
data are integrated over the azimuthal angle. 
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Figure 32. Left: Q 2 -dependence of the differential -f* p — > 7p cross section measured 
by HI and Zeus. Right: as left panel, but here the W^-dependence is shown. The 
figures are taken from [215]. 



The <3 2 -dependence of the differential 7*p — > 7p cross section is shown in the 
left panel of figure 32, extracted from Hera-I and Hera-II data of HI [215] and the 
Hera-I data of Zeus [216], the latter based on a somewhat larger event sample. The 
corresponding W dependence, where W is the invariant mass of the system of virtual 
photon and proton, is displayed in the right panel. The virtuality appears high enough 
to assign the observed rise with W to the nature of DVCS as a hard process, as increasing 
W implies decreasing x B , where the parton densities in the proton show a rapid rise. 
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Figure 33. Hera-I data shown in the previous figure in comparison to various Colour 
Dipole Model calculations (see text). The figures are taken from [217]. 



In figure 33, various Colour Dipole Model calculations [218-220] are compared to 
only the Hera-I data of both experiments. Presently it is impossible to discrimate 
between the various model versions. 
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Figure 34. The —t dependence of the differential 7*p — > 7^1 cross section, extracted 
from HI Hera II data at 3 values of (Q 2 ). The figure is taken from [215]. 



A commonly used ansatz for the t dependence of quark GPDs is an exponential 
decrease e~ fe '*L The slope parameter b was first measured in DVCS by HI [215], the 
results being shown in the two left panels of figure 34 for (Q 2 ) = 8, 15.5 and 25 
GeV 2 integrated over W, and for (W) = 40, 70 and 100 GeV integrated over Q 2 . 
Parameterizing a possible Q 2 dependence as b{Q 2 ) — A (1 — B • log(Q 2 /2 GeV 2 )), the 
fit shows that it is indeed significant with A = 6.98 ± 0.54 Gev 2 and B — 0.12 ± 0.03. 
At Q 2 = 8 GeV 2 the result is b = 5.45 ± 0.19 stat ± 0.34 syst . Following [135,221], this t 
slope can be converted to an average impact parameter which describes the transverse 
extension of partons. The HI result at (x) = 1.2 ■ 10~ 3 , where sea quarks and gluons 
dominate, is \J (r^) = 0.65 ± 0.02 fm. This can be compared to the transverse charge 
radius of the proton [222]: 

VVtU = J-4 ■^F 1 (Q 2 )\ Q 2 =0 = 0.65 ± 0.01 fm. (88) 

Recently, the convergence of the perturbative treatment of DVCS was studied in 
the low-£ region typical for collider kinematic conditions, in a theoretical study using 
the predictive power of conformal symmetry [223]. Beyond NLO, further terms in the 
perturbation series are found to be small for an input scale of a few GeV 2 . Concerning 
the scale dependence, larger corrections due to evolution may appear only for values of 
£ below 5 • 10 -3 , similar to the situation in DIS [224]. 

Based on this approach, very recently a simultaneous fit was performed to the 
DVCS and DIS collider data, for the first time up to NNLO accuracy [186]. The CS 
scheme was used which differs from the more commonly used MS scheme only in the 
skewness dependence of the conformal moments in the Mellin-Barnes representation, the 
latter being used to evaluate the Compton form factors. In figure 35 the fit results are 




Figure 35. Simultaneous fit to the DVCS and DIS data in the CS scheme at NNLO. 
Upper left panel: DVCS cross section for Q 2 = 4 GeV 2 and W = 71 GeV (circles, 
dashed line) as well as Q 2 = 8 GeV 2 and W = 82 GeV (squares, solid line) [217]. 
Upper right panel: DVCS cross section (— t < —1 GeV 2 ) versus Q 2 for W = 82 
GeV (HI, circles, dashed line) and W = 89 GeV (ZEUS, triangles, dash-dotted 
line) [216]. Lower left panel: DVCS cross section versus W for Q 2 = 4 GeV 2 (HI, 
circles, dashed line), Q 2 = 8 GeV 2 (HI, squares, solid line), and Q 2 = 9.6 GeV 2 
(ZEUS, triangles, dash-dotted line). Lower right panel: i*2(xB,Q 2 ) versus Q 2 for 
x B =8 - 10 -3 ,3.2 • 10~ 3 ,1.3 • 10" 3 ,5 • 1CT 4 [225]. The figure is taken from [186]. 



compared to the data, showing fair agreement. It is argued in [186] that this should be 
considered satisfactory, as several aspects of the theoretical approach are not completely 
settled yet, in particular regarding the skewness dependence of quarks GPDs. 

4-7.2. DVCS cross sections at low energies The helicity-dependent DVCS cross 
sections d A a~~* /dQ 2 dx B dtd<j) and d 4 o~*~~ /dQ 2 dx B dtd<f) were measured at JLab with 
5.75 GeV electron beams of positive and negative helicity, respectively [226]. Here 
the symbol ---» («---) designates positive (negative) beam helicity. In a three-arm 
kinematically complete experiment, the scattered electron was measured by the high- 
resolution spectrometer, the produced real photon by an electromagnetic calorimeter 
and the recoil proton by a scintillator array. Three photon virtualities were selected 
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with x B = 0.36 fixed: Q 2 = 1.5, 1.9, 2.3 GeV 2 . For each of them, four —t values were 
chosen: 0.17, 0.23, 0.28, 0.33 GeV. 

From the two measured single-helicity cross sections, the polarized and unpolarized 
DVCS cross sections were obtained: 

d 4 £ 1 / dV- ciV" 

2 \dQ 2 dx B dtd(fi dQ 2 dx B dtd<p / 

1 / <Pa~* d A a*- 

2 \dQ 2 dx B dtd<p dQ 2 dx B dtd<p , 



dQ 2 dx B dtd(p 
d 4 a 



oc Im I 



oc Re I, 



(89) 



(90) 



dQ 2 dx B dtd(f) 

where I is the interference term that consists of linear combinations of Compton Form 
Factors. Cross section measurements involve additional experimental systematics, as no 
cancellations of systematic uncertainties occur as in asymmetry measurements. Such 
measurements have the advantage that the cross sections in (89) and (90) are directly 
proportional to Im I and Re I, respectively, while in case of asymmetries unwanted 
Bethe-Heitler terms appear in the denominator. 
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Figure 36. JLab data and fit to e? 4 E / 'dQ 2 dx B dtd<f) and d 4 a/dQ 2 dx B dtd(f>, shown as 
a function of the azimuthal angle <j> at Q 2 = 2.3, —t = —0.28 at (x B ) = 0.36. Only 
statistical errors are shown. The solid line shows the total fit with the band representing 
one standard deviation. The dot-dot-dashed line is the |BH| 2 squared contribution, 
twist-2 (twist-3) interaction terms are shown by a long-dashed (dot-dashed) line. The 
figure is taken from [226]. 

The data for the polarized and unpolarized DVCS cross sections are shown in 
figure 36 together with Fourier decompositions according to (62). The polarized cross 
section exhibits the expected sin0 behaviour. For the unpolarized one, the main 
contributions stem from the |BH| 2 term and the twist-2 part of the interference term, 
while | DVCS | 2 and twist-3 contributions are found to be negligible. 
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4-7.3. Beam-spin asymmetry It was shown in section 4.5 that in hard exclusive 
leptoproduction of a real photon, the interference of the Bethe-Heitler and Deeply 
Virtual Compton Scattering processes allows the extraction of a wealth of information on 
quark GPDs. For the case of an unpolarized (U) proton target, the beam-spin asymmetry 
(BSA) for a longitudinally (L) polarized beam is given by: 

da"\(j)) -da*~(<f>) 



A LU ((j>) 



(91) 



da- + 0) + da^-{<p)' 

Here ---> (<---) denotes beam spin parallel (antiparallel) to the beam direction. 

The interference term I can be Fourier expanded in the azimuthal angle <j). 
Evaluating the above asymmetries to leading power in 1/Q in each contribution to 
this expansion and to leading order in as, only the sin0 (cos0) term remains in 
the numerator of the beam-spin (beam-charge) asymmetry. To the extent that the 
leading BH-term dominates the denominator, the azimuthal dependence of the beam- 
spin (beam-charge) asymmetry is reduced to sin <f) (cos 0) : 

A LU ((j)) oc Im M v sin 0. (92) 

The BSA is sensitive to the linear combination Mjj of Compton Form Factors, which 
describes an unpolarized proton target [178] (see (63)). Generally, the quark GPDs H* 
are expected to dominate expression (63), because i) the second term is suppressed by 
at least a factor of 10, as £ is usually not larger than 0.2 even in fixed-target kinematics, 
and the unpolarized contribution 7i is expected to dominate the polarized one 7i, in 
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Figure 37. Azimuthal dependence of the beam-spin asymmetry. Left: Hermes 
proton data taken with 27.6 GeV positrons. Right: Clas proton data taken with 4.25 
GeV electrons. The figures are taken from [227,228]. 



analogy to the forward case; ii) the third term is t-suppressed, by about a factor of 25 
for typical t- values of about 0.15 GeV 2 . For scattering on the proton, the quark GPD 
H u will yield the major contribution to Mu because u quarks i) are more abundant in 
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the proton and ii) more readily absorb virtual photons because of their 4 times larger 
charge-squared factor. 

The first published GPD-related experimental results were beam-spin asymmetries 
measured in DVCS on the proton by the fixed-target experiments Hermes at 
Hera [227] with a positron beam and by Clas at Jefferson Laboratory [228] with 
an electron beam (see figure 37). Note that opposite beam charges imply opposite signs 
of the measured beam spin asymmetries (this is not immediately apparent in figure 37 
because of different 0-ranges shown for Hermes and Clas). In both experiments 
significant sinusoidal modulations have been observed, as expected from the theoretical 
considerations that were briefly explained above. 

DVCS on the proton was recently measured at JLab with a 5.77 GeV polarised 
electron beam in a kinematically complete experiment using the large-acceptance Clas 
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Figure 38. Beam-spin asymmetry measured at JLab as a function of t and shown 
at <j> = 90° in several (x B , Q 2 ) bins. Systematic uncertainties and bin limits are shown 
by the (grey) band in the lowest left panel. (Black) circles are from [229], the (red) 
square is earlier data [228] and the open (green) triangle is the cross section data [226] 
described in the previous section. The curves are described in the text. The figure is 
taken from [229]. 

spectrometer. Beam-spin asymmetry results were shown in a simultaneous binning in 
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x, t and Q 2 (62 bins in total) (see figure 38 [229]). The background from undetected 
exclusive 7r°'s decaying asymmetrically was estimated using measured symmetric decays 
and their relative acceptance function determined by Monte Carlo techniques. It was 
subtracted in every kinematic bin, over which it varies from 1 to 25% (5% on average). 
Calculations using the earlier described VGG model including only the contribution of 
the quark GPD H, shown as two neighbouring lines, solid for twist-2 and dashed for 
twist-3, clearly overshoot most of the data. Calculations in a Regge model [230], shown 
for comparison in a few panels as the lower dashed curve, describe the data well in 
certain bins but not in others. 

For a few particular bins, the azimuthal dependence of the same data was compared 
in [213] to the VGG and the dual-parameterization models that were described in 
section 4.6. For both models, the same unfactorised ansatz (77) was used. For simplicity, 




Figure 39. Azimuthal dependence of the beam-spin asymmetry in DVCS in several 
kinematic bins. The black bands in the two middle panels are JLab data from [226], 
while the data points are JLab data from [229]. The curves are explained in the text. 
The figure is taken from [213]. 
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only the BSA is shown in figure 39, although in [213] also the polarised and unpolarised 
cross sections are discussed. The VGG calculation, for which b = 1 was used as profile 
parameter, overshoots the data as already discussed in connection with the previous 
figure. The dual-parameterization calculation, done with only the forward function 
Qo, also overshoots the data, although less severely. Adding some estimate for the 
non-forward function Q 2 does not improve the situation. Using the values 0, ±| for the 
D-term form factor V(t)\ t=0 that was mentioned after (87) does not yield a simultaneous 
consistent description of the <p dependence of the unpolarised and polarised cross sections 
and the BSA [213]. Such a description can not be achieved by any of the presently 
existing GPD models. A similar conclusion is reached in [231] after attempting a 
combined fit of the two JLab data sets [226, 229] to variants of the VGG model. It 
has not been tried to estimate the possible influence of higher-order effects for either of 
the two models. 



4-7.4- Beam-charge asymmetry The measurement of a beam-charge asymmetry (BCA) 

where the superscripts + and — denote the lepton beam charge, requires data for both 
beam charges. Alternatively, a BCA can also be derived from the cos0 dependence of 
the cross section for one beam charge in an hermetic detector of known efficiency. Using 
the same approximations as in the previous section, the azimuthal dependence of the 
beam-charge asymmetry is reduced to cos</>: 

Ac{4>) oc Re Mu cos0. (94) 

The BCA is sensitive to the real part of the same linear combination Mu (63), of which 
the imaginary part determines the BSA described in the previous section. 

Hera was the only multi-GeV accelerator providing both electron and positron 
beams. It offered the additional flexibility of inverting every few weeks, for the same 
charge of the beam, the direction of its polarization to reduce systematic effects. The 
first measurement of a beam-charge asymmetry was published by Hermes [232] , based 
on the analysis of parts of their electron and positron beam data sets. Somewhat later, 
BCA results were published based on the analysis of another much larger Hermes 
data sample, taken with transverse target polarisation [233]. Here, beam-charge and 
transverse-target-spin asymmetries were extracted simultaneously (see section 4.7.6). 
In figure 40, the experimental data on the azimuthal asymmetry amplitude A™ 5 ^ are 
compared to calculations using the model of [125,198] that is briefly described in section 
4.6. The model results are calculated at Hermes kinematics, each at the average 
kinematics of each individual bin. The best description within the VGG model is given 
by the version using an unfactorised Regge-inspired t dependence, no D-term and profile 
parameters b va i = oo (no skewing) and b sea = 1 (some skewing). 
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Figure 40. Beam-charge asymmetry measured on the proton by Hermes: azimuthal 
amplitude describing the dependence of the interference term on the beam charge 
(Ac). The triangles (shifted right for visibility) represent previous results [232], while 
most recent data [233] are represented by squares. The error bars (bands) represent 
the statistical (systematic) uncertainties. The curves labelled 'DD' are calculations of 
variants of a double-distribution GPD model [125, 198] using b v = oo and b s = 1 as 
profile parameters for valence and sea quarks. This figure is extracted from figure 5 
of [233]. 



4-7.5. Longitudinal target-spin asymmetry Experimental access to quark GPDs other 
than H' becomes possible when using an unpolarized beam (U) and a polarized target. 
For longitudinal (L) target polarization, the same sin0 dependence is expected as for 
the BSA, although the relevant twist-2 Compton amplitude Ml (64) involves a different 
combination of Compton Form Factors than does Mu (63). 
The longitudinal target-spin asymmetry is defined as: 
a (a\ do-^((p) - do-^{(j>) 

AuM = da^mTda^ify (95) 
where <= (=>) denotes target spin antiparallel (parallel) to the beam direction. This 
asymmetry is expected to be most sensitive to a combination of H' and Iff , because 
the kinematic suppression of the first term in (64), relative to the second, may 
compensate the usual dominance of the unpolarized quark GPDs fff over their polarized 
counterparts fff . Hence they might be disentangled by combining this measurement 
with BSAs. For not too small values of t, there exists also some sensitivity to (££), 
which is written in this way because £ itself is defined to be inversely proportional to £ 
(see e.g. [187]). 

Experimental results on longitudinal target-spin asymmetries were obtained at 
JLab [234] and Hermes [235]. Both experiments observe a sin0 modulation of the 
yield from the entire acceptance, and scaling of the sinusoidal amplitude with \f—i, 
in agreement with expectations. In figure 41, the JLab data [234] are shown. The 
continous curve represents a fit with the function a sin <p + (3 sin 2<p with the parameters 
a = 0.252 ± 0.042(stat) ± 0.020(syst) and (3 = -0.022 ± 0.045(stat) ± 0.021(syst). 
The sin 20 term is consistent with zero, indicating that higher-twist contributions in 
this kinematic domain are negligible. The two other curves show calculations using the 
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Figure 41. Azimuthal dependence of the longitudinal target-spin asymmetry, as 
measured at JLab [234] with 5.7 GeV electrons and a longitudinally polarized 
NH3 target, shown after subtraction of the ir° background. Error bars represent 
the statistical uncertainties and the band at the bottom represents the systematic 
uncertainties. The curves are explained in the text. The figure is taken from [234]. 

double-distribution model of [198], based on MRST02 PDFs, with the £ dependence 
Kai = b sea = 1, including leading- twist terms only, and target-mass corrections applied. 
The nucleon-spin-flip quark GPDs E and E are set to zero. The dashed (dotted) line 
shows the asymmetry when the polarized quark GPD H is included (excluded), so that 
the necessity of H for the description of the longitudinally target-spin asymmetry is 
clearly demonstrated. 

4-7.6. Transverse target-spin asymmetry In the case of target polarization transverse 
(T) to the direction of the virtual photon, the polarization can be resolved into two 
independent components, 'normal' to and 'sideways' in the production plane [236], 
defined by the directions of the virtual and real photons. The DVCS cross section 
depends only on Mt, the normal component of the twist-2 Compton amplitude (see 
(65)). This is the only amplitude that embodies a combination of CFFs in which the 
contributions of the quark GPDs E* are not kinematically suppressed relative to those 
of H* . Hence DVCS measurements on a transversely polarized proton target appear to 
be crucial for the evaluation of J u and Jd through the Ji relation (37). A complication 
inherent to the exploitation of this relation lies in the fact that both quark GPDs 
and E$ need to be determined in the limit t — > 0, whereas at small t the relevant 
asymmetry is suppressed by a factor of y/t — t/2m when extracting the GPDs and 
even by a factor of (t — t)/ 4m 2 when extracting the GPDs E* (which is apparent from 
a comparison of (63) and (65)). The extraction of mo del- dependent constraints on J u 
and Jd from experimental data is discussed in section 4.8. 

From data taken with both beam charges and transverse target polarization iSj_, 
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the beam-charge asymmetry Ac{4>) (see (93)) can be extracted simultaneously with the 
transverse target-spin asymmetries arising from the |T£>ycs| 2 an d interference terms (see 
(61)), respectively: 

1 da + ((p, 4> s ) - da + (0, 4> s + 7r) + da" (0, <p s ) - da~ (0, <p s + 7r) 



yjWCS 



l C7T 







[/TV 



5, 



6', 



+ Ax+( 
-d<r+( 



TTj 
7Tl 



<icx~ 



+ do- 
+ da~ 



>5 + 7T) 
5»S + tt) 



(96) 
.(97) 



da+ (0, S ) + c/cr+ (0, 0s + 7r) + da~ (0, 0s) + c/a - (0, 0s + 7r) 
>5 of the target polarization vector was explained in section 3.4. 



The azimuthal angle 

The subscripts on the A's represent the dependence on beam Charge (C) or Transverse 
(T) target polarization, with an Unpolarized (U) beam, and the superscripts ± stand 
for the lepton beam charge. All three asymmetries are expanded in terms of various 
harmonics in and 05 which respective amplitudes are hereafter called '(effective) 
azimuthal (asymmetry) amplitudes'. 
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Figure 42. Asymmetry amplitudes describing the dependence of the squared DVCS 
amplitude (circles. Ajjt,dvcs) and the interference term (squares, Ajjtj) on the 
transverse target polarization. The circles (squares) are shifted right (left) for visibility. 
The error bars represent the statistical uncertainties, while the top (bottom) bands 
denote the systematic uncertainties for Ajjtj {Ajjt.dvcs), excluding an 8.1 % scale 
uncertainty from the target polarization measurement. The curves in the upper panel 
are calculations of the GPD model variant (Reg, no D) with three different values for 
the u-quark total angular momentum J u and fixed <i-quark total angular momentum 
Jd = [237] . This figure is based on figure 6 of [233] . 

Only a few azimuthal amplitudes are expected to show substantial sensitivity to 
J u and Jd [237]. At Hermes all harmonics of physical significance were measured 
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with unpolarized beam and a transversely polarized proton target [233] . The kinematic 
dependences of two asymmetry amplitudes of particular interest are shown in figure 42 
in comparison with calculations based on two different types of GPD models. The data 
is clearly able to discriminate among these GPD models. The model variants that agree 
best with the data, are used in the next section to derive constraints from the data on 
the total angular momentum of u quarks. Sensitivity exists, as can be seen from the 
top panel of the figure. 

4-8. Quark total angular momenta 

In the theoretical description of the DVCS process, experimental observables like 
cross section differences or asymmetries are related to various combinations of the 
Compton Form Factors 7i, S, H., S, which can be calculated from the respective GPDs 
,H-t ,E* ,E-f (see section 4.5). Parameterizations of the spin-flip GPD E* embody 
explicitly or implicitly the quark total angular momenta J u and Jd [125]. Hence model- 
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Figure 43. Model-dependent constraints on u-quark total angular momentum J u vs d- 
quark total angular momentum J^, obtained by comparing DVCS experimental results 
and theoretical calculations using the double-distribution GPD model of [125,198]. 
The constraints based on the Hermes data [233] for the azimuthal amplitudes 
A si„W- fe )cos0 and A *m£*-M are labelled Hermes DD (VGG). Those based on 

the JLab neutron cross section data [238] are labelled JLab DD (VGG). Also 
shown as small (overlapping) rectangles are results from lattice gauge theory by the 
QCDSF [239] and LHPC [165] collaborations, as well as a result for only the valence- 
quark contribution (DFJK) based on zero-skewness GPDs extracted from nuclear 
form factor data [193,240]. The sizes of the small rectangles represent the statistical 
uncertainties of the lattice results and the parameter range for which a good DFJK fit 
to the nucleon form factor data was achieved. Theoretical uncertainties are unavailable. 
The figure is taken from [233]. 

dependent constraints on these two parameters can be derived by fitting them to 
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experimental observables [237]. The double-distribution model [125,198] was thereby 
fitted to azimuthal amplitudes of beam-charge and transverse target-spin asymmetries 
measured on the proton at Hermes (see previous section). The resulting constraint [233] 
is shown in figure 43 as sloped band, which in units of h~ can be represented as 

J u + J d /2.8 = 0.49 ± 0.17(exp tot ). (98) 

Combining DVCS measurements on deuteron and proton targets, the beam-spin 
difference for the neutron was determined at JLab [238]. Also in this case sensitivity to 
the GPDs can be found, as for the neutron the prefactor Fi(t) in (63) is small but 
F 2 (t) is not. Similarly as described above, a mo del- dependent constraint on (J u , Jd) was 
derived by comparison to a calculation of the double-distribution model. The slope of 
the resulting band differs from that of the band derived from proton data (see figure 43) 
due to the different quark contents of proton and neutron: 

J u + J d /5.0 = 0.18 ± 0.14(exp tot ). (99) 

These results show that existing data have the potential to provide quantitative 
information about the quark total angular momenta J u and J d . However, as discussed 
in the previous sections, the double-distribution model of [125, 198] cannot explain all 
existing DVCS data. 



5. The spin budget of the nucleon 

5.1. Decomposing the Nucleon Spin 

While the total angular momentum of an isolated system is uniquely defined, ambiguities 
arise when decomposing the total angular momentum of an interacting multi-constituent 
system into contributions from various constituents. Moreover, in a gauge theory, 
switching the gauge may result in shuffling angular momentum between matter and 
gauge degrees of freedom. In the context of nucleon structure, this gives rise to subtleties 
in defining these quantities that are more fundamental than those subtleties associated 
with the choice of factorization scheme, mentioned in section 2.4. 

In section 4.1, it was discussed how GPDs can be used to determine the quark total 
angular momentum contribution Jf to the nucleon spin (see (37)) in a gauge invariant 
decomposition of the nucleon's spin |. The expectation value of the z component of 
the quark orbital angular momentum operator can then be evaluated indirectly as the 
difference, 

1 



J f - - (Aq f + Aq f ) 



J q - -AS, (100) 



2 



which is otherwise presently inaccessible experimentally. Similarily one can indirectly 
evaluate the expectation value of the operator representing the gluon total angular 
momentum as the difference between the nucleon's spin 1 and J q : 

Jq=\- J q - (101) 
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One might consider going even further and use experimental results for AG to determine 
the gluon orbital angular momentum as the difference between Jq and AG, providing 
what appears to be a complete understanding for the nucleon spin budget. However, 
the difference between these two operators has no fundamental connection to orbital 
angular momentum other than being the difference between Jq and AG. 
Equation (100) is a special case of Ji's decomposition [19] 

■j= x; (£/+£/)+ (102) 
/ 

of the expectation value of the angular momentum operator. Since this decomposition 
of the angular momentum density M fluX = T fJ-u x x — T^ x x v , where T^ u = T Ufl is 
the energy-momentum tensor, in terms of manifestly gauge invariant local operators 
(see footnote on page 64) is based on the decomposition of the M 0y '-component 
G {x, y,z}), it is independent of the quantization axis. Here Sf is defined through 
the expectation value of the axial vector current t^ljjs'i^f, Lf through the expectation 

value of ip^rx (j> — gjtj i/jf, and Jq through the expectation value of fx Tr yE x B^j, 

where E and B are the QCD color electric and magnetic fields. A further decomposition 
of Jq into spin and orbital components using manifestly gauge invariant local operators 
has not been found and may be impossible. In general, spin (or 'intrinsic' angular 
momentum) is identified with terms in the operator that do not depend explicitly on r, 
in contradistinction to the orbital (or 'extrinsic') piece [241]. 

Being based on manifestly gauge invariant local operators, all three terms in (102) 
can be calculated in lattice gauge theory. As was discussed in previous sections, the 
expectation values Sj and Jf can be accessed experimentally, and Jq can, at least 
in principle, be accessed in exclusive deeply virtual heavy meson production. Beyond 
leading order in a s , all three terms in (102) depend on Q 2 , but the Q 2 dependence 
cancels in their sum. 

The main disadvantage of decomposition (102) is that for no choice of the 
quantization axis is there a partonic interpretation as a difference between number 
densities for the quark orbital angular momentum Lf or the gluon total angular 
momentum Jq. Only the quark intrinsic angular momentum Sf has such an 
interpretation (for the z quantization axis). This is because the operator definitions for 
Lf and Jq contain interactions (e.g. Lf contains the gluon field through the gauge- 
covariant derivative D = d — igA). With quarks and gluons thereby inextricably 
intertwined, it is impossible to unambigiously identify these terms with quark orbital 
angular momentum, and gluon angular momentum respectively [242]. 

An alternative 'light-cone decomposition' was proposed by Jaffe and Manohar [243] 
(see also [244,245]). While the Ji-decomposition is based on the M 0y component of the 
angular momentum density, the Jaffe-Manohar decomposition is based on the M +xy 
component and applies only to the z component of J. It can be written as: 

l - = J z = ^AS + AG + C (103) 
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The quantity AS appears also in (100), the special case of (102) in the helicity basis. The 
quark-intrinsic piece is the same in both components: the operator that measures this 
piece in M +12 is e +12 ~ip Jf, j-'-f5ip and the operator that measures it in M 012 is e 0123 ip^'~f3'j5ip, 
i.e. both operators are components of the axial vector current and therefore their 
expectation values are proportional to the axial charge, so that AS appears in both 
decompositions. This is not true for any other term in the decompositions! One can 
understand this point by analogy with the energy momentum tensor T^ v [242]: The 
energy of a single particle state can be expressed in terms of the matrix element of 
T 00 and decomposed into a quark and gluon piece. The 'momentum' in the infinite- 
momentum frame corresponds to the matrix element of T ++ and can also be decomposed 
into quark and gluon pieces. The quark and gluon pieces of these two objects are 
unrelated. 

In decomposition (103), each term has a partonic interpretation. The gluon spin 
contribution AG appears explicitly. It is experimentally accessible (see section 2.8) 
and can be defined as the expectation value of a (nonlocal) manifestly gauge invariant 
operator. In light-cone gauge, this operator collapses to a local operator (and its 
expectation value has a partonic interpretation). No direct experimental access to the 
parton orbital angular momentum £ has been identified. Its value can be obtained only 
by subtracting the quark and gluon spin contributions from the nucleon spin. Both 
AG and £ can be defined through matrix elements of local operators only in light-cone 
gauge A + = 0. Explicit definitions for the operators appearing in both decompositions 
can be found in [242]. Since neither one can be represented as the matrix element of 
a manifestly gauge invariant local operator, they cannot be analytically continued to 
Euclidean space and are thus inaccessible for lattice QCD. 

One can decompose £ further into contributions from gluons and each quark flavor 



where for example, the contribution £f from quarks with flavour / can be obtained from 
the expectation value of (xp y — yp x ) tpf + (with vpf + being the dynamical or 'good' 
component of the quark field operator ipf in that framework) for a nucleon polarized 
along the z direction. Each of the terms in (104) has a partonic interpretation. Even a 
decomposition with respect to the momentum fraction x is possible, but for none of these 
terms has experimental access been identified nor are they accessible in lattice QCD, 
since their manifestly gauge invariant forms are nonlocal in Minkowski space [246]. 
Nevertheless, their expectation values can be easily related to the (light-cone) wave 
functions of hadrons, which play an important role in phenomenology [247] . 

Summarising the above dichotomy, the two decompositions described above have 
only AS in common, which is also manifestly gauge invariant and experimentally 
accessible. The experimentally accessible quantities Jf (and, by subtraction, Lf and 
Jg) appear in one decomposition, while AG (and, by subtraction, £) appear in the 
other. Combining terms from these incommensurate decompositions would be a sterile 
exercise. For example, after determining Jg = \ ~ J2f Jf by evaluating the right hand 




(104) 
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side using the Ji relation for Jf, it would be fruitless to 'determine' L G by subtracting 
AG from J G (or equivalently by subtracting L q from £). In fact, the last step would 
yield a term that is neither directly experimentally accessible nor could it be represented 
by an operator other than that defined as the difference of the two operators. Defining 
gluon orbital angular momentum through such a subtraction procedure has no predictive 
power (not even in principle). 

Note that an analogous dichotomy (multichotomy) appears when one considers the 
motion of a charged particle in an electromagnetic field. While the canonical momentum 
p = — has a simple interpretation in terms of only the charged particle's wave 
function, it is not gauge invariant. Evaluation of the gauge invariant combination p — eA 
requires in addition knowledge about the vector potential. This combination is directly 
accessible experimentally, e.g. on a macroscopic scale by studying the trajectory of the 
particle in a magnetic field. Similarly, while the canonical orbital angular momentum 
r x p = —ir x V has a simple interpretation as a property of the charged particle's 
wave function alone, it is unfortunately not gauge invariant. On the other hand, 
the gauge invariant combination r x (p — eA) contains the vector potential A and 
its interpretation is therefore more subtle. Furthermore, only the gauge noninvariant 
orbital angular momentum r xp satisfies the familiar angular momentum commutation 
relations [Li,Lj] = ie 1 ^ L\~. In the end, both definitions of the angular momentum 
of the charged particle have their merit, and depending on the context either may be 
preferable. 

Analogously, in the decomposition of the nucleon spin, it also turns out that 
both schemes have their respective merits. In Ji's decomposition both J q and J G 
are manifestly gauge invariant and calculable in lattice gauge theory (see section 
4.4). Furthermore, using (37) Jf can be measured, which allows determination of 
Lf = Jf — | (Agj + Acjf). However, since it is even more difficult to relate J G to 
measurable quantities (e.g. using (38)) than is the case for Jf, the relation can currently 
not be 'tested', while it can nevertheless be used to determine J G = \ — J q , which is 
interpreted as gluon total angular momentum in the equal time quantized framework 
only. 

Besides a partonic interpretation of all terms, the Jaffe-Manohar decomposition 
(103) has the advantage that it contains the experimentally accessible quantity AG 
(see section 2.8). Unfortunately, also here no experimental observable has been found 
to relate to £/ or C G even in principle, so the light-cone based spin decomposition 
cannot be tested either, but one can still use the sum rule to determine the net 
parton orbital angular momentum C. Another disadvantage of the light-cone based 
decomposition (103) is that because the operators used to define its terms (except AS) 
are not manifestly gauge invariant, its terms cannot be calculated in lattice gauge theory. 
Nevertheless, although a simple (i.e. local) operator expression for AG exists only in 
light-cone gauge, it can be defined through a manifestly gauge invariant operator. In 
all other gauges, AG is described by matrix elements of a non-local operator rendering 
its calculation or interpretation exceedingly difficult. AG has a physical interpretation 
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as the gluon intrinsic angular momentum in the light-cone framework only. 

The two decompositions of the nucleon spin discussed here are not the only 
possibilities. For example, another alternative has been proposed in [248,249]. However, 
it neither has a partonic interpretation nor does its orbital angular momentum seem to 
be experimentally accessible and we will not discuss this decomposition here. 

An analogue may illustrate the dangers in mixing components from different 
decompositions: suppose hypothetically that DIS data on deuterium were available only 
at Q 2 = 100 GeV 2 and proton data were available only at Q 2 = 10 GeV 2 . The normal 
procedure for extracting neutron data would be to evolve to a common Q 2 value first. 
However, again hypothetically, suppose there existed no evolution equations. One may 
then ask oneself whether it would make sense (without first evolving the results from 
one Q 2 to the other value of Q 2 ) to subtract the distribution ofu quarks in the proton 
at Q 2 = 10 GeV 2 from that in the deuteron at Q 2 = 100 GeV 2 and use the results to 
define the distribution of u quarks in the neutron. The main weakness of this analogy 
is of course that this definition is falsifiable because experiments may be performed at 
different values of Q 2 or QCD evolution may be discovered. 

5.2. Implications for quark angular momenta 

5.2.1. Results using the Ji relation Due to the difficulties in extracting GPDs from 
DVCS data, no model-independent experimental constraints on the Ji relation yet exist. 
Lattice-based constraints on Jf have been discussed in section 4.4 and indicate that 
L u + Ld is consistent with zero, while L u — « —0.4 ± 0.05. In that subsection 
it was already emphasized that this result is counterintuitive, where the intuition 
derives from comparison to constituent quark models and phenomenology for single-spin 
asymmetries. Given the complementary and incommensurate definitions of quark orbital 
angular momentum, as discussed in the previous section, the question arises whether 
the orbital angular momentum in these models should be identified with Lf in (102) or 
with Cf in (104). If they were to be identified with Lf, there would be additional, and 
presently unknown, corrections added to account for the difference before the orbital 
angular momentum in the models could be compared to the lattice QCD results. 

5.2.2. Results in the light-cone decomposition The existing information on AG from 
global fits of experimental data suggests that its first moment is small [65]. Neglecting 
'theoretical' and model uncertainties such as those associated with the choice of 
functional forms used to represent PDFs in the fits while interpreting table 2 leads 
to AG ~ —0.1 ± 0.1(expt.) for A% 2 = 1. Using the most precise result for AS = 
0.33 ± 0.03(£ot.) [47] the net parton orbital angular momentum is obtained from (103) 
as: 




(105) 
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which would be more than twice the intrinsic contribution from quarks. It remains open 
how gluons and quarks share their orbital contributions. 

6. Conclusions and Outlook 

One measure of the vitality of a field is the rapidity with which its conceptual framework 
and focus evolve. By this measure, the study of the internal spin structure of the proton 
is indeed vital. Its explosive birth in the late 1980's was triggered by the development of 
experimental technologies for Deeply Inelastic Scattering (DIS) of polarized high energy 
leptons on targets containing polarized protons. The pivotal EMC measurement [17, 18] 
suggested that the helicities of the quarks combine to make a net contribution to the 
helicity of the proton that is small compared to naive expectations. This finding 
inspired a flood of theoretical papers and several major new experimental efforts at 
three laboratories, which confirmed the original finding with much higher precision, and 
extended the measurements to the neutron target to help distinguish the contributions 
of the various quark flavours. While u quarks were found to be strongly polarized in 
the same direction as the parent proton, both d quarks and even s quarks were found 
to be negatively polarized. These measurements of inclusive DIS asymmetries were 
extended to semi-inclusive measurements where the type of hadron produced in the 
fragmentation of the struck quark tends to 'tag' the flavour of that quark. The resulting 
data indicate that the helicity densities of sea quarks are small in the measured x range 
of 0.02 < x < 0.4, but still consistent with fits to inclusive data yielding a negative 
first moment for the strange sea, which receives a significant negative contribution from 
x < 0.02. 

Attention soon focused on a possibly substantial role of the helicity of the gluons, 
which is experimentally difficult to probe in DIS. While some of the existing new 
experiments struggled to extract constraints on the gluon polarization, yet more new 
efforts were launched that were designed specifically to target this problem. At the 
time of writing, these are beginning to bear fruit, suggesting that the fractional gluon 
helicity contribution is small, and together with the quark helicities, do not exhaust the 
'spin budget' of the proton. This leaves the question of the contribution of parton 
orbital angular momentum. In the light-cone framework for defining this quantity 
with a partonic interpretation in the axial gauge, one can estimate it as the difference 
£ = ~ — |AE — AG. Although AE is precisely known from measurements of gf, the 
first moment AG is still so poorly known that there is still little direct experimental 
constraint on the contribution of parton orbital angular momentum to the nucleon spin. 

In the 1990's, dramatic theoretical progress reshaped the field, identifying new 
experimental signals for nucleon spin structure. One key issue has been quark 
transversity, which differs from helicity because of the relativistic internal motion 
of quarks. It was realized that single-spin asymmetries in semi-inclusive scattering 
offered access to this elusive property, through a dependence of the fragmentation 
process on the transverse polarization of the struck quark — the Collins effect. While 
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initial experimental evidence for this phenomenon was ambiguous, eventually SIDIS 
measurements with transverse target polarization in combination with studies of the jet 
structure of hadrons from e + e~ collisions provided a data set with a clear interpretation, 
leading to the first extraction of transversity distributions. As expected, they are found 
to share the features of the helicity distributions in that they are positive for u quarks, 
and negative and considerably smaller for d quarks. Their magnitudes are found be 
be about half of the Soffer positivity bounds, which are expressed in terms of the 
unpolarized and helicity densities. 

One of the ambiguities that complicated the interpretation of early measurements 
of single-spin asymmetries was a possible contribution by the Sivers effect, which 
arises from a 'naively T-odd' correlation between the nucleon spin and the intrinsic 
transverse momentum of quarks. This developed into a fascinating topic of its own 
when it was realized that such distributions raise new questions about key theorems 
for factorization of hard and soft dynamics of processes, and the importance of initial 
or final-state interactions in the understanding of parton distributions that depend 
on intrinsic transverse momentum. The above-mentioned SIDIS measurements with 
transverse target polarization also provided clear evidence of the Sivers effect. Future 
high-luminosity measurements of the Drell-Yan process with transversely polarized 
protons are now anticipated as a test of the new prediction that the Sivers effect will 
appear there with the opposite sign. 

The pursuit of parton orbital angular momentum also led to dramatic theoretical 
progress in the late 1990's. It was realized that experimental access to this quantity 
might be provided by hard exclusive processes that involve hard interactions, yet leave 
the target nucleon intact. Most prominent among these processes is deeply virtual 
Compton scattering, the production of a single energetic photon, because of its fruitful 
interference with the well- understood Bethe-Heitler process, providing access to the 
phase of the DVCS amplitude. Such exclusive processes are now interpreted within the 
framework of generalized parton distributions that depend on three kinematic variables 
and subsume both elastic form factors as moments and ordinary parton distributions 
as special cases. The excitement centered around the finding that a second moment of 
certain GPDs provides the total angular momentum of partons, including orbital angular 
momentum. Although existing DIS experimental facilities were not designed to measure 
exclusive processes, they managed to make some pioneering measurements, which have 
been compared to calculations based on the few available GPD parameterizations. 
While the basic features of the data are described, none of the existing GPD models 
are able to simultaneously explain all available DVCS data, including cross sections 
and their differences or asymmetries with respect to beam charge and helicity, and 
target polarization. Nevertheless, the data for transverse target polarization were 
shown to be sensitive to quark total angular momentum, of which the contribution 
from quark intrinsic spin has already been determined as mentioned above. Hence the 
development of more successful parameterizations can be expected to thereby result in 
constraints on quark orbital angular momentum. Promising efforts already underway 
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based on 'conformal moments' offer a practical means of fitting DVCS data with 
parameterizations at NLO in a s . Some success with this approach has been achieved 
with DVCS cross sections at collider energies, and work continues to extend this to 
charge and polarization dependences at fixed-target energies. 

A decade ago, progress in the understanding of nonperturbative QCD was 
made mainly by comparing experimental data with models that were constructed to 
approximately embody some important features of QCD. Since then, rapid progress 
has been made in improving the scope and precision of numerical simulations of QCD 
on a Euclidean lattice. Improved algorithms and techniques together with rapidly 
increasing computational power have combined to bring this field to the point that 
results for a large variety of observables can now be taken seriously as implications of 
QCD. An example is the spin budget of the proton, where it is computed that the 
only substantial contribution by quarks is from the u flavour, in which the positive 
contribution from quark intrinsic spin dominates the negative contribution from orbital 
angular momentum. While several types of important observables will remain inacessible 
to this approach for the forseeable future (such as the gluonic contribution to the 
spin budget), others that are are calculable are so difficult to measure experimentally 
that lattice calculations will be the only source of information about them for a 
considerable time. Examples here include chiral-odd generalized parton distributions. 
Hence calculations on the lattice play a role that is highly complementary to both 
experiment and model building. 

The existing experimental facilities for polarized DIS have provided a wealth of 
new data, much of it unanticipated when those experiments were conceived. For 
example, they have provided high quality data for inclusive measurements of double- 
spin asymmetries, determining the spin structure function gx of the proton and neutron 
over the kinematic range of fixed-target experiments. However, other important avenues 
to nucleon spin structure such as semi-inclusive DIS and exclusive processes have been 
exploited with only enough precision to inspire and guide rapid theoretical progress, 
but not to provide detailed information. Also, the experience from the HERA e — p 
collider taught us that the vastly wider kinematic range available at such a collider 
can lead to unanticipated insights, in part because of the much improved access to 
the region of small x. One clear motivation for planning for a high-luminosity e — p 
collider with polarized beams is to obtain data for g± over a wide range in hard scale 
Q 2 to allow the precise inference of gluon polarization from the Q 2 evolution of g±, in 
analogy to the success in extracting the unpolarized gluon density from HERA data. 
The rapidly expanding horizons of the field in the last decade have provided other 
compelling motivations. Only with high luminosity free of target dilution by unpolarized 
components and a broad kinematic range can precise information be produced in the 
study of transversity, transverse-momentum-dependent parton distributions and quark 
orbital momentum through the constraint of generalized parton distributions. The clear 
need for a new accelerator facility to make this all possible has inspired specific proposals 
and design studies [250]. Luminosities of order 10 33 nucleons/ (cm 2 s) are considered 
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possible in the energy range y^i = 15-100 GeV. The continuing vitality of this exciting 
field depends on the determined pursuit of this goal. While the community works toward 
this new facility, the present rate of rapid conceptual progress in the field will depend 
on a continuing flow of new data from COMPASS [89], the only remaining high-energy 
lepton DIS facility, as well as from Jefferson Laboratory, especially after it achieves its 
12 GeV upgrade [251]. 
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Table 3. Recent published values of first moments of hclicity densities, separated into experimental evaluations based on the first moment 
of gid, and evaluations from NLO QCD fits of various overlapping data sets. The DSSV values for As + Asarc computed from their fit 
result for As using their assumption As = As. The fit uncertainties for A^ 2 = 1 account only for experimental statistical and systematic 
uncertainties added in quadrature, and not for 'theoretical' or model uncertainties. The DSSV uncertainties apply only to the partial 
moments over the range x > 0.001, and hence should be interpreted with care. All evaluations are in the MS scheme. 



Analysis 




year 


Q 2 

(GeV 2 ) 


AS 


As + As 


AG 


Direct Experimental evaluations 


Hermes [47] (d target) 




2006 


5 


0.330 ± 0.011 (theor.) 
±0.025(exp.) ± 0.028(evol.) 


-0.085 ± 0.013(theor.) 
±0.008(exp.) ± 0.009(evol.) 




Compass [72] 




U 


3 


0.35±0.03(stat.) 
±0.05(syst.) 


-0.08 ± O.Ol(stat) 
±0.02(syst.) 




NLO QCD fits 


AAC [57] (Ax 2 = 12.65) 




2006 


1 


0.27 ±0.07 


-0.10 ±0.02 


0.31 ±0.32 


LSS [53] (Ax 2 = 1) 




2006 


1 


0.17. ..0.31 


-0.068... -0.047 


-0.6... 0.3 


DSSV [65] (Ax 2 = 1) 




2008 


1 


242+ - 015 


-0.114 ±0.02 


-o.o843:Sj 


DSSV [65] (Ax 2 = 0.02, x 2 


= 7.85) 


2008 


1 


n 242+ - 042 
u.z^z_ 062 


-0.114 ±0.06 


-0.0841°;™ 



